[考试反思]1018csp-s模拟测试78(lrd day2) ：规律
作者：互联网
zkt没素质果然考炸了！
但是他考炸了和我一个分
这场的状态是真的不好，T3比较简单但没有做，一直干T2结果还是跪了
T1的哈希写挂了，模数比int大了结果一乘就炸long long了。
调了一个小时，傻逼哈希啊调了一个小时。。。心态初崩。
然后看到T2感觉不是特别难就写出来了。少考虑了一种情况挂成45分。
考后没多久就改过来了。
然后T2弄了一个小时40分钟左右，还剩40分钟看了一下T3的部分分感觉很不好打，不如写好T2，就完了。
注意时间分配。不要硬肝一道题。
不要怂大样例，对于DAG手动DFS一下其实可能很小。
哈希模数不要开到long long级别。
T1：串串香
暴力哈希。
1 #include<cstdio>
2 int max(int a,int b){return a>b?a:b;}
3 #define Hash_Mod 33514539
4 int n,m;char s[1000005];long long hsh[1000005],pw[1000005];
5 int main(){//freopen("ex_ccx2.in","r",stdin);
6 pw[0]=1;
7 for(int i=1;i<=1000000;++i)pw[i]=pw[i-1]*29%Hash_Mod;
8 int t;scanf("%d",&t);
9 while(t--){
10 scanf("%d%d%s",&n,&m,s+1);int max_match=0;
11 for(int i=1;i<=m;++i)hsh[i]=(hsh[i-1]*29+s[i]-'A')%Hash_Mod;
12 for(int i=1;i<m;++i)if(hsh[i]==((hsh[m]-hsh[m-i]*pw[i]%Hash_Mod)+Hash_Mod)%Hash_Mod)
13 if(n==1||hsh[m-i]==((hsh[m]-hsh[i]*pw[m-i]%Hash_Mod)+Hash_Mod)%Hash_Mod)max_match=i;
14 if(n==1&&max_match==0)puts("-1");
15 else printf("%lld\n",(n-1ll)*m+max_match);
16 }
17 }
View Code
T2：迷糊图
代码注释写的很清楚。不懂的在评论区里问就好了。
1 //所谓新建了一条一次性边，其实就是在用到这条边时重新设定了先后手与起点位置，其余不变
2 //分类讨论是否经过了新建边即可
3 #include<cstdio>
4 #include<iostream>
5 using namespace std;
6 int n,m,s,fir[100005],l[300005],to[300005],cnt,out[100005],mx_win_p,mn_win_p;
7 int _fir[100005],_l[300005],_to[300005],_cnt;
8 double win[100005],reach[2][100005],lrd_ans1,lrd_ans2,ans1,ans2;
9 double tot_win,mx_win=-0.1,mn_win=1.1,mx2_win,mn2_win;
10 void link(int a,int b){l[++cnt]=fir[a];fir[a]=cnt;to[cnt]=b;out[b]++;}
11 void _link(int a,int b){_l[++_cnt]=_fir[a];_fir[a]=_cnt;_to[_cnt]=b;}
12 int main(){
13 scanf("%d%d%d",&n,&m,&s);reach[0][s]=1;
14 for(int i=1,x,y;i<=m;++i)scanf("%d%d",&x,&y),link(y,x),_link(x,y);
15 scanf("%lf%lf",&lrd_ans1,&lrd_ans2);
16 //因为题目保证输入的边满足a<b，所以1～n就是一个合法的拓扑序
17 //求在这个点的先手胜率，需要逆向拓扑
18 for(int i=n;i;--i)for(int j=fir[i];j;j=l[j])win[to[j]]+=(1-win[i])/out[to[j]];
19 //求从s开始用奇数/偶数步到达每一个点的概率，正向拓扑
20 for(int i=1;i<=n;++i)for(int j=_fir[i];j;j=_l[j])
21 reach[1][_to[j]]+=reach[0][i]/out[i],reach[0][_to[j]]+=reach[1][i]/out[i];
22 for(int i=1;i<=n;++i){
23 tot_win+=win[i];
24 if(mx_win<=win[i])mx2_win=mx_win,mx_win_p=i,mx_win=win[i];
25 if(mn_win>=win[i])mn2_win=mn_win,mn_win_p=i,mn_win=win[i];
26 }
27 //最大胜率，就是建边把自己引向胜率最大的点或把对方引向胜率最小的点
28 for(int i=1;i<=n;++i){//枚举起点
29 double mx=i==mx_win_p?mx2_win:mx_win,mn=i==mn_win_p?mn2_win:mn_win;//不能是自环，所以如果撞了就要取次大值
30 ans1=max(ans1,max(
31 (win[s]-reach[0][i]*win[i]-reach[1][i]*(1-win[i]))//根本没路过起点的胜率
32 +reach[0][i]/(out[i]+1)*(1-mx)//lrd走了新边
33 +reach[1][i]/(out[i]+1)*mx//对手走了新边
34 +reach[0][i]*out[i]/(out[i]+1)*win[i]//路过起点然后轮到达哥继续走
35 +reach[1][i]*out[i]/(out[i]+1)*(1-win[i])//u路过起点轮到对手继续走
36 ,
37 (win[s]-reach[0][i]*win[i]-reach[1][i]*(1-win[i]))//根本没路过起点的胜率
38 +reach[0][i]/(out[i]+1)*(1-mn)//lrd走了新边
39 +reach[1][i]/(out[i]+1)*mn//对手走了新边
40 +reach[0][i]*out[i]/(out[i]+1)*win[i]//路过起点然后轮到达哥继续走
41 +reach[1][i]*out[i]/(out[i]+1)*(1-win[i])//u路过起点轮到对手继续走
42 ));
43 }
44 //平均胜率，就是建到其它任意点的平均胜率，还是分先后手
45 for(int i=1;i<=n;++i){
46 double tot=(tot_win-win[i])/(n-1);
47 ans2+=(1-(reach[0][i]+reach[1][i])/(out[i]+1))*win[s]//没用新边
48 +reach[0][i]/(out[i]+1)*(1-tot)//lrd走了新边
49 +reach[1][i]/(out[i]+1)*tot;//对手走了新边
50 }
51 printf("%.3lf %.3lf\n",lrd_ans1>-0.5?lrd_ans1:ans1,lrd_ans2>-0.5?lrd_ans2:ans2/n);
52 }
注释版(2.2k)
1 #include<cstdio>
2 #include<iostream>
3 using namespace std;
4 int n,m,s,fir[100005],l[300005],to[300005],c,O[100005],_fir[100005],_l[300005],_to[300005],_c;
5 double w[100005],r[100005],R[100005],A1,A2,Tot,Mx,Mn=1,p,tot;
6 void L(int a,int b){l[++c]=fir[a];fir[a]=c;to[c]=b;O[b]++;}
7 void _L(int a,int b){_l[++_c]=_fir[a];_fir[a]=_c;_to[_c]=b;}
8 int main(){
9 scanf("%d%d%d",&n,&m,&s);r[s]=1;
10 for(int i=1,x,y;i<=m;++i)scanf("%d%d",&x,&y),L(y,x),_L(x,y);
11 for(int i=n;i;--i)for(int j=fir[i];j;j=l[j])w[to[j]]+=(1-w[i])/O[to[j]];
12 for(int i=1;Tot+=w[i],Mx=max(Mx,w[i]),Mn=min(Mn,w[i]),i<=n;++i)for(int j=_fir[i];j;j=_l[j])R[_to[j]]+=r[i]/O[i],r[_to[j]]+=R[i]/O[i];
13 for(int i=1;p=r[i]<R[i]?Mx:Mn,i<=n;++i)A1=max(A1,((w[s]-r[i]*w[i]-R[i]*(1-w[i]))+r[i]/(O[i]+1)*((1-p)+O[i]*w[i])+R[i]/(O[i]+1)*(p+O[i]*(1-w[i]))));
14 for(int i=1;tot=(Tot-w[i])/(n-1),i<=n;++i)A2+=(1-(r[i]+R[i])/(O[i]+1))*w[s]+r[i]/(O[i]+1)*(1-tot)+R[i]/(O[i]+1)*tot;
15 printf("%.3lf %.3lf\n",A1,A2/n);
16 }
究级压码长版(969B,但其实单行不是很长所以并不算太压行)
样例3简化图：
T3：木叶下
维护子树内的贡献，父亲方向上的贡献，对于每个点p其父亲f除却p以外的儿子的贡献。
树上倍增，还需要维护子树最大，第二大，第三大。
1 #include<cstdio>
2 int max(int a,int b){return a>b?a:b;}
3 int fi[200005],se[200005],th[200005],fir[200005],l[400005],to[400005],son[200005],ANS;
4 int n,m,cnt,fip[200005],sep[200005],fav[200005][20],f[200005][20],dep[200005],F[200005];
5 void link(int a,int b){l[++cnt]=fir[a];fir[a]=cnt;to[cnt]=b;}
6 void dfs(int p,int fa){
7 f[p][0]=fa;dep[p]=dep[fa]+1;
8 for(int i=1;i<=18;++i)f[p][i]=f[f[p][i-1]][i-1];
9 for(int i=fir[p];i;i=l[i])if(to[i]!=fa){
10 dfs(to[i],p);int mx=son[to[i]]+1;
11 if(mx>fi[p])th[p]=se[p],se[p]=fi[p],sep[p]=fip[p],fi[p]=mx,fip[p]=to[i];
12 else if(mx>se[p])th[p]=se[p],se[p]=mx,sep[p]=to[i];
13 else if(mx>th[p])th[p]=mx;
14 if(son[p]<mx)son[p]=mx;
15 }
16 for(int i=fir[p];i;i=l[i])if(to[i]!=fa){
17 if(fip[p]!=to[i])fav[to[i]][0]=fi[p];
18 else if(sep[p]!=to[i])fav[to[i]][0]=se[p];
19 else fav[to[i]][0]=th[p];
20 }ANS=max(ANS,fi[p]+se[p]);
21 }
22 void DFS(int p,int fa){
23 for(int i=1;i<=18;++i)fav[p][i]=max(fav[p][i-1],fav[f[p][i-1]][i-1]);
24 for(int i=fir[p];i;i=l[i])if(to[i]!=fa){
25 F[to[i]]=F[p];
26 if(fip[p]==to[i])F[to[i]]=max(F[to[i]],se[p]);
27 else F[to[i]]=max(F[to[i]],fi[p]);
28 F[to[i]]++;
29 DFS(to[i],p);
30 }
31
32 }
33 int main(){//freopen("ex_myx2.in","r",stdin);freopen("my.out","w",stdout);
34 //freopen("r.out","r",stdin);freopen("t3.out","w",stdout);
35 scanf("%d",&n);
36 for(int i=1,x,y;i<n;++i)scanf("%d%d",&x,&y),link(x,y),link(y,x);
37 dfs(1,0);DFS(1,0);
38 scanf("%d",&m);int a,b,ans,sub,lca,ra,rb;
39 while(m--){
40 scanf("%d%d",&a,&b);ans=0;//printf("%d %d\n",a,b);
41 if(a==b){printf("%d\n",(ANS>>1)+1);continue;}
42 sub=dep[a]-dep[b];if(sub<0)a^=b^=a^=b,sub=-sub;ra=a;rb=b;
43 for(int i=18;~i;--i)if(sub&1<<i)a=f[a][i];
44 if(a==b)lca=a,ans=son[ra];
45 else{
46 for(int i=18;~i;--i)if(f[a][i]!=f[b][i])a=f[a][i],b=f[b][i];
47 lca=f[a][0];ans=max(son[ra],son[rb]);
48 }a=ra;b=rb;
49 sub=dep[a]-dep[lca]-1;
50 for(int i=18;~i;--i)if(sub&1<<i)ans=max(ans,fav[a][i]),a=f[a][i];
51 sub=dep[b]-dep[lca]-1;
52 if(sub>0)for(int i=18;~i;--i)if(sub&1<<i)ans=max(ans,fav[b][i]),b=f[b][i];
53 if(!(fip[lca]==a||fip[lca]==b))ans=max(ans,fi[lca]);
54 else if(!(sep[lca]==a||sep[lca]==b))ans=max(ans,se[lca]);
55 else ans=max(ans,th[lca]);
56 ans=max(ans,F[lca]);
57 printf("%d\n",ans);
58 }//printf("%d\n",fav[6][0]);
59 }
View Code
标签：fir,cnt,lrd,int,1018csp,day2,100005,win,200005 来源： https://www.cnblogs.com/hzoi-DeepinC/p/11699730.html