[考试反思]0306省选模拟39：延续
俩原题。对于看视频课太匆忙的我来说，还是算俩新题吧。
为啥感觉最近总是考第$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 }
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 }
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 }
没有调到暴毙还是非常开心的，写的时间比调的时间长。
