码迷,mamicode.com
首页 > 其他好文 > 详细

[考试反思]0306省选模拟39:延续

时间:2020-03-07 09:32:38      阅读:69      评论:0      收藏:0      [点我收藏+]

标签:std   三元组   space   形式   展开   之间   次数   class   img   

技术图片

 技术图片

俩原题。对于看视频课太匆忙的我来说,还是算俩新题吧。

 

 为啥感觉最近总是考第$9$啊。最近这种名次的频率好高。。

(好菜啊第$9$还学个啥呢啊脑子呢???这还翻个什么盘啊)

然而其实第$9$也是混来的。

最开始先看的$T1$发现好像是原题,但是印象不深,大概还记得第一步怎么做。

接着往后看,$T2$看起来像弱智题。$T3$的话又是原题,思路极其简单但是是$exSAM+LCT$写出来会暴毙。于是扔在最后。

然后回$T1$尝试回忆。然后就没有然后了。折腾来折腾去想了两个小时硬是没想到$dp$定义。

果断换题。然后一看$T2$,好像的确有点弱智啊。

思路大概是对的,但是在枚举三元环上卡了一下。然后就真的忘记了怎么$O(m \sqrt{m})$枚举三元环了。

然后的确就一直没有想起来。又耗了一个小时,最后利用大小点卡常数,理论上时间复杂度还是$O(n^2)$

但是实际上常数非常小,我自己貌似也构造出什么数据卡满我自己。然而时间也不多了,调着调着就交了。

最后剩了十几分钟给$T3$。尝试爆写$exSAM$最后还是放弃了。

我为啥又把暴跳父亲忘了???如果想起来的话这次就能有两道题都是$n^2$过$2 \times 10^5$了!

然而既然考场上没写。。。那考后就老老实实写正解吧。。。

强烈谴责RNB考后暴力暴跳父亲AC的无脸行为

话说我自己数据结构能力是真的弱啊

 

T1:gift

大意:两数列$a,b$。有些位置为空,问有多少种方案填数使之成为排列,满足通过至少$k$次交换$a$变成$b$。对所有$k \in [0,n-1]$求解。$n \le 2000$

如果序列已经填满了数,那么最少交换次数可以这么表示:所有$a_i \rightarrow b_i$连边。然后每个环需要交换$size-1$次。所以交换次数就是$n-$环个数。

现在有些位置为$0$了,那么所有的边分四类$0-x,x-y,x-0,0-0$。而有些边已经有公共点了($0$除外)。把它们缩成一条链

对于$x-y$这一类,我们发现在缩成链之后两端点都未消失,证明它们都只出现了一次。

我们可以把它们两个当成同一个数看待,可以发现这对最后的图没有影响,在最后的时候展开就行了。

所以边只剩下三类,我们分别讨论它们之间的关系:

我们发现,$0-x$这种边它的$x$一定只出现了一次,所以一定与$0$相连了,那么要么就是接上了一个$0-0$,要么就是接上了另一个$0-x$,或者说自己连上自己变成自环。

对于第一种情况,新的两个端点就是$0-0$了,那么相当于这条边变成$0-0$了,而同时$0-0$的总量加一减一后并没有改变。

对于第二种情况,那么就是合并成了一个新的$0-x$。对最后的环数答案并没有什么影响。

所以我们只在意$0-x$这一类边自己形成了多少个环,以及方案数。

直接算不好算,考虑容斥,得到$i$个环的方案数是$g_i =\sum\limits_{j=i}^{cnt_{0-x}} \binom{cnt_{0-x}}{j} \begin{bmatrix} j \\ i \end{bmatrix} (cnt_{0-0}+cnt_{0-x}-j)^{\underline{cnt_{0-x}-j}}$

意思就是说,枚举有多少条$0-x$边在环里,对于剩下的边,我们在$cnt_{0-0}$个$0-0$边和$cnt_{0-x}$个$0-x$边中任意选择一条接在后面。方案数是下降幂。

后面乱选的部分也有可能形成环,所以这个$g$是至少的形式。用一个简单的二项式反演就能得到恰好的形式。

对于所有的$x-0$边,同理。它和$0-x$边如果要产生关系,一定要借助$0-0$而上面的第一种情况说了只要用到$0-0$就会变成$0-0$所以只需要考虑$x-0$与自身以及$0-0$的关系。

这两部分像个背包。直接$OGF$卷起来就好了。得到的是把所有的$0-x$边以及$x-0$边连成环或者变味$0-0$的方案数了。

所以现在只需要继续考虑$cnt_{0-0}$条$0-0$边就好了。

这个的方案数比较好算,环的形状一共是第二类斯特林数种,然后你再给它们赋予$cnt_{0,0}$的先后顺序,也就是个阶乘。

一起卷起来就得到答案了。

$Dy$讲的时候真实啥也没听懂,自己回放的时候其实也不明白,说啥都没有自己做一遍有效。

技术图片
 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define mod 998244353
 4 int mo(int a){return a>=mod?a-mod:a;}
 5 #define S 2002
 6 int n,a[S],b[S],al[S],c00,c01,c10,cir,in[S],out[S],C[S][S],s[S][S],A[S][S],G[S],g[S],H[S],h[S],f[S],ans[S];
 7 int main(){
 8     scanf("%d",&n);
 9     for(int i=1;i<=n;++i)scanf("%d",&a[i]);
10     for(int i=1;i<=n;++i)scanf("%d",&b[i]);
11     for(int i=1;i<=n;++i)in[i]=out[i]=-1;
12     for(int i=1;i<=n;++i)out[a[i]]=b[i],in[b[i]]=a[i];
13     for(int i=1;i<=n;++i)if(in[i]==-1&&out[i]!=-1){
14         int p=i;al[p]=1;
15         while(p&&out[p]!=-1)p=out[p],al[p]=1;
16         if(!p)c10++;
17     }
18     for(int i=1;i<=n;++i)if(in[i]==0){
19         int p=i;al[p]=1;
20         while(p&&out[p]!=-1)p=out[p],al[p]=1;
21         if(!p)c00++;else c01++;
22     }
23     for(int i=1;i<=n;++i)c00+=a[i]==0&&b[i]==0;
24     for(int i=1;i<=n;++i)if(in[i]>0&&out[i]>0&&!al[i]){
25         int p=i;cir++;
26         while(!al[p])al[p]=1,p=out[p];
27     }
28     for(int i=0;i<=n;++i)for(int j=C[i][0]=1;j<=i;++j)C[i][j]=mo(C[i-1][j-1]+C[i-1][j]);
29     for(int i=s[0][0]=1;i<=n;++i)for(int j=1;j<=i;++j)s[i][j]=(s[i-1][j-1]+s[i-1][j]*(i-1ll))%mod;
30     for(int i=0;i<=n;++i)for(int j=A[i][0]=1;j<=i;++j)A[i][j]=A[i][j-1]*(i+1ll-j)%mod;
31     for(int i=0;i<=c01;++i)for(int j=i;j<=c01;++j)g[i]=(g[i]+1ll*C[c01][j]*s[j][i]%mod*A[c01+c00-j][c01-j])%mod;
32     for(int i=0;i<=c01;++i)for(int j=i;j<=c01;++j)G[i]=(G[i]+(j-i&1?mod-1ll:1ll)*g[j]%mod*C[j][i])%mod;
33     for(int i=0;i<=c10;++i)for(int j=i;j<=c10;++j)h[i]=(h[i]+1ll*C[c10][j]*s[j][i]%mod*A[c10+c00-j][c10-j])%mod;
34     for(int i=0;i<=c10;++i)for(int j=i;j<=c10;++j)H[i]=(H[i]+(j-i&1?mod-1ll:1ll)*h[j]%mod*C[j][i])%mod;
35     for(int i=0;i<=c01;++i)for(int j=0;j<=c10;++j)f[i+j]=(f[i+j]+1ll*G[i]*H[j])%mod;
36     for(int i=0;i<=c01+c10;++i)for(int j=0;j<=c00;++j)ans[i+j]=(ans[i+j]+1ll*f[i]*s[c00][j]%mod*A[c00][c00])%mod;
37     for(int i=0;i<n;++i)printf("%d ",n>=i+cir?ans[n-i-cir]:0);
38 }
View Code

 

T2:girls

大意:$[0,n)$。定义合法三元组为三者之间均无连边且$i<j<k$。贡献为$Ai+Bj+Ck$。求总贡献。$n,m \le 2\times 10^5$

这题也就是经典的容斥了,全部情况加起来减去至少有一条边加上至少有两条边再去掉三元环

别的都好说,三元环的问题大概就是按照度数排序那一套复杂度$O(m^{1.5})$

做法也就是把所有的点按照度数排序之后,度数小的向度数大的连边。

对于点$i$枚举所有出边打标记,再枚举一次出边$j$,再枚举$j$的儿子$k$。若$k$有标记则是三元环。

联赛前做过,现在忘掉了,每个三元环只在度数最小的点初被枚举,具体复杂度证明忘记了,感性理解一下吧。

考场上思维僵化只会严格按照$\sqrt{n}$为阈值分大小点了,然而一小点两大点的情况不能很好地统计,于是又退化成$O(n^2)$了。

但是不难发现这样做的常数很小很小。所以$O(n^2)$过$200000$没什么大问题。

技术图片
 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define S 200005
 4 #define ull unsigned long long
 5 ull A,B,C,ans,tot[S],deg[S]; int n,m;
 6 vector<int>v[S],in[S],vs[S],vb[S],inb[S],Bi;
 7 unordered_map<int,bool>M[S];
 8 ull cal(int x){return x*(x-1ull)>>1;}
 9 int main(){
10     cin>>n>>m>>A>>B>>C; int sq=sqrt(m); 
11     for(int i=0,a,b;i<m;++i){
12         scanf("%d%d",&a,&b);
13         if(a>b)swap(a,b); deg[a]++; deg[b]++; M[a][b]=M[b][a]=1;
14         v[a].push_back(b); in[b].push_back(a); tot[a]+=b;
15         ans-=cal(a)*A+a*(a*B+b*C) + (cal(b)-cal(a+1))*B+(b-a-1)*(a*A+b*C) + (cal(n)-cal(b+1))*C+(n-b-1)*(a*A+b*B);
16     }
17     for(int i=0;i<n;++i)sort(v[i].begin(),v[i].end()),sort(in[i].begin(),in[i].end());
18     for(int i=0;i<n;++i)ans+=i*(cal(n-1-i)*A+B*i*(n-1-i)+C*cal(i));
19     for(int i=0;i<n;++i){
20         ull x=0,cnt=0;
21         for(auto j:v[i])ans+=i*cnt*A+x*B+j*cnt*C,x+=j,cnt++;
22     }
23     for(int i=0;i<n;++i){
24         ull x=0,cnt=0;
25         for(auto j:in[i])ans+=i*cnt*C+x*A+j*cnt*B,x+=j,cnt++;
26     }
27     for(int i=0;i<n;++i)for(auto j:v[i])ans+=v[j].size()*(i*A+j*B)+tot[j]*C;
28     
29     for(int i=0;i<n;++i)for(auto j:v[i])if(deg[j]>=sq)vb[i].push_back(j);else vs[i].push_back(j);
30     for(int i=0;i<n;++i)for(auto j:in[i])if(deg[j]>=sq)inb[i].push_back(j);
31     for(int i=0;i<n;++i)if(deg[i]>=sq&&vb[i].size()>1)for(auto j:vb[i])for(auto k:vb[j])if(M[i][j]&&M[j][k]&&M[i][k])ans-=i*A+j*B+k*C;
32     for(int i=0;i<n;++i)if(deg[i]<sq&&vs[i].size()>1)for(auto j:vs[i])for(auto k:vs[j])if(M[i][j]&&M[j][k]&&M[i][k])ans-=i*A+j*B+k*C;
33     for(int i=0;i<n;++i)if(deg[i]<sq&&vs[i].size()>1)for(auto j:vs[i]){
34         for(auto k:vb[i])if(M[i][j]&&M[j][k]&&M[i][k])ans-=i*A+min(j,k)*B+max(j,k)*C;
35         for(auto k:inb[i])if(M[i][j]&&M[j][k]&&M[i][k])ans-=k*A+i*B+j*C;
36     }
37     
38     for(int i=0;i<n;++i)if(deg[i]<sq){
39         for(auto j:inb[i])Bi.push_back(j);
40         for(auto j:vb[i])Bi.push_back(j);
41         for(int j=0;j<Bi.size();++j)for(int k=j+1;k<Bi.size();++k){
42             ull x=i,y=Bi[j],z=Bi[k];
43             if(x>z)swap(x,z);
44             if(x>y)swap(x,y);
45             if(M[x][y]&&M[x][z]&&M[y][z])ans-=A*x+B*y+C*z;
46         }Bi.clear();
47     }
48     cout<<ans<<endl;
49 }
View Code

 

T3:string

原题链接。题号12.

思路很简单,代码很难写。

干就是了。

技术图片

技术图片
 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define S 888888
 4 int n,rt=1,pc=1,lst[22],c[11][S],len[S],ans,m,la,f[S],t,w[22][S],lz[22][S];char s[20000003];
 5 vector<int>v[22],P[22];
 6 struct LCT{
 7     int c[2][S],f[S],s[S];
 8     void down(int p){for(int i=1;i<=n;++i)if(lz[i][p])w[i][c[0][p]]+=lz[i][p],w[i][c[1][p]]+=lz[i][p],lz[i][c[0][p]]+=lz[i][p],lz[i][c[1][p]]+=lz[i][p],lz[i][p]=0;} 
 9     bool nr(int p){return c[0][f[p]]==p||c[1][f[p]]==p;}
10     void spin(int p){
11         int F=f[p],G=f[F],d=c[1][F]==p,B=c[!d][p];
12         if(nr(F))c[c[1][G]==F][G]=p; c[!d][p]=F; c[d][F]=B;
13         f[f[f[B]=F]=p]=G;
14     }
15     void splay(int p){
16         int F=p,t=0; while(nr(F))s[++t]=F,F=f[F]; s[++t]=F; while(t)down(s[t--]);
17         for(;F=f[p],nr(p);spin(p))if(nr(F))spin(c[1][f[F]]==F^c[1][F]==p?p:F);
18     }
19     void access(int p){for(int r=0,s=p;s;s=f[r=s])splay(s),c[1][s]=r;splay(p);}
20     void link(int a,int b){splay(a);f[a]=b;}
21     void cut(int p){access(p);f[c[0][p]]=0;c[0][p]=0;}
22 }T;
23 void add(int p,int op){T.access(p);for(int i=1;i<=n;++i)lz[i][p]+=op*w[i][p];}
24 void Fa(int p,int fa){T.link(p,fa);add(p,1);f[p]=fa;}
25 void extend(int I,int C){
26     int p=lst[I],np,q,nq;
27     if(q=c[C][p]){
28         if(len[q]==len[p]+1){lst[I]=q;T.access(q);lz[I][q]++;w[I][q]++;return;}
29         nq=++pc; len[nq]=len[p]+1; for(int i=0;i<10;++i)c[i][nq]=c[i][q];
30         w[I][nq]++;
31         Fa(nq,f[q]); add(q,-1); T.cut(q); Fa(q,nq);
32         for(;c[C][p]==q;p=f[p])c[C][p]=nq; lst[I]=nq;
33     }else{
34         lst[I]=np=++pc; len[np]=len[p]+1; w[I][np]=1;
35         for(;p&&!c[C][p];p=f[p])c[C][p]=np;
36         if(!p){Fa(np,1);ans+=len[np]-len[f[np]];return;}
37         if(len[q=c[C][p]]==len[p]+1){Fa(np,q);ans+=len[np]-len[f[np]];return;}
38         nq=++pc; len[nq]=len[p]+1; for(int i=0;i<10;++i)c[i][nq]=c[i][q];
39         Fa(nq,f[q]); Fa(np,nq); add(q,-1); T.cut(q); Fa(q,nq);
40         for(;c[C][p]==q;p=f[p])c[C][p]=nq; ans+=len[np]-len[f[np]];
41     }
42 }
43 int main(){//freopen("string1.in","r",stdin);
44     scanf("%d%d",&n,&t);
45     for(int i=1;i<=n;++i){
46         scanf("%s",s);lst[i]=1;
47         for(int j=0;s[j];++j)extend(i,s[j]-48);
48         v[i].push_back(0);P[i].push_back(lst[i]);
49     }scanf("%d",&m);
50     for(int i=1;i<=m;++i){
51         int op,x,y,z;scanf("%d",&op);
52         if(op==1){
53             scanf("%d%d",&x,&y);
54             if(t)y^=la,y%=10;
55             extend(x,y);
56             v[x].push_back(i);P[x].push_back(lst[x]);
57         }if(op==2){
58             scanf("%d%d%d",&x,&y,&z);
59             int p=P[x][upper_bound(v[x].begin(),v[x].end(),y)-1-v[x].begin()];
60             T.splay(p);printf("%d\n",la=w[z][p]);
61         }if(op==3)printf("%d\n",ans);
62         if(op==4){
63             int p=1;la=0;scanf("%s",s);
64             for(int j=0;s[j];++j)p=c[s[j]-48][p];
65             if(p){T.splay(p);for(int j=1;j<=n;++j)la=max(la,w[j][p]);} 
66             printf("%d\n",la);
67         }
68     }
69 }
View Code

没有调到暴毙还是非常开心的,写的时间比调的时间长。

 

[考试反思]0306省选模拟39:延续

标签:std   三元组   space   形式   展开   之间   次数   class   img   

原文地址:https://www.cnblogs.com/hzoi-DeepinC/p/12431632.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!