AtCoder Grand Contest
一句话题解
QwQ主要是因为这篇文章写的有点长……有时候要找某一个题可能不是很好找，所以写了这个东西。
具体的题意、题解和代码可以再往下翻。_(:з」∠)_
AGC 001
C：枚举中点/中边。
D：构造。
E：根据组合数意义转为$DP$。
F：拓扑排序，线段树优化连边。
AGC 002
C：水题，看是否有a[i]+a[i+1]>=k。
D：并查集上倍增，二分答案。
E：博弈（坑）
F：模型转化然后$DP$。
AGC 003
C：一个数到自己应到位置距离为奇数的个数/2。
D：数学，质因数分解，素数筛。
E：转化为递增序列，拆分（建议看代码理解……）
F：难的一比的矩阵快速幂。
AGC 004
C：构造。
D：贪心，满了就断。
E：$DP$，设$f[i][j][k][l]$表示取完了$(i,j)(k,l)$这个矩形内的机器人。
F：分类讨论，模型转化。
AGC 005
C：构造。
D：建立二分图，转化后$DP$。
E：贪心？博弈？
F：$NTT$。
AGC 006
C：差分，倍增。
D：二分塔顶的数，把底层$01$转化。
E：结论题……
F：三元环，对染色过程分类讨论。
AGC 007
C：期望。不会。
D：线段树优化$DP$。
E：二分，合并二元组。
F：用折线覆盖变换的过程，需要操作的步数就是折线的个数。
AGC 008
C：讨论奇偶性。
D：按出现位置$sort$一下从前往后填就好了。
E：没看懂题解。（坑）
F：$DP$求最短路和次短路。
AGC 009
C：设f[i]表示第i个元素放到Y集合的方案数。
D：神仙性质题？
E：神仙$DP$。
AGC 010
C：贪心。
D：博弈，根据偶数个数的奇偶来讨论一下。
AGC 001
C - Shorten Diameter
题意：
给定一颗节点个数为n的树，问最少删除多少个节点能使得直径小于等于k
题解：
若k为偶数，枚举作为直径中点的点。若k为奇数，枚举作为直径中点的边。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<queue> 5 #define N (2000+100) 6 using namespace std; 7 8 struct Edge{int to,next;}edge[N<<1]; 9 int n,k,u[N],v[N],dis[N],ans=2001; 10 int head[N],num_edge; 11 queue<int>q; 12 13 void add(int u,int v) 14 { 15 edge[++num_edge].to=v; 16 edge[num_edge].next=head[u]; 17 head[u]=num_edge; 18 } 19 20 int main() 21 { 22 scanf("%d%d",&n,&k); 23 for (int i=1; i<=n-1; ++i) 24 { 25 scanf("%d%d",&u[i],&v[i]); 26 add(u[i],v[i]); add(v[i],u[i]); 27 } 28 29 if (k%2==0) 30 { 31 for (int i=1; i<=n; ++i) 32 { 33 int cnt=1; 34 memset(dis,-1,sizeof(dis)); 35 dis[i]=0; q.push(i); 36 while (!q.empty()) 37 { 38 int x=q.front(); q.pop(); 39 for (int j=head[x]; j; j=edge[j].next) 40 if (dis[edge[j].to]==-1) 41 { 42 dis[edge[j].to]=dis[x]+1; 43 if (dis[edge[j].to]<=k/2) cnt++; 44 q.push(edge[j].to); 45 } 46 } 47 ans=min(ans,n-cnt); 48 } 49 } 50 else 51 { 52 for (int i=1; i<=n-1; ++i) 53 { 54 int cnt=2; 55 memset(dis,-1,sizeof(dis)); 56 dis[u[i]]=dis[v[i]]=0; q.push(u[i]); q.push(v[i]); 57 while (!q.empty()) 58 { 59 int x=q.front(); q.pop(); 60 for (int j=head[x]; j; j=edge[j].next) 61 if (dis[edge[j].to]==-1) 62 { 63 dis[edge[j].to]=dis[x]+1; 64 if (dis[edge[j].to]<=k/2) cnt++; 65 q.push(edge[j].to); 66 } 67 } 68 ans=min(ans,n-cnt); 69 } 70 } 71 printf("%d",ans); 72 }
D - Arrays and Palindrome
题意：
给出一个数字之和为N的顺序可变的数组，构造一个数字和为N的数组，满足对于任意满足前两个条件的字符串，一定满足第三个条件:
1.从头开始的a1个字符，a1个后紧接着的a2个、a3个……分别为回文串。
2.从头开始的b1个字符，b1个后紧接着的b2个、b3个……分别为回文串。
3.该字符串仅由同一字符构成。
题解：
一个鬼畜有趣的构造题……还TM被卡输出空格和回车了……→QAQ
1 #include<iostream> 2 #include<cstdio> 3 #define N (100000+1000) 4 using namespace std; 5 6 int a[N],b[N],cnt,n,m; 7 8 int main() 9 { 10 scanf("%d%d",&n,&m); 11 for (int i=1; i<=m; ++i) 12 { 13 scanf("%d",&a[i]); 14 if (a[i]%2) b[++cnt]=i; 15 } 16 if (cnt>2) {puts("Impossible"); return 0;} 17 if (m==1) 18 { 19 if (a[1]==1) printf("1\n1\n1\n"); 20 else printf("%d\n2\n1 %d\n",a[1],a[1]-1); 21 return 0; 22 } 23 if (b[1]) swap(a[1],a[b[1]]); 24 if (b[2]) swap(a[m],a[b[2]]); 25 26 for (int i=1; i<=m; ++i) 27 printf("%d%c",a[i]," \n"[i==m]); 28 a[1]++; a[m]--; 29 printf("%d\n",a[m]==0?m-1:m); 30 for (int i=1; i<=m; ++i) 31 if (a[i]!=0) printf("%d%c",a[i]," \n"[i==m]); 32 }
E - BBQ Hard
题意：
询问$\sum_{1 \leq i,j \leq n,i \neq j} {a_i+a_j \choose a_i+a_j+b_i+b_j}$
题解：
一个比较巧妙的题，对于任意i,j，可以发现${a_i+a_j \choose a_i+a_j+b_i+b_j}$可以被描述成平面上从$(-a_i,-b_i)$到$(a_j,b_j)$的方案数
所以只需要把平面上所有点$(-a_i,-b_i)$标记成1，然后暴力dp就可以了：$f_{i,j}=f_{i-1,j}+f_{i,j-1}$，注意要把从$(-a_i,-b_i)$到$(a_i,b_i)$的方案数减去
1 #include<iostream> 2 #include<cstdio> 3 #define N (200000+100) 4 #define MOD (1000000007) 5 using namespace std; 6 long long n,a[N],b[N],f[5005][5005],ans; 7 long long Inv[N],Fac[N],FacInv[N]; 8 9 void Init() 10 { 11 Fac[0]=1; Inv[1]=1; FacInv[0]=1; 12 for (int i=1; i<=10000; ++i) 13 { 14 if (i!=1) Inv[i]=(MOD-MOD/i)*Inv[MOD%i]%MOD; 15 Fac[i]=Fac[i-1]*i%MOD; FacInv[i]=FacInv[i-1]*Inv[i]%MOD; 16 } 17 } 18 19 long long C(long long n,long long m) 20 { 21 if (m>n) return 0; 22 return Fac[n]*FacInv[m]%MOD*FacInv[n-m]%MOD; 23 } 24 25 int main() 26 { 27 Init(); 28 scanf("%lld",&n); 29 for (int i=1; i<=n; ++i) 30 { 31 scanf("%lld%lld",&a[i],&b[i]); 32 f[-a[i]+2000][-b[i]+2000]++; 33 } 34 for (int i=0; i<=4000; ++i) 35 for (int j=0; j<=4000; ++j) 36 { 37 if (i) (f[i][j]+=f[i-1][j])%=MOD; 38 if (j) (f[i][j]+=f[i][j-1])%=MOD; 39 } 40 for (int i=1; i<=n; ++i) 41 (ans+=f[a[i]+2000][b[i]+2000])%=MOD; 42 for (int i=1; i<=n; ++i) 43 ans=((ans-C(2*a[i]+2*b[i],2*b[i]))%MOD+MOD)%MOD; 44 printf("%lld",(ans*Inv[2])%MOD); 45 }
F - Wide Swap
题意：
给你一个序列，你可以任意交换两个距离大于等于k，且差的绝对值为1的数的位置，求可以实现的字典序最小的序列。
题解：
一个神题……首先我们设$pos_{a_i}=i$，显然要令a字典序最小，只要让pos字典序最小就好了。因为“权值小的尽量靠前”和“前面的权值尽量小”是一个意思。
那么对于pos数组，相邻且权值差>=k的可以交换顺序。也就是说，i与后面所有的j比，若$abs(pos_i-pos_j)<k$，那么$pos_i$和$pos_j$的相对顺序就是确定不变的，连一条$pos_i$到$pos_j$的边即可。这显然是一个DAG，连完边跑拓扑排序。
但这样的话显然边数是$n^2$级别的，发现有一些边是没有用的，所以考虑优化掉一些的边。
例如三个点x,y,z。x->y, y->z, x->z，可以发现x->z这条边是没有用的。
对于$pos_i$，我们是向$(pos_i-k+1,pos_i-1)$，$(pos_i+1,pos_i+k-1)$这两个区间连边的。我们可以从n~1枚举，查询$(1,pos_i-1)$和$(pos_i+1,n)$两个区间内的两个最小元素x和y，$pos_i$向$pos_x$和$pos_y$连边即可(可行的原因可以考虑一下拓扑序小的在前的性质)。查询最小值可以用线段树。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<queue> 5 #define N (500000+1000) 6 #define INF (10000000) 7 using namespace std; 8 9 struct Edge{int to,next;}edge[N<<2]; 10 int n,k,a[N],pos[N],Segt[N<<2]; 11 int head[N],num_edge,Ind[N]; 12 priority_queue<int,vector<int>,greater<int> >q; 13 14 void add(int u,int v) 15 { 16 edge[++num_edge].to=v; 17 edge[num_edge].next=head[u]; 18 head[u]=num_edge; 19 Ind[v]++; 20 } 21 22 void Build(int now,int l,int r) 23 { 24 if (l==r){Segt[now]=INF; return;} 25 int mid=(l+r)>>1; 26 Build(now<<1,l,mid); Build(now<<1|1,mid+1,r); 27 Segt[now]=min(Segt[now<<1],Segt[now<<1|1]); 28 } 29 30 int Query(int now,int l,int r,int l1,int r1) 31 { 32 if (r<l1 || l>r1) return INF; 33 if (l1<=l && r<=r1) return Segt[now]; 34 int mid=(l+r)>>1; 35 return min(Query(now<<1,l,mid,l1,r1),Query(now<<1|1,mid+1,r,l1,r1)); 36 } 37 38 void Update(int now,int l,int r,int x,int k) 39 { 40 if (l==r) {Segt[now]=k; return;} 41 int mid=(l+r)>>1; 42 if (x<=mid) Update(now<<1,l,mid,x,k); 43 else Update(now<<1|1,mid+1,r,x,k); 44 Segt[now]=min(Segt[now<<1],Segt[now<<1|1]); 45 } 46 47 void Toposort() 48 { 49 int cnt=0; 50 for (int i=1; i<=n; ++i) if (!Ind[i]) q.push(i); 51 while (!q.empty()) 52 { 53 int x=q.top(); q.pop(); a[x]=++cnt; 54 for (int i=head[x]; i; i=edge[i].next) 55 { 56 Ind[edge[i].to]--; 57 if (!Ind[edge[i].to]) q.push(edge[i].to); 58 } 59 } 60 for(int i=1; i<=n; ++i) printf("%d\n",a[i]); 61 } 62 63 int main() 64 { 65 scanf("%d%d",&n,&k); 66 for (int i=1; i<=n; ++i) 67 scanf("%d",&a[i]),pos[a[i]]=i;; 68 Build(1,1,n); 69 for (int i=n; i>=1; --i) 70 { 71 int x=Query(1,1,n,pos[i]+1,min(pos[i]+k-1,n)); 72 if (x!=INF) add(pos[i],pos[x]); 73 int y=Query(1,1,n,max(pos[i]-k+1,1),pos[i]-1); 74 if (y!=INF) add(pos[i],pos[y]); 75 Update(1,1,n,pos[i],i); 76 } 77 Toposort(); 78 }
AGC 002
C - Knot Puzzle
题意：
有n段绳子，将他们打结形成一段有n-1个绳结的绳子。每次可以选择一段长度大于k的绳子解开它的任意一个绳结。问将所有绳结解开的顺序。
题解：
水题，看是否有a[i]+a[i+1]>=k，有的话以这两个为中心两头拆。
1 #include<iostream> 2 #include<cstdio> 3 #define N (100000+1000) 4 using namespace std; 5 6 int n,k,a[N]; 7 8 int main() 9 { 10 scanf("%d%d",&n,&k); 11 for (int i=1; i<=n; ++i) 12 scanf("%d",&a[i]); 13 for (int i=1; i<=n-1; ++i) 14 if (a[i]+a[i+1]>=k) 15 { 16 puts("Possible"); 17 for (int j=1; j<=i-1; ++j) printf("%d\n",j); 18 for (int j=n-1; j>=i+1; --j) printf("%d\n",j); 19 printf("%d\n",i); 20 return 0; 21 } 22 puts("Impossible"); 23 }
D - Stamp Rally
题意：
给出一个无向连通图。有q次询问，每次询问给出两个不同的点x和y，还有一个参数z。求一个最小的k使得只考虑编号不超过k的边时，x的连通块和y的连通块中的点的并集不小于z。
题解：
并查集+倍增好题……首先我们可以二分一个答案mid，那么问题就转化成了求一个点走的边的编号不超过mid的最大连通块。
加第i条边的时候，新创建一个权值为i的点new，然后把x,y两个点并查集的根连到new。可以发现这样做原图中的点都是叶子节点。那么加入第i条边的时候new点子树的叶子数量就是连通块大小。还发现点权值是随着点的深度减小而增大的，所以就可以用倍增了。
具体做法就是让x,y两个点跳到权值<=mid的最高的地方。如果两个点跳到一起就说明他们两个在可以用编号<=mid的边的时候处于同一个连通块，连通块大小就是这个点的size，否则就是两个点分别跳到的点的size之和。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N (200000+1000) 5 using namespace std; 6 7 int n,m,x,y,z,q,tot,Father[N],f[N][20],Size[N],Val[N]; 8 9 int Find(int x) 10 { 11 if (x==Father[x]) return x; 12 return Father[x]=Find(Father[x]); 13 } 14 15 int Get(int x,int v) 16 { 17 for (int i=16; i>=0; --i) 18 if (f[x][i] && Val[f[x][i]]<=v) x=f[x][i]; 19 return x; 20 } 21 22 int check(int x,int y,int v) 23 { 24 int xx=Get(x,v), yy=Get(y,v); 25 if (xx==yy) return Size[xx]; 26 return Size[xx]+Size[yy]; 27 } 28 29 int main() 30 { 31 scanf("%d%d",&n,&m); tot=n; 32 for (int i=1; i<=n*2; ++i) 33 Father[i]=i,Size[i]=(i<=n); 34 for (int i=1; i<=m; ++i) 35 { 36 scanf("%d%d",&x,&y); 37 if (Find(x)==Find(y)) continue; 38 int fx=Find(x), fy=Find(y); 39 tot++; Val[tot]=i; 40 Size[tot]=Size[fx]+Size[fy]; 41 Father[fx]=Father[fy]=f[fx][0]=f[fy][0]=tot; 42 } 43 for (int j=1; j<=16; ++j) 44 for (int i=1; i<=2*n; ++i) 45 f[i][j]=f[f[i][j-1]][j-1]; 46 scanf("%d",&q); 47 for (int i=1; i<=q; ++i) 48 { 49 scanf("%d%d%d",&x,&y,&z); 50 int l=1, r=m, ans=-1; 51 while (l<=r) 52 { 53 int mid=(l+r)>>1; 54 if (check(x,y,mid)>=z) ans=mid,r=mid-1; 55 else l=mid+1; 56 } 57 printf("%d\n",ans); 58 } 59 }
E - Candy Piles
是个博弈论好像……不会博弈论先坑着
F - Leftmost Ball
题意：
有n种颜色的球，每种m个。现在把这n*m个球排成一排，然后把每种颜色的最左边的球染成第n+1种颜色。求最终的颜色序列有多少种，对1000000007取模。
题解：
感觉是一道思维难度很高的题目……不过看了题解就发现这个思路并不难懂……
QAQ这个模型转化可以说是非常巧妙了 我怎么什么都不会啊……我怎么这么菜啊……
1 #include<iostream> 2 #include<cstdio> 3 #define N (5000000+1000) 4 #define MOD (1000000007) 5 using namespace std; 6 7 long long n,k,fac[N],inv[N],facinv[N],f[2010][2010]; 8 9 void Init() 10 { 11 fac[0]=inv[1]=facinv[0]=1; 12 for (int i=1; i<=5000000; ++i) 13 { 14 if (i!=1) inv[i]=(MOD-MOD/i)*inv[MOD%i]%MOD; 15 fac[i]=fac[i-1]*i%MOD; 16 facinv[i]=facinv[i-1]*inv[i]%MOD; 17 } 18 } 19 20 long long C(long long n,long long m) 21 { 22 if (m>n) return 0; 23 return fac[n]*facinv[m]%MOD*facinv[n-m]%MOD; 24 } 25 26 int main() 27 { 28 Init(); 29 scanf("%lld%lld",&n,&k); 30 if (k==1){printf("1"); return 0;} 31 f[0][1]=1; 32 for (int i=2; i<=n; ++i) 33 f[0][i]=f[0][i-1]*C(i*(k-1)-1,k-2)%MOD; 34 for (int i=1; i<=n; ++i) 35 for (int j=i; j<=n; ++j) 36 f[i][j]=(f[i-1][j]+f[i][j-1]*C(i+j*(k-1)-1,k-2)%MOD)%MOD; 37 printf("%lld",f[n][n]*fac[n]%MOD); 38 }
AGC 003
C - BBuBBBlesort!
题意：
有两种操作，第一种是将相邻的两个元素反转，第二种是将相邻的三个元素反转，尽可能多的使用第二种操作将序列排序，输出最小的第一种操作的使用次数。
题解：
可以发现二操作的本质就是随意交换下标奇偶性相同的位置。那么只需要统计一个数到自己应该到的位置的距离为奇数的个数最后除二就好了。
1 #include<iostream> 2 #include<cstdlib> 3 #include<cstdio> 4 #include<map> 5 #include<algorithm> 6 #define N (100000+1000) 7 using namespace std; 8 9 int n,a[N],b[N],ans; 10 map<int,int>pos; 11 12 int main() 13 { 14 scanf("%d",&n); 15 for (int i=1; i<=n; ++i) 16 scanf("%d",&a[i]),b[i]=a[i]; 17 sort(b+1,b+n+1); 18 for (int i=1; i<=n; ++i) pos[b[i]]=i; 19 for (int i=1; i<=n; ++i) 20 if (abs(i-pos[a[i]])%2) 21 ans++; 22 printf("%d",ans>>1); 23 }
D - Anticube
题意：
给你一个序列，让你尽可能选择多的数，让被选择的数两两间的乘积不为立方数。
题解：
感觉最近降智有点严重啊……
首先很容易想到的是，对于一个数，我们可以将其质因数分解然后将所有质因子的质数模3。
假设在指数模3下，一个数是$p_1^2*p_2^1*p_3^2$，设它的补数为$p_1^1*p_2^2*p_3^1$，显然一个数只有与其唯一的补数相乘才可以得到立方数。所以重点是如何将一个数质因子指数取模后的数及其补数求出来，最后答案只需要统计一个数和其补数的出现次数的max就好了。
显然若一个数的某个立方质因子的指数>=3，那么这个立方质因子的范围一定小于$10^{\frac{10}{3}}$，直接化简即可。至于一个数的补数，可以用$10^{\frac{10}{3}}$以内的素数去分解这个数，分解完后剩下的质因子个数肯定小于等于2，判断一下是否是完全平方数就好了。
1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<map> 5 #define N (100000+1000) 6 using namespace std; 7 8 long long vis[N],prime[N],cnt,n,x,a[N],b[N]; 9 map<long long,long long>Map; 10 11 void Init(int n) 12 { 13 for (int i=2; i<=n; ++i) 14 { 15 if (!vis[i]) prime[++cnt]=i; 16 for (int j=1; j<=cnt && prime[j]*i<=n; ++j) 17 { 18 vis[prime[j]*i]=true; 19 if (i%prime[j]==0) break; 20 } 21 } 22 } 23 24 long long sqr(long long x){return x*x;} 25 26 int main() 27 { 28 Init(3000); 29 scanf("%lld",&n); 30 for (int i=1; i<=n; ++i) 31 { 32 scanf("%lld",&x); 33 for (int j=1; j<=cnt; ++j) 34 while (x%(prime[j]*prime[j]*prime[j])==0) 35 x/=prime[j]*prime[j]*prime[j]; 36 Map[x]++; a[i]=x; 37 long long y=1; 38 for (int j=1; j<=cnt; ++j) 39 if (x%prime[j]==0) 40 { 41 y*=x%(prime[j]*prime[j])==0?prime[j]:prime[j]*prime[j]; 42 while (x%prime[j]==0) x/=prime[j]; 43 } 44 if (sqr((long long)sqrt(x))==x) y*=(long long)sqrt(x); 45 else y*=x*x; 46 b[i]=y; 47 } 48 long long ans=0; 49 if (Map[1]) ans++; Map[1]=0; 50 for (int i=1; i<=n; ++i) 51 { 52 ans+=max(Map[a[i]],Map[b[i]]); 53 Map[a[i]]=Map[b[i]]=0; 54 } 55 printf("%lld",ans); 56 }
E - Sequential operations on Sequence
题意：
有一个数字串S，初始长度为n，是1 2 3 4 …… n。
有m次操作，每次操作给你一个正整数a[i]，你先把S无穷重复，然后把前a[i]截取出来成为新的S。
求m次操作后，每个数字在S中出现的次数。
题解：
心态崩了
发现如果$a_i>=a_{i+1}$的话那么$a_i$显然是无效的。所以可以把m次操作处理一下成为一个递增序列。之前看的题解大小于号写反了让我懵逼了半天……
设m个操作处理后剩下了cnt个，设f[i]表示第i次操作后形成的序列在最终序列中出现了几次，显然f[cnt]=1
设初始长度为3，有两个操作7,15（如下图）。为了方便计算，将3也算入一个操作。
可以发现，长度变成15后的结果=长度为7时候的结果*2+余下的一点边角料。而长为x边角料显然与总序列前x是相等的。（如下图）
当边角料大于最小的操作的时候，显然是可以不停拆分的。就这样不停拆分直到边角料小于最小的操作，这样的话就可以直接将其加入对最终答案的贡献了。 实在不理解可以举几个例子对着代码模拟，然后感性理解一下反正我就是这么干的
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #define N (100000+1000) 5 using namespace std; 6 7 long long n,q,x,a[N],cnt,ans[N],f[N]; 8 9 int main() 10 { 11 scanf("%lld%lld",&n,&q); 12 a[++cnt]=n; 13 for (int i=1; i<=q; ++i) 14 { 15 scanf("%lld",&x); 16 while (cnt && x<=a[cnt]) --cnt; 17 a[++cnt]=x; 18 } 19 f[cnt]=1; 20 for (int i=cnt; i>=1; --i) 21 { 22 long long k=a[i]; 23 while (k>a[1]) 24 { 25 long long t=lower_bound(a+1,a+cnt+1,k)-a-1; 26 f[t]+=f[i]*(k/a[t]); 27 k%=a[t]; 28 } 29 ans[k]+=f[i]; 30 } 31 for (int i=n; i>=1; --i) ans[i]+=ans[i+1]; 32 for (int i=1; i<=n; ++i) printf("%lld\n",ans[i]); 33 }
F - Fraction of Fractal
题意：
题解：
我 不会 抄别人的题解
不过思路还是很妙的
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N (1000+100) 5 #define MOD (1000000007) 6 using namespace std; 7 8 long long n=2,H,W,K,a[N][N],S,A,B; 9 char st[N]; 10 11 struct Matrix 12 { 13 long long m[11][11]; 14 void clear(){memset(m,0,sizeof(m));} 15 }V,G; 16 Matrix operator * (Matrix a,Matrix b) 17 { 18 Matrix ans; ans.clear(); 19 for (int i=1; i<=n; ++i) 20 for (int j=1; j<=n; ++j) 21 for (int k=1; k<=n; ++k) 22 (ans.m[i][j]+=a.m[i][k]*b.m[k][j])%=MOD; 23 return ans; 24 } 25 26 Matrix Mat_Qpow(Matrix a,long long p) 27 { 28 Matrix ans; ans.clear(); 29 for (int i=1; i<=n; ++i) ans.m[i][i]=1; 30 while (p) 31 { 32 if (p&1) ans=ans*a; 33 a=a*a; p>>=1; 34 } 35 return ans; 36 } 37 38 long long Qpow(long long a,long long p) 39 { 40 long long ans=1,base=a; 41 while (p) 42 { 43 if (p&1) ans=ans*base%MOD; 44 base=base*base%MOD; p>>=1; 45 } 46 return ans; 47 } 48 49 int Get() 50 { 51 int ansLR=0,ansUD=0;//左右,上下接口的个数 52 int anslr=0,ansud=0;//左右，上下相邻为黑色的位置个数 53 for (int i=1; i<=H-1; ++i) 54 for (int j=1; j<=W; ++j) 55 ansud+=(a[i][j]==a[i+1][j] && a[i][j]==1); 56 for (int i=1; i<=H; ++i) 57 for (int j=1; j<=W-1; ++j) 58 anslr+=(a[i][j]==a[i][j+1] && a[i][j]==1); 59 60 for (int i=1; i<=H; ++i) 61 ansLR+=(a[i][1]==a[i][W] && a[i][1]==1); 62 for (int i=1; i<=W; ++i) 63 ansUD+=(a[1][i]==a[H][i] && a[1][i]==1); 64 65 if (ansLR && ansUD) return 2; 66 if (!ansLR && !ansUD) return 0; 67 if (ansLR) A=anslr,B=ansLR; 68 else A=ansud,B=ansUD; 69 return 1; 70 } 71 72 void Solve() 73 { 74 G.m[1][1]=A; G.m[1][2]=B; 75 V.m[1][1]=S; V.m[1][2]=0; V.m[2][1]=A; V.m[2][2]=B; 76 V=Mat_Qpow(V,K-2); G=G*V; 77 long long ans=Qpow(S,K-1)-G.m[1][1]; 78 printf("%lld",ans<=0?ans+MOD:ans); 79 } 80 81 int main() 82 { 83 scanf("%lld%lld%lld",&H,&W,&K); 84 for (int i=1; i<=H; ++i) 85 { 86 scanf("%s",st); 87 for (int j=1; j<=W; ++j) 88 a[i][j]=st[j-1]=='#'?1:0,S+=a[i][j]; 89 } 90 if (K<=1) {printf("1"); return 0;} 91 92 int g=Get(); 93 if (g==0) {printf("%lld",Qpow(S,K-1)); return 0;} 94 if (g==2) {printf("1"); return 0;} 95 96 Solve(); 97 }
AGC 004
C - AND Grid
题意：
给你一个网格图，有一些地方涂了色。构造两个涂色后的四连通的网格图，让他们涂色的方格的并集等于给定网格图。给定的网格图第一行，第一列，第n行，第m列一定没有涂色。
题解：
设给定网格图为s，构造的为a和b。初始值要设成a=s且b=s。然后a涂第一行，b涂最后一行。最后a涂奇数列，b涂偶数列。正确性显然。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N (601) 5 using namespace std; 6 int n,m; 7 char s[N][N],a[N][N],b[N][N]; 8 int main() 9 { 10 scanf("%d%d",&n,&m); 11 for (int i=1; i<=n; ++i) 12 scanf("%s",s[i]+1); 13 memcpy(a,s,sizeof(s)); memcpy(b,s,sizeof(s)); 14 for (int i=1; i<=m; ++i) 15 a[1][i]=b[n][i]='#'; 16 for (int j=2; j<m; ++j) 17 for (int i=2; i<n; ++i) 18 if (j%2) a[i][j]='#'; 19 else b[i][j]='#'; 20 for (int i=1; i<=n; ++i) 21 { 22 for (int j=1; j<=m; ++j) 23 printf("%c",a[i][j]); 24 printf("\n"); 25 } 26 printf("\n"); 27 for (int i=1; i<=n; ++i) 28 { 29 for (int j=1; j<=m; ++j) 30 printf("%c",b[i][j]); 31 printf("\n"); 32 } 33 }
D - Teleporter
题意：
一个n个点的基环树（环套树），更改尽量少的出边，使得任意点在走k步后都正好到达1节点。
题解：
一开始sb了从上往下贪结果只能过三分之二……首先很容易发现1的出边必须连向自己，否则一定无法满足条件，这样就形成了一棵树。接下来只需要从下往上贪心，若一个点最长能往下延伸k-1条边，就将这个点断掉连向1节点。
1 #include<iostream> 2 #include<cstdio> 3 #define N (100000+1000) 4 using namespace std; 5 6 struct Edge{int to,next;}edge[N<<1]; 7 int n,k,a[N],head[N],num_edge,ans; 8 9 void add(int u,int v) 10 { 11 edge[++num_edge].to=v; 12 edge[num_edge].next=head[u]; 13 head[u]=num_edge; 14 } 15 16 int Dfs(int x) 17 { 18 int maxn=0; 19 for (int i=head[x]; i; i=edge[i].next) 20 maxn=max(maxn,Dfs(edge[i].to)+1); 21 if (maxn==k-1 && x!=1 && a[x]!=1){maxn=-1; ans++;} 22 return maxn; 23 } 24 25 int main() 26 { 27 scanf("%d%d%d",&n,&k,&a[1]); 28 for (int i=2; i<=n; ++i) 29 scanf("%d",&a[i]),add(a[i],i); 30 Dfs(1); 31 printf("%d",ans+(a[1]!=1)); 32 }
E - Salvage Robots
题意：
一个n*m的网格，上面有一些机器人，还有一个出口。每次可以将所有机器人向一个方向一起移动一格，机器人出界就挂了。问最多多少个机器人能逃离。
题解：
我再也不做边界处理这么恶心的题了……
可以换一种想法，把移动机器人转换成移动出口。设f[i][j][k][l]表示取完了(i,j)(k,l)这个矩形内的机器人，然后分别向四周扩展即可。开short不会MLE。
1 #include<iostream> 2 #include<cstdio> 3 #define N (100+5) 4 using namespace std; 5 6 short n,m,line[N][N],list[N][N],f[N][N][N][N],sx,sy,ans; 7 char a[N][N]; 8 9 int main() 10 { 11 cin>>n>>m; 12 for (short i=1; i<=n; ++i) 13 scanf("%s",a[i]+1); 14 for (short i=1; i<=n; ++i) 15 for (short j=1; j<=m; ++j) 16 { 17 line[i][j]=line[i][j-1]+(a[i][j]=='o'); 18 list[i][j]=list[i-1][j]+(a[i][j]=='o'); 19 if (a[i][j]=='E') sx=i,sy=j; 20 } 21 for (short i=sx; i>=1; --i) 22 for (short j=sy; j>=1; --j) 23 for (short k=sx; k<=n; ++k) 24 for (short l=sy; l<=m; ++l) 25 { 26 if (i>1 && k-sx<i-1) 27 f[i-1][j][k][l]=max((int)f[i-1][j][k][l],f[i][j][k][l]+line[i-1][min((int)l,m-sy+j)]-line[i-1][max(j-1,l-sy)]); 28 if (k<n && sx+k<n+i) 29 f[i][j][k+1][l]=max((int)f[i][j][k+1][l],f[i][j][k][l]+line[k+1][min((int)l,m-sy+j)]-line[k+1][max(j-1,l-sy)]); 30 if (j>1 && l-sy<j-1) 31 f[i][j-1][k][l]=max((int)f[i][j-1][k][l],f[i][j][k][l]+list[min((int)k,n-sx+i)][j-1]-list[max(i-1,k-sx)][j-1]); 32 if (l<m && sy+l<m+j) 33 f[i][j][k][l+1]=max((int)f[i][j][k][l+1],f[i][j][k][l]+list[min((int)k,n-sx+i)][l+1]-list[max(i-1,k-sx)][l+1]); 34 short now=max(max(f[i-1][j][k][l],f[i][j][k+1][l]),max(f[i][j-1][k][l],f[i][j][k][l+1])); 35 ans=max(ans,now); 36 } 37 cout<<ans; 38 }
F - Namori
题意：
给你一颗全白的树或环套树，每次可以选择一条连接两个同色点的边，将两个端点反色。 问变成全黑的最小步数，可能无解。
题解：
这个题好神啊……说白了就是我不会
树的做法
首先要知道，树是一个二分图。
不妨设深度为奇数（根的深度为1）的点放着一枚硬币，深度为偶数的点是空位。这样就将树染成了一个二分图
可以发现，一次操作相当于将一个硬币运到相邻的空位。
当初始硬币处都是空位，初始空位处都是硬币的时候，对应着原图全被染成黑色。
显然只有硬币数=空位数才有解。
设硬币权值为1，空位权值为-1,1为树的根。
size对应子树和，当size[1]=0的时候就对应着有解。
因为size[i]是一定要从i这个子树运出或运入的硬币数，这个值代表边(i,fa[i])被操作的次数。
$\sum_{i=1}^nabs(size_i)$显然是答案下界。我们尽可能让子树内配对就可以达到这个答案下界了。
奇环的做法
任意断开环上一条边(x,y)，让基环树变成一棵树。
x点和y点一定属于二分图同一侧（要么都是空位，要么都是硬币）。
操作一次边(x,y)，就是在x点和y点各增加或减少一枚硬币。
计算出变成树后多余或者缺少的硬币数量，是奇数就无解，否则就用操作(x,y)边让硬币数=空缺数。
同时修改x和y的祖先的size值。因为操作边(x,y)平衡硬币数和空位数的时候，他们的祖先对应的操作次数也会改变。
最后按照树的做法做一遍就可以了。
偶环的做法
任意断开一边(x,y)，设(x,y)这条边被操作了k次(可以为负数)，那么x上增加了k个硬币，y上减少了k个。
这k次操作显然只影响路径x-lca(x,y)-y的操作数(也就是size)。
为了最小化答案，显然是对这个路径上每一个位置的size加上或减去一个数，最小化这条路径上点的size和。
这个加上/减去的数显然就是路径上所有size的中位数。
1 #include<iostream> 2 #include<cstdlib> 3 #include<cstdio> 4 #include<algorithm> 5 #define N (200000+1000) 6 using namespace std; 7 8 struct Edge{int to,next;}edge[N<<1]; 9 int n,m,u,v,K[N],size[N],sum,cir_num; 10 int Depth[N],Father[N],stx,sty; 11 int head[N],num_edge,Q[N],tot; 12 13 void add(int u,int v) 14 { 15 edge[++num_edge].to=v; 16 edge[num_edge].next=head[u]; 17 head[u]=num_edge; 18 } 19 20 void Dfs(int x,int fa) 21 { 22 Depth[x]=Depth[fa]+1; 23 Father[x]=fa; 24 size[x]=(Depth[x]%2==1)?1:-1; 25 for (int i=head[x]; i; i=edge[i].next) 26 if (edge[i].to!=fa) 27 { 28 if (Depth[edge[i].to]) 29 { 30 if (cir_num) continue; 31 stx=x; sty=edge[i].to; 32 cir_num=Depth[x]-Depth[edge[i].to]+1; 33 continue; 34 } 35 Dfs(edge[i].to,x); 36 size[x]+=size[edge[i].to]; 37 } 38 } 39 40 int LCA(int x,int y) 41 { 42 while (Depth[x]>Depth[y]) x=Father[x]; 43 while (Depth[y]>Depth[x]) y=Father[y]; 44 while (x!=y) x=Father[x], y=Father[y]; 45 return x; 46 } 47 48 int main() 49 { 50 scanf("%d%d",&n,&m); 51 for (int i=1; i<=m; ++i) 52 { 53 scanf("%d%d",&u,&v); 54 add(u,v); add(v,u); 55 } 56 Dfs(1,0); 57 58 if (m==n-1)//tree 59 { 60 if (size[1]!=0) printf("-1"); 61 else 62 { 63 long long ans=0; 64 for (int i=1; i<=n; ++i) ans+=abs(size[i]); 65 printf("%lld",ans); 66 } 67 } 68 else if (cir_num%2==1)//odd 69 { 70 if (abs(size[1])%2==1) printf("-1"); 71 else 72 { 73 int sz=size[1]; 74 long long ans=abs(sz)/2; 75 for (int u=stx; u; u=Father[u]) size[u]-=sz/2; 76 for (int u=sty; u; u=Father[u]) size[u]-=sz/2; 77 for (int i=1; i<=n; ++i) ans+=abs(size[i]); 78 printf("%lld",ans); 79 } 80 } 81 else//even 82 { 83 if (size[1]!=0) printf("-1"); 84 else 85 { 86 int lca=LCA(stx,sty); 87 for (int u=stx; u!=lca; u=Father[u]) 88 Q[++tot]=size[u],Depth[u]=0; 89 for (int u=sty; u!=lca; u=Father[u]) 90 Q[++tot]=-size[u],Depth[u]=0; 91 Q[++tot]=0; 92 sort(Q+1,Q+tot+1); 93 int sz=Q[tot/2]; 94 95 long long ans=0; 96 for (int i=1; i<=n; ++i) if (Depth[i]) ans+=abs(size[i]); 97 for (int i=1; i<=tot; ++i) ans+=abs(Q[i]-sz); 98 printf("%lld",ans); 99 } 100 } 101 }
AGC 005
C - Tree Restoring
题意：
给你每个点到离它最远的点的距离，问是否可以构造出这样的一颗树。
题解：
首先先把直径搞出来，看看能否满足直径，不能直接impossible。如果能满足直径的话，就判断一下剩下的点往直径上连边是否可行。（代码应该比较好懂）
1 #include<algorithm> 2 #include<cstdio> 3 #define N (1001) 4 using namespace std; 5 6 int n,a[N],Keg[N],maxn; 7 8 int main() 9 { 10 scanf("%d",&n); 11 for (int i=1; i<=n; ++i) 12 { 13 scanf("%d",&a[i]); 14 Keg[a[i]]++; maxn=max(maxn,a[i]); 15 } 16 for (int i=(maxn+1)/2; i<=maxn; ++i) 17 { 18 if (maxn%2==0 && i==(maxn+1)/2) Keg[i]-=1; 19 else Keg[i]-=2;//从桶里取点往直径上放。 20 if (Keg[i]<0){puts("Impossible"); return 0;}//桶里没点可放就无解 21 } 22 for (int i=1; i<=(maxn+1)/2; ++i)//显然若还存在最远距离在这个范围内的点，这个点是无法连到直径上的 23 if (Keg[i]) {puts("Impossible"); return 0;} 24 puts("Possible"); 25 }
D - ~K Perm Counting
题意：
给定n和k，a是n的一个排列。求满足$abs(a_i-i)!=k$的n的排列的个数。
题解：
这个题用到的算法都是基础难度但的确很妙啊_(:з」∠)_
首先令$g_i$表示有i个位置不合法的排列总数，简单容斥一下可以发现
$ans=n!+\sum_{i=0}^ng_i*(n-i)!*(-1)^i$。那么现在的问题在于怎么求$g_i$。
考虑一个n=6, k=2的情况，可以建立一个二分图，若i位置不能放数字x，就左边的i连右边的x。(如下)
发现这个二分图是由若干条链构成的？！
我们可以把这些链都提取出来，任意顺序摆成一条链的形状。
很显然我们选择一条红边就相当于选了一个不合法的情况，且不能连续选择两个连续的红边。
所以可以设f[i][j][0/1]表示前i个点选了j条红边，且第i个点不选/选红边。
DP的时候特判一下没有红边的地方即可。最后g[i]=f[2*n][i][0]+f[2*n][i][1]
1 #include<iostream> 2 #include<cstdio> 3 #define N (4009) 4 #define MOD (924844033) 5 using namespace std; 6 7 long long n,k,vis[N],tot,g[N],f[N][N][2],fac[N]; 8 9 int main() 10 { 11 scanf("%lld%lld",&n,&k); 12 for (int t=1; t<=2; ++t) 13 for (int i=1; i<=k; ++i) 14 for (int j=i; j<=n; j+=k) 15 { 16 tot++; 17 if (j!=i) vis[tot]=true; 18 } 19 f[0][0][0]=1; 20 for (int i=1; i<=2*n; ++i) 21 for (int j=0; j<=n; ++j) 22 { 23 (f[i][j][0]+=f[i-1][j][0]+f[i-1][j][1])%=MOD; 24 if (vis[i]) (f[i][j][1]+=f[i-1][j-1][0])%=MOD; 25 } 26 for (int i=0; i<=n; ++i) 27 g[i]=(f[2*n][i][0]+f[2*n][i][1])%MOD; 28 29 fac[0]=1; 30 for (int i=1; i<=n; ++i) fac[i]=fac[i-1]*i%MOD; 31 long long ans=fac[N]; 32 for (int i=0,j=1; i<=n; ++i,j=-j) 33 (ans+=j*g[i]*fac[n-i]%MOD)%=MOD; 34 printf("%lld\n",(ans+MOD)%MOD); 35 }
E - Sugigma: The Showdown
题意：
一个图有红边和蓝边，红边和蓝边各自构成一颗树。A走红边，B走蓝边。
两人相遇游戏结束，A想最大化时间，B想最小化时间。若B永远抓不到A输出-1
题解：
先将蓝树以B的起点为根建立出来，可以发现A在蓝树上的运动轨迹：要么往子树走，要么往兄弟走，要么往父亲走。那么B只需要一步一步往下沉底去抓A就可以了。
什么时候无解呢？当一条红边连接的两个点，在蓝树上距离超过2，且A能到达两点之一，那么A就可以通过反复横跳来让B永远抓不到。（这个地方的check函数写法我觉得挺巧妙的……）
否则就让A尽可能往蓝树的更深的点走，到一个最深的地方弃疗等B来抓它就好了。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cstdio> 5 #define N (400000+1000) 6 using namespace std; 7 8 struct edge{int to,next;}edge[N<<1]; 9 int Father[N],Depth[N],Dfn[N],Low[N]; 10 int n,x,y,S1,S2,cnt,ans,u[N],v[N]; 11 int head[N],num_edge; 12 13 void add(int u,int v) 14 { 15 edge[++num_edge].to=v; 16 edge[num_edge].next=head[u]; 17 head[u]=num_edge; 18 } 19 20 bool check(int x,int y) 21 { 22 if (Dfn[x]>Dfn[y]) swap(x,y); 23 if (Low[x]>=Low[y]) 24 return (Depth[y]-Depth[x])>2; 25 return Father[x]!=Father[y]; 26 } 27 28 void Dfs1(int x) 29 { 30 Dfn[x]=++cnt; Depth[x]=Depth[Father[x]]+1; 31 for (int i=head[x]; i; i=edge[i].next) 32 if (edge[i].to!=Father[x]) 33 { 34 Father[edge[i].to]=x; 35 Dfs1(edge[i].to); 36 } 37 Low[x]=++cnt; 38 } 39 40 void Dfs2(int x,int fa,int k) 41 { 42 ans=max(ans,Depth[x]-1); 43 if (Depth[x]<=k) return; 44 for (int i=head[x]; i; i=edge[i].next) 45 if (edge[i].to!=fa) 46 { 47 if (check(x,edge[i].to)) puts("-1"),exit(0); 48 Dfs2(edge[i].to,x,k+1); 49 } 50 } 51 52 int main() 53 { 54 scanf("%d%d%d",&n,&S1,&S2); 55 for (int i=1; i<=n-1; ++i) 56 scanf("%d%d",&u[i],&v[i]); 57 for (int i=1; i<=n-1; ++i) 58 scanf("%d%d",&x,&y),add(x,y),add(y,x); 59 Dfs1(S2); 60 memset(head,0,sizeof(head)); num_edge=0; 61 for (int i=1; i<=n-1; ++i) add(u[i],v[i]),add(v[i],u[i]); 62 Dfs2(S1,-1,1); 63 printf("%d\n",ans<<1); 64 }
F - Many Easy Problems
题意：
给你一颗n个点的树和一个整数k。设S为树上某些点的集合，定义f(S)为包含S的最小的联通子图大小。
n个点选k个点有$C_n^k$种方案，求所有方案的f(S)的和。
题解：
没仔细读题一开始质数设的998244353……然而应该是924844033。第一个原根是3，第二个原根是5……
如果直接对于每个k直接求解答案好像不可做……就分开考虑每个点对答案的贡献。
容斥一下可以得到i点对于一个k的答案的贡献为$C_n^k-\sum_{fa[j]=i}C_{size[j]}^k-C_{n-size[i]}^k$
显然每个子树只会对答案贡献一次，所以最终答案为
$ans[k]=n*C_{n}^{k}-\sum_{i=0}^{n}cnt[i]*C_{i}^{k}$，其中cnt[i]表示大小为i的子树个数。
将组合数拆开后发现是$\frac{i!}{(i-k)!k!}$，将其中的$i!$和cnt[i]放到一起，$k!$提到最前面，就成了
$n*C_{n}^{k}-k!*\sum_{i=0}^{n}cnt[i]*i!*\frac{1}{(i-k)!}$，发现那个sigma可以通过翻转数组的套路卷积求得。
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #define N (800000+1000) 5 #define MOD (924844033) 6 #define LL long long 7 using namespace std; 8 9 struct Edge{int to,next;}edge[N]; 10 LL n,u,v,fn=1,l,r[N],a[N],b[N]; 11 LL fac[N],inv[N],facinv[N],ans[N]; 12 LL head[N],num_edge,size[N],cnt[N]; 13 14 void add(int u,int v) 15 { 16 edge[++num_edge].to=v; 17 edge[num_edge].next=head[u]; 18 head[u]=num_edge; 19 } 20 21 LL Qpow(LL a,LL b) 22 { 23 LL ans=1,base=a; 24 while (b) 25 { 26 if (b&1) ans=ans*base%MOD; 27 base=base*base%MOD; b>>=1; 28 } 29 return ans; 30 } 31 32 void NTT(LL n,LL *a,LL opt) 33 { 34 for (int i=0; i<n; ++i) 35 if (i<r[i]) swap(a[i],a[r[i]]); 36 for (int k=1; k<n; k<<=1) 37 { 38 LL wn=Qpow(5,(MOD-1)/(k<<1)); 39 for (int i=0; i<n; i+=(k<<1)) 40 { 41 LL w=1; 42 for (int j=0; j<k; ++j, w=w*wn%MOD) 43 { 44 LL x=a[i+j], y=w*a[i+j+k]%MOD; 45 a[i+j]=(x+y)%MOD; a[i+j+k]=(x-y+MOD)%MOD; 46 } 47 } 48 } 49 if (opt==-1) 50 { 51 reverse(a+1,a+n); 52 LL g=Qpow(n,MOD-2); 53 for (int i=0; i<n; ++i) a[i]=a[i]*g%MOD; 54 } 55 } 56 57 void Init() 58 { 59 fac[0]=1; inv[1]=1; facinv[0]=1; 60 for (int i=1; i<=n; ++i) 61 { 62 if (i!=1) inv[i]=(MOD-MOD/i)*inv[MOD%i]%MOD; 63 fac[i]=fac[i-1]*i%MOD; facinv[i]=facinv[i-1]*inv[i]%MOD; 64 } 65 } 66 67 LL C(LL n,LL m) 68 { 69 if (n<m) return 0; 70 return fac[n]*facinv[m]%MOD*facinv[n-m]%MOD; 71 } 72 73 void Dfs(int x,int fa) 74 { 75 size[x]=1; 76 for (int i=head[x]; i; i=edge[i].next) 77 if (edge[i].to!=fa) 78 { 79 Dfs(edge[i].to,x); 80 cnt[size[edge[i].to]]++; 81 size[x]+=size[edge[i].to]; 82 } 83 cnt[n-size[x]]++; 84 } 85 86 int main() 87 { 88 scanf("%lld",&n); 89 for (int i=1; i<=n-1; ++i) 90 { 91 scanf("%lld%lld",&u,&v); 92 add(u,v); add(v,u); 93 } 94 Init(); Dfs(1,-1); 95 while (fn<=n+n) fn<<=1, l++; 96 for (int i=0; i<=fn; ++i) 97 r[i]=(r[i>>1]>>1)|((i&1)<<(l-1)); 98 99 for (int i=0; i<=n; ++i) a[i]=cnt[i]*fac[i]%MOD; 100 for (int i=0; i<=n; ++i) b[i]=facinv[n-i]; 101 102 NTT(fn,a,1); NTT(fn,b,1); 103 for (int i=0; i<=fn; ++i) a[i]=a[i]*b[i]%MOD; 104 NTT(fn,a,-1); 105 106 for (int i=1; i<=n; ++i) ans[i]=-facinv[i]*a[n+i]%MOD+MOD; 107 for (int i=1; i<=n; ++i) ans[i]=(ans[i]+n*C(n,i))%MOD; 108 for (int i=1; i<=n; ++i) printf("%lld\n",ans[i]); 109 }
AGC 006
C - Rabbit Exercise
题意：
数轴上有n只兔子，第i只兔子的坐标为xi。 有一组操作，这组操作的第i个操作是要让第ai只兔子等概率的跳到自己关于第ai+1或第ai-1只兔子的对称点。 进行K组操作，求每只兔子最后坐标的期望值。
题解：
操作编号为x的时候期望结果为$a_x=\frac{1}{2}(2a_{x-1}-a_x)+\frac{1}{2}(2a_{x+1}-a_x)=a_{x-1}+a_{x+1}-a_x$。
感性理解一下又由于期望的线性发现可以直接用期望到达的点位置替换当前点位置。
对原数组$a_1,a_2,a_3$做一个差分，就成了$a_1-a_0,a_2-a_1,a_3-a_2$。
对原数组中间的求期望，结果是$a_1,a_1+a_3-a_2$，对应到差分数组上就是$a_1-a_0,a_3-a_2,a_2-a_1$。
可以发现操作i位置的时候就相当于swap差分数组的i和i+1位置。
然后这里有一个非常巧妙的方法，就是用类似快速幂的倍增思想。
可以用一个点跳一组后的位置求跳两组后的位置……然后将k二进制分解就好了。感觉看代码可能好懂一点
1 #include<iostream> 2 #include<cstdio> 3 #define N (100000+1000) 4 using namespace std; 5 6 long long n,m,k,x,a[N],d[N],p[N],ans[N],tmp[N],sum; 7 8 void Qpow() 9 { 10 for (int i=1; i<=n; ++i) ans[i]=i; 11 while (k) 12 { 13 if (k&1) 14 { 15 for (int i=1; i<=n; ++i) tmp[i]=ans[p[i]]; 16 for (int i=1; i<=n; ++i) ans[i]=tmp[i]; 17 } 18 for (int i=1; i<=n; ++i) tmp[i]=p[p[i]]; 19 for (int i=1; i<=n; ++i) p[i]=tmp[i]; 20 k>>=1; 21 } 22 for (int i=1; i<=n; ++i) tmp[i]=d[ans[i]]; 23 for (int i=1; i<=n; ++i) d[i]=tmp[i]; 24 } 25 26 int main() 27 { 28 scanf("%lld",&n); 29 for (int i=1; i<=n; ++i) 30 scanf("%lld",&a[i]),d[i]=a[i]-a[i-1],p[i]=i; 31 scanf("%lld%lld",&m,&k); 32 for (int i=1; i<=m; ++i) 33 scanf("%lld",&x),swap(p[x],p[x+1]); 34 Qpow(); 35 for (int i=1; i<=n; ++i) 36 sum+=d[i],printf("%lld\n",sum); 37 }
D - Median Pyramid Hard
题意：
初始给你一个左边这样的塔，从下往上填。一个空位置的取值等于他下面三个格子的中位数，问塔顶的数是多少。
题解：
首先二分塔顶的数，把最底层大于的转化为1，小于等于的转化为0。
画个图可以发现底层连续的0或者1一定是能延续到达塔顶的。如果两者只存在一个那好说，存在啥塔顶就是啥。
如果连续的0和连续的1都存在，就看谁离中间更近，近的那个能够占领塔顶。
注意判断没有连续0和连续的1的情况。
1 #include<iostream> 2 #include<cstdio> 3 #define N (200000+1000) 4 using namespace std; 5 6 int n,a[N],l,r,ans,mid; 7 int g(int x){return x>mid;} 8 9 int check(int x) 10 { 11 for (int i=1; i<=n-1; ++i) 12 { 13 if (g(a[n-i+1])==g(a[n-i])) return g(a[n-i]); 14 if (g(a[n+i-1])==g(a[n+i])) return g(a[n+i]); 15 } 16 return g(a[1]); 17 } 18 19 int main() 20 { 21 scanf("%d",&n); 22 for (int i=1; i<=2*n-1; ++i) 23 scanf("%d",&a[i]); 24 l=1; r=2*n-1; ans=-1; 25 while (l<=r) 26 { 27 mid=(l+r)>>1; 28 if (!check(mid)) ans=mid,r=mid-1; 29 else l=mid+1; 30 } 31 printf("%d\n",ans); 32 }
E - Rotate 3x3
题意：
给你一个3*n的网格图，网格内乱序有一些数。有一个操作是选择一个3*3的矩形将其旋转180°。问是否能将其操作成为第一列为1,2,3，第二列为4,5,6……的网格。
题解：
对这种奇怪的结论题不是很感冒……感觉这篇题解写的挺清楚的……不过感觉复杂度有点假……也不知道是我感性理解错了还是难造数据……
1 #include<iostream> 2 #include<cstdlib> 3 #include<cstdio> 4 #define N (100000+1000) 5 using namespace std; 6 7 int n,a[4][N],tar[N],R[2]; 8 9 int main() 10 { 11 scanf("%d",&n); 12 for (int i=1; i<=3; ++i) 13 for (int j=1; j<=n; ++j) 14 scanf("%d",&a[i][j]); 15 for (int i=1; i<=n; ++i) 16 { 17 tar[i]=(a[2][i]+1)/3; 18 if (abs(a[2][i]-a[1][i])!=1 || abs(a[2][i]-a[3][i])!=1 || abs(i-tar[i])%2) 19 { 20 puts("No"); 21 return 0; 22 } 23 R[i&1]^=(a[1][i]>a[3][i]); 24 } 25 for (int i=1; i<=n; ++i) 26 while (tar[i]!=i) 27 R[i&1^1]^=1,swap(tar[i],tar[tar[i]]); 28 if (R[0] || R[1]) puts("No"); 29 else puts("Yes"); 30 }
F - Blackout
题意：
给你一个n*n的网格，初始给你m个点让你把他们染黑，同时若点(x,y)和(y,z)为黑色，则(z,x)也被染成黑色。问最后有多少个黑色的格子。
题解：
将模型转化到图上，题目就成了若存在边x->y, y->z，则会产生一条z->x的边。这显然是一个三元环。
我们可以发现两个弱连通分量之间肯定不会产生边，所以可以对于每个弱连通分量单独考虑。
对于当前弱连通分量，我们可以进行三染色染成0,1,2。若染色过程中：
1、颜色没有全用到。那么一个点最多只能往后走一步，肯定走不出来一个三元环，所以这个弱连通分量对答案的贡献是这个分量的边数。
2、成功染色。任意0点可以连任意1点，任意1点可以连任意2点，任意2点可以连任意0点。
3、染色失败，即一个点可能存在多种颜色。这样的话图中肯定是存在环的（自环也考虑）。画画图可以发现这个图可以被连成一个完全图，两两之间任意连边即可。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N (100000+1000) 5 using namespace std; 6 7 struct Edge{int to,next,len;}edge[N<<1]; 8 int head[N],num_edge; 9 long long n,m,cnt[4],flag,col[N]; 10 long long ans,sum,u,v; 11 12 void add(int u,int v,int l) 13 { 14 edge[++num_edge].to=v; 15 edge[num_edge].next=head[u]; 16 edge[num_edge].len=l; 17 head[u]=num_edge; 18 } 19 20 void Dfs(int x) 21 { 22 cnt[col[x]]++; 23 for (int i=head[x]; i; i=edge[i].next) 24 { 25 if (edge[i].len==1) sum++; 26 if (col[edge[i].to]==-1) 27 { 28 col[edge[i].to]=(col[x]+edge[i].len)%3; 29 Dfs(edge[i].to); 30 } 31 else if (col[edge[i].to]!=(col[x]+edge[i].len)%3) 32 flag=1; 33 } 34 } 35 36 int main() 37 { 38 scanf("%lld%lld",&n,&m); 39 for (int i=1; i<=m; ++i) 40 { 41 scanf("%lld%lld",&u,&v); 42 add(u,v,1); add(v,u,2); 43 } 44 memset(col,-1,sizeof(col)); 45 for (int i=1; i<=n; ++i) 46 if (col[i]==-1) 47 { 48 sum=cnt[0]=cnt[1]=cnt[2]=flag=0; 49 col[i]=0; Dfs(i); 50 if (flag) 51 { 52 long long t=cnt[0]+cnt[1]+cnt[2]; 53 ans+=t*t; continue; 54 } 55 if ((!cnt[0]) || (!cnt[1]) || (!cnt[2])) 56 { 57 ans+=sum; 58 continue; 59 } 60 ans+=cnt[0]*cnt[1]+cnt[1]*cnt[2]+cnt[2]*cnt[0]; 61 } 62 printf("%lld\n",ans); 63 }
AGC 007
C - Pushing Balls
题意：
平面上有n+1个洞，洞两两之间有一个球，球洞之间的距离是一个等差数列。每次随机选择一个不在洞中的球，再随机选择一个方向推动，球将会进入这个方向上距离它最近的没有球的洞。滚动过程中若洞内有球则会滚过这个洞，求球都进洞后球滚动的期望总距离。
题解：
期望题 不会 抄题解
1 #include<cstdio> 2 double n,d,x,a; 3 int main() 4 { 5 scanf("%lf%lf%lf",&n,&d,&x); 6 for (d+=(n-0.5)*x; n; --n)a+=d,d+=d/n; 7 printf("%.10f",a); 8 }
D - Shik and Game
题意：
有一个在数轴上进行的游戏。最初，玩家位于位置0，拥有N颗糖果，出口在位置E。游戏中还有N只熊。第i熊位于$x_i$。玩家的移动速度是1，而熊不移动。
当玩家给熊一个糖果时，熊会在单位时间T后提供一个硬币。更具体地说，如果第i只熊在时刻t得到一个糖果，它将在时刻t+T在它的位置放一个硬币。这个游戏的目的是给所有的熊糖果，拿起所有的硬币，然后去出口。注意，玩家只能给每只熊一个糖果。此外，每个熊只会产生一次硬币。如果玩家访问一个有硬币的位置，玩家可以捡起硬币。硬币不会消失，直到被玩家收集。
计算收集所有硬币到出口的最短时间。
题解：
挺有意思的一个题……比C简单多了
首先观察下数据范围和部分分，基本DP没跑了。显然最优情况一定是每次放置一段连续的糖果然后取走这一段产生的硬币。
那么$n^2$的DP就比较容易了，设f[i]表示取完前i个位置的硬币的最小花费时间，转移也显然。
$f[i]=min(f[i],f[j-1]+(a[i]-a[j])*3+(a[j]-a[j-1])+max(0,t-2*(a[i]-a[j])))$
这个转移表示连续在熊$[j,i]$放了糖果并取走了他们的硬币。行走路径为$j->i->j->i$。后面的max项是为了防止出现从$i->j$后j位置的熊硬币还没有产出的情况。
考虑如何优化呢？发现式子里面有和i相关的项有和j相关的项，$i$相关的项可以在计算$f[i]$的时候直接加进去，而$j$相关的项我们可以插入到线段树中，更新的时候直接区间查询最小值。
还有一个问题，因为DP式子里有一个max，所以没法直接插入到线段树中。但是可以发现，当计算$f[i]$的时候，前面的答案一定能分成两部分，前半部分1~？是离$a[i]$足够远所以计算的时候没有max的，后半部分?+1~i-1是离i太近而计算的时候考虑max的。而这个分界线又是一直向右挪的，所以用个指针记录一下，当一个位置离当前转移位置足够远的时候就把max的影响消除掉。（注释掉的是部分分代码）
1 #include<iostream> 2 #include<cstdio> 3 #define N (100000+1000) 4 #define LL long long 5 using namespace std; 6 7 LL n,e,t,f[N],a[N],Segt[N<<2],p; 8 9 void Update(LL now,LL l,LL r,LL x,LL v) 10 { 11 if (l==r) {Segt[now]+=v; return;} 12 LL mid=(l+r)>>1; 13 if (x<=mid) Update(now<<1,l,mid,x,v); 14 else Update(now<<1|1,mid+1,r,x,v); 15 Segt[now]=min(Segt[now<<1],Segt[now<<1|1]); 16 } 17 18 LL Query(LL now,LL l,LL r,LL l1,LL r1) 19 { 20 if (l>r1 || r<l1) return 1e17; 21 if (l1<=l && r<=r1) return Segt[now]; 22 LL mid=(l+r)>>1; 23 return min(Query(now<<1,l,mid,l1,r1),Query(now<<1|1,mid+1,r,l1,r1)); 24 } 25 26 int main() 27 { 28 scanf("%lld%lld%lld",&n,&e,&t); 29 for (int i=1; i<=n; ++i) 30 scanf("%lld",&a[i]); 31 for (int i=1; i<=n; ++i) 32 { 33 while (p<=i-1 && 2*(a[i]-a[p])>=t) 34 Update(1,1,n,p,-2*a[p]),p++; 35 LL temp=1e17; 36 temp=min(temp,Query(1,1,n,1,p-1)); 37 if (p<=i-1) temp=min(temp,Query(1,1,n,p,i-1)+t-2*a[i]); 38 f[i]=min(temp+3*a[i],f[i-1]+a[i]-a[i-1]+t); 39 Update(1,1,n,i,-a[i-1]+f[i-1]); 40 } 41 printf("%lld\n",f[n]+e-a[n]); 42 } 43 /* 44 #include<iostream> 45 #include<cstdio> 46 using namespace std; 47 48 long long n,e,t,f[10001],a[10001]; 49 50 int main() 51 { 52 scanf("%lld%lld%lld",&n,&e,&t); 53 for (int i=1; i<=n; ++i) 54 scanf("%lld",&a[i]); 55 for (int i=1; i<=n; ++i) 56 { 57 f[i]=f[i-1]+(a[i]-a[i-1])+t; 58 for (int j=i-1; j>=1; --j) 59 f[i]=min(f[i],f[j-1]+(a[i]-a[j])*3+(a[j]-a[j-1])+max(1ll*0,t-2*(a[i]-a[j]))); 60 } 61 printf("%lld\n",f[n]+e-a[n]); 62 } 63 */
E - Shik and Travel
题意：
给定a full binary tree，确定一个叶子排序，最小化从叶子到叶子的路径的最大值。
题解：
心态崩了 我码力怎么越来越弱了啊
可以发现这么一个性质，当我们进入一颗子树u后，我们只有遍历完u的所有叶子才能出来。
对于每一个节点，我们存储所有的二元组$(a,b)_u$，表示可以从点u向下走距离a到达一个叶子，遍历完叶子节点后可以从结束的地方向上走b距离回到点u。
可以发现对于每一个点我们存储的二元组$(a,b)$必须是a单增b单减的，其他的都是无用状态。
考虑如何用点u的两个儿子s0,s1合并出u应该有的二元组呢？先二分一个答案mid，然后就成了判断可行性了。
假设s0有一个二元组$(a1,b1)$,s1有一个二元组$(a2,b2)$
显然新二元组走法为$a1->b1$ , $b1->a2$ , $a2->b2$,只要保证$b1+a2+dis(s0,u)+dis(s1,u)<=ans$，那么这两个二元组就可以合并成$(a1,b2)$
枚举子树大小小的那个儿子的二元组，在另一个儿子里二分找a2最大的且满足可以和b1合并的，因为a2尽量大那么b2就是尽量小的了。
感觉写代码的时候脑子有点乱写的非常鬼畜……
1 #include<bits/stdc++.h> 2 #define N (140000) 3 #define LL long long 4 using namespace std; 5 6 struct Node 7 { 8 LL a,b; 9 bool operator < (const Node A) const 10 { 11 if (a!=A.a) return a<A.a; 12 return b<A.b; 13 } 14 }s; 15 struct Edge{int to,next,len;}edge[N<<1]; 16 LL n,a[N],v[N],size[N],dis[N],l,r,ans,mid; 17 LL a1,b1,a2,b2,L,R,MID,ANS; 18 LL head[N],num_edge; 19 vector<Node>A[N],T; 20 21 void add(LL u,LL v,LL l) 22 { 23 edge[++num_edge].to=v; 24 edge[num_edge].next=head[u]; 25 edge[num_edge].len=l; 26 head[u]=num_edge; 27 } 28 29 void Update(LL x,LL s0,LL s1) 30 { 31 if (size[s0]>size[s1]) swap(s0,s1); 32 LL len0=A[s0].size(), len1=A[s1].size(); 33 34 T.clear(); 35 for (int i=0; i<len0; ++i) 36 { 37 a1=A[s0][i].a; b1=A[s0][i].b; 38 L=0; R=len1-1; ANS=-1; 39 while (L<=R) 40 { 41 MID=(L+R)>>1; 42 a2=A[s1][MID].a; b2=A[s1][MID].b; 43 if (b1+a2+dis[s0]+dis[s1]<=mid) ANS=MID,L=MID+1; 44 else R=MID-1; 45 46 } 47 if (ANS==-1) continue; 48 a2=A[s1][ANS].a; b2=A[s1][ANS].b; 49 s.a=a1+dis[s0]; s.b=b2+dis[s1]; T.push_back(s); 50 s.a=b2+dis[s1]; s.b=a1+dis[s0]; T.push_back(s); 51 } 52 53 sort(T.begin(),T.end()); 54 LL maxn=1ll*N*N, len=T.size(); 55 for (int i=0; i<len; ++i) 56 if (T[i].b<maxn) 57 A[x].push_back(T[i]),maxn=T[i].b; 58 } 59 60 void Dfs(int x,int fa) 61 { 62 size[x]=1; 63 LL s0=0, s1=0; 64 for (int i=head[x]; i; i=edge[i].next) 65 { 66 if (edge[i].to==fa) continue; 67 if (!s0) s0=edge[i].to; 68 else s1=edge[i].to; 69 dis[edge[i].to]=edge[i].len; 70 Dfs(edge[i].to,x); 71 size[x]+=size[edge[i].to]; 72 } 73 if (!s0) s.a=0, s.b=0,A[x].push_back(s); 74 else Update(x,s0,s1); 75 } 76 77 bool check(LL x) 78 { 79 for (int i=1; i<=n; ++i) A[i].clear(); 80 Dfs(1,0); return A[1].size()>0; 81 } 82 83 int main() 84 { 85 scanf("%lld",&n); 86 for (int i=1; i<=n-1; ++i) 87 { 88 scanf("%lld%lld",&a[i],&v[i]); 89 add(i+1,a[i],v[i]); add(a[i],i+1,v[i]); 90 } 91 92 l=0, r=1ll*N*N, ans=-1; 93 while (l<=r) 94 { 95 mid=(l+r)>>1; 96 if (check(mid)) ans=mid,r=mid-1; 97 else l=mid+1; 98 } 99 printf("%lld\n",ans); 100 }
F - Shik and Copying String
题意：
求字符串S->T经过几次变化得到。字母改变的方式有两种:S[i][j] = S[i-1][j] || S[i][j-1]
题解：
这题为什么这么神仙啊QwQ
变换的过程发现可以用折线覆盖起来（如下图S=aabab，T=aaaba，需要3步）
而对于一段连续相互影响的折线，他们需要的操作步数就是折线的个数（如上图就有三条折线相互影响）
具体实现就是开个双端队列……然后T串每一个连续相同的串对应一条折线……感觉说也说不明白还是直接上代码吧……QAQ 反正我就是对着题解的代码感悟的
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N (1000009) 5 using namespace std; 6 7 int n,p,ans,head=1,tail,q[N]; 8 char s[N],t[N]; 9 10 int main() 11 { 12 scanf("%d%s%s",&n,s+1,t+1); p=n; 13 if (!strcmp(s+1,t+1)){puts("0"); return 0;} 14 for (int i=n; i>=1; --i) 15 { 16 if (t[i]==t[i-1]) continue; 17 p=min(p,i); 18 while (p && t[i]!=s[p]) p--; 19 if (!p) {puts("-1"); return 0;} 20 21 while (head<=tail) 22 if (q[head]-i-(tail-head+1)>=0) head++; 23 else break; 24 q[++tail]=p; 25 if (i!=p) ans=max(ans,tail-head+1); 26 } 27 printf("%d\n",ans+1); 28 }
AGC 008
C - Tetromino Tiling
题意：
给你一些奇形怪状的面积为4的图形，问最多能拼起来一个面积多大的2*2K的矩形。
题解：
可以发现能用到的只有I,O,J,L四种，而O种是可以直接计算进入答案的。分类讨论一波I,J,L的奇偶性就好了。
1 #include<cstdio> 2 long long ans,a[8]; 3 int main() 4 { 5 for (int i=1; i<=7; ++i) 6 scanf("%lld",&a[i]); 7 ans=a[2]; 8 if (!a[1]||!a[4]||!a[5]) 9 ans+=a[1]/2*2+a[4]/2*2+a[5]/2*2; 10 else 11 { 12 ans+=a[1]+a[4]+a[5]; 13 if (!(a[1]%2==a[4]%2&&a[4]%2==a[5]%2)) ans--; 14 } 15 printf("%lld",ans); 16 }
D - K-th K
题意：
给出数组x[1..n]，要求构造一个长度为n*n的数组a，使得a中1到n每个数恰好出现n次，且对于任意的i，有从左到右第i个出现的i的下标为x[i]。
题解：
将所有数字按出现的位置从小到大sort一下，从前往后把一个数字前面应该有的填了，填完了再从后往前把一个数字后面应该有的填了。看代码应该非常好懂……
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #define N (250000+1000) 5 using namespace std; 6 7 int n,ans[N],now; 8 struct Node{int num,pos;}a[N]; 9 bool cmp(Node a,Node b){return a.pos<b.pos;} 10 11 int main() 12 { 13 scanf("%d",&n); 14 for (int i=1; i<=n; ++i) 15 { 16 scanf("%d",&a[i].pos); 17 a[i].num=i; ans[a[i].pos]=i; 18 } 19 sort(a+1,a+n+1,cmp); 20 21 now=1; 22 for (int i=1; i<=n; ++i) 23 { 24 for (int j=1; j<a[i].num; ++j) 25 { 26 while (ans[now]) now++; 27 ans[now]=a[i].num; 28 } 29 if (now>a[i].pos){puts("No"); return 0;} 30 } 31 now=n*n; 32 for (int i=n; i>=1; --i) 33 { 34 for (int j=1; j<=n-a[i].num; ++j) 35 { 36 while (ans[now]) now--; 37 ans[now]=a[i].num; 38 } 39 if (now<a[i].pos){puts("No"); return 0;} 40 } 41 puts("Yes"); 42 for (int i=1; i<=n*n; ++i) 43 printf("%d ",ans[i]); 44 }
E - Next or Nextnext
没看懂题解所以弃了
F - Black Radius
题意：
给定一棵树，树上有一些点是可选点。可以选中某个可选点和一个半径d，把与选定点距离不超过距离d的点全染黑，求染色方案数。
题解：
瞻仰了一发国家队队长的代码QvQ……果然神仙的代码思路就是清晰
每个点单独考虑，对于以每个点为中心，我们卡出ta的一个d的上界和d的下界，保证对于每个可选点的可选范围取并集能做到方案不重不漏。为了方便计算我们中途不取全集，最后再把全集的答案加上。以下讨论均把当前点当作根来看
（1）如果当前点是可选点，那么下界L肯定是0，即只选中本身。上界R怎么卡呢？可以发现如果以当前点为中心的话，半径肯定是和最深深度次大的子树相关联的。如下图
最深深度次大的子树为2，显然在半径d不超过子树2的最大深度的时候，染色染出来的连通块的直径的中心肯定一直是当前点（根）。
（2）如果当前点不是可选点，我们就要看哪棵子树里面有可选点。如下图
假设1子树有可选点，下界显然是d刚好等于1的最深深度的时候，也就是1子树内的可选点不停向周围扩散，填满了1子树且把2和3子树染了d。上界和上一个情况同理。
所以我们要求的就是一个点的最短路和次短路，DP就可以了。
1 #include<iostream> 2 #include<cstdio> 3 #define N (300000+1000) 4 using namespace std; 5 6 struct Edge{int to,next;}edge[N<<1]; 7 int head[N],num_edge,Maxdep[N],size[N],sum,n,u,v; 8 long long ans; 9 char s[N]; 10 11 void add(int u,int v) 12 { 13 edge[++num_edge].to=v; 14 edge[num_edge].next=head[u]; 15 head[u]=num_edge; 16 } 17 18 void Dfs1(int x,int fa) 19 { 20 size[x]=(s[x]=='1'); 21 for (int i=head[x]; i; i=edge[i].next) 22 if (edge[i].to!=fa) 23 { 24 Dfs1(edge[i].to,x); 25 size[x]+=size[edge[i].to]; 26 Maxdep[x]=max(Maxdep[x],Maxdep[edge[i].to]+1); 27 } 28 } 29 30 void Dfs2(int x,int fa,int maxdis) 31 { 32 int l=(s[x]=='1'?0:(1<<30)), r, max1=maxdis, max2=0; 33 if (sum-size[x]) l=min(l,maxdis); 34 for (int i=head[x]; i; i=edge[i].next) 35 if (edge[i].to!=fa) 36 { 37 if (Maxdep[edge[i].to]+1>max2) max2=Maxdep[edge[i].to]+1; 38 if (max2>max1) swap(max1,max2); 39 if (size[edge[i].to]) l=min(l,Maxdep[edge[i].to]+1); 40 } 41 if (max1==max2) r=max2-1; 42 else if (max1==max2+1) r=max2; 43 else r=max2+1; 44 if (l<=r) ans+=r-l+1; 45 for (int i=head[x]; i; i=edge[i].next) 46 if (edge[i].to!=fa) 47 { 48 if (Maxdep[edge[i].to]+1!=max1) Dfs2(edge[i].to,x,max1+1); 49 else Dfs2(edge[i].to,x,max2+1); 50 } 51 } 52 53 int main() 54 { 55 scanf("%d",&n); 56 for (int i=1; i<n; ++i) 57 { 58 scanf("%d%d",&u,&v); 59 add(u,v); add(v,u); 60 } 61 scanf("%s",s+1); 62 for (int i=1; i<=n; ++i) sum+=(s[i]=='1'); 63 Dfs1(1,0); Dfs2(1,0,0); 64 printf("%lld\n",ans+1); 65 }
AGC 009
C - Division into Two
题意：
给n个不同的数，划分成两个集合X,Y。满足X中的数两两之间的差>=A，y中的数两两之间的差>=B
题解：
假设A<B,那么若存在s[i]-s[i-2]<A则无解，因为i,i-1,i-2的三个数两两间不能放到同一集合里。
设f[i]表示第i个元素放到Y集合的方案数。若f[j]可以转移到f[i]，则一定满足s[i]-s[j]>=B 且[j+1,i-1]的数可以一起被放到X集合里。
也就是说可转移的范围是f[?,j]。而f[?,j]这个区间是随着i右移而右移的，所以用一个变量sum时刻存储一下这段区间的和就好了。
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #define N (100009) 5 #define MOD (1000000007); 6 using namespace std; 7 long long n,a,b,r,sum,s[N],f[N]; 8 int main() 9 { 10 scanf("%lld%lld%lld",&n,&a,&b); 11 if (a>b) swap(a,b); 12 for (int i=1; i<=n; ++i) 13 scanf("%lld",&s[i]); 14 for (int i=3; i<=n; ++i) 15 if (s[i]-s[i-2]<a){puts("0"); return 0;} 16 f[0]=1; s[0]=-b; sum=0; r=0; 17 for (int i=1; i<=n; ++i) 18 { 19 while (r<i && s[i]-s[r]>=b) sum=(sum+f[r++])%MOD; 20 f[i]=sum; 21 if (s[i]-s[i-1]<a) sum=0, r=i-1; 22 } 23 for (int i=r; i<=n; ++i) 24 sum=(sum+f[i])%MOD; 25 printf("%lld\n",sum); 26 }
D - Uninity
题意：
给定一棵树，让你自己确定点分治的分治中心，问点分树的最小树高是多少
题解：
这题真是神仙题，托爷比赛的时候都没A掉……
首先由点分治的性质可以发现，若我们给每个点编号为它在点分树上的深度，那么若x点编号为k，y点编号为k，则x->y的路径上肯定存在一个点的编号小于k。这个很容易由点分治的过程得出。为了方便计算，后面将叶子标号为0，根标号为最大深度。
下面考虑原树，对于一个点来说，若它的子树里面有两个编号为k的，那么这个点的编号就至少为k+1。又因为如果按点分治的流程走，标号最大为log，所以我们可以开一个数组f[x][i]来存储子树x内有多少个编号为i的点。
若一个点标号为lim，那么它的子树里所有标号小于lim的点都合法了，因为这些点一定是合法的了。看代码可能更好懂……？
1 #include<iostream> 2 #include<cstdio> 3 #define N (100000+1000) 4 using namespace std; 5 6 struct Edge{int to,next;}edge[N<<1]; 7 int n,u,v,dep[N],f[N][22],ans; 8 int head[N],num_edge; 9 10 void add(int u,int v) 11 { 12 edge[++num_edge].to=v; 13 edge[num_edge].next=head[u]; 14 head[u]=num_edge; 15 } 16 17 void Dfs(int x,int fa) 18 { 19 for (int i=head[x]; i; i=edge[i].next) 20 if (edge[i].to!=fa) 21 { 22 Dfs(edge[i].to,x); 23 for (int j=0; j<=20; ++j) 24 f[x][j]+=f[edge[i].to][j]; 25 } 26 int lim=0; 27 for (int i=20; i>=0; --i) 28 if (f[x][i]>1) {lim=i+1; break;} 29 while (f[x][lim]) lim++; 30 dep[x]=lim; 31 for (int i=0; i<=lim-1; ++i) f[x][i]=0; 32 f[x][lim]++; ans=max(ans,dep[x]); 33 } 34 35 int main() 36 { 37 scanf("%d",&n); 38 for (int i=1; i<=n-1; ++i) 39 { 40 scanf("%d%d",&u,&v); 41 add(u,v); add(v,u); 42 } 43 Dfs(1,0); 44 printf("%d\n",ans); 45 }
E - Eternal Average
题意：
给你n个0，m个1，和一个k。每次操作你选择k个数，擦去这k个数并加入他们的平均数（1个），问最后会有多少种不同的实数。
题解：
我也就指望着题解苟活一下了……
首先转换一下模型，假设现在有一个k叉树，有n+m个叶子，n个叶子为0，m个叶子为1。非叶子节点的权值等于它儿子权值的平均数。设n个为0的点深度分别为$a_1…a_n$，m个为1的点深度为$b_1…b_m$，那么可以列出式子$\sum_{i=1}^nk^{-ai}+\sum_{i=1}^mk^{-bi}=1$
这是显然正确的，因为这个式子就相当于把0看成1求平均数，最后平均数一定为1.
而又显然根节点权值为$ans=\sum_{i=1}^mk^{-bi}$，那么ans可以表示成m个k的多少次幂的和，1-ans可以表示成n个k的多少次幂的和。不妨用k进制来表示ans
设$ans=0.s_1s_2...s_l$
如果不考虑进位的话，所有s的和就是m，那么如果考虑进位呢？
因为是k进制所以只要所有s的和%(k-1)==m%(k-1)就可以了。1-ans同理。
接下来我们只需要构造一个ans，使得ans和1-ans满足条件就好了。
用$f[i][j][0/1]$表示DP到小数点后第i位，前i位的和为j，第i位是否为0。转移比较显然，可以用前缀和优化。
1 #include<iostream> 2 #include<cstdio> 3 #define N (4009) 4 #define MOD (1000000000+7) 5 using namespace std; 6 long long n,m,k,sum[N],f[N][N][2],ans; 7 int main() 8 { 9 scanf("%lld%lld%lld",&n,&m,&k); m--; k--; 10 f[0][0][0]=1; 11 for (int i=1; i<=n+m; ++i) 12 { 13 for (int j=0; j<=n; ++j) 14 sum[j+1]=(sum[j]+f[i-1][j][0]+f[i-1][j][1])%MOD;//因为可能有sum[-1]所以sum数组整体右移了一位 15 for (int j=0; j<=n; ++j) 16 { 17 f[i][j][0]=(sum[j+1]-sum[j]+MOD)%MOD; 18 f[i][j][1]=(sum[j]-sum[max(j-k,0ll)])%MOD; 19 } 20 for (int j=0; j<=n; ++j) 21 if (j%k==n%k && (i*k-j)%k==m%k && i*k-j<=m) 22 ans=(ans+f[i][j][1])%MOD; 23 } 24 printf("%lld\n",(ans+MOD)%MOD); 25 }
AGC 010
C - Cleaning
题意：
给你一颗树，每个节点上有一些石子，你可以选中两个叶子节点，并将路径上的节点的石子数量减一。必须保证选中路径上每个节点都有石子。问是否可以把石子删光。
题解：
一个exciting的贪心应该是……
假设现在只有一个两层的树，包含一个根节点和若干个叶子。要想使叶子节点全部变成0，需要满足：
1、叶子权值和一定≥父亲权值和，因为叶子节点减1的时候父亲节点肯定也要减1.。
2、叶子权值和≤2*父亲权值和，因为叶子节点的权值每次是减2的，而父亲节点是减1。
3、叶子权值的最大值应该小于等于父亲节点。
然后从下往上推就好了……
1 #include<iostream> 2 #include<cstdio> 3 #define N (100009) 4 using namespace std; 5 6 struct Edge{int to,next;}edge[N<<1]; 7 int n,u,v,a[N],Ind[N],flag=true; 8 int head[N],num_edge; 9 10 void add(int u,int v) 11 { 12 edge[++num_edge].to=v; 13 edge[num_edge].next=head[u]; 14 head[u]=num_edge; 15 } 16 17 void Dfs(int x,int fa) 18 { 19 if (!flag) return; 20 int maxn=0; 21 long long sum=0; 22 for (int i=head[x]; i; i=edge[i].next) 23 if (edge[i].to!=fa) 24 { 25 Dfs(edge[i].to,x); 26 sum+=a[edge[i].to]; 27 maxn=max(maxn,a[edge[i].to]); 28 } 29 if (!maxn) return; 30 if (sum<a[x] || sum>2*a[x] || maxn>a[x]) flag=false; 31 a[x]=a[x]*2-sum; 32 } 33 34 int main() 35 { 36 scanf("%d",&n); 37 for (int i=1; i<=n; ++i) 38 scanf("%d",&a[i]); 39 for (int i=1; i<=n-1; ++i) 40 { 41 scanf("%d%d",&u,&v); 42 add(u,v); add(v,u); 43 Ind[u]++; Ind[v]++; 44 } 45 46 if (n==2) 47 { 48 if (a[1]==a[2]) puts("YES"); 49 else puts("NO"); 50 return 0; 51 } 52 53 int root=1; 54 while (Ind[root]<=1) root++; 55 Dfs(root,0); 56 if (a[root]!=0) flag=false; 57 if (flag) puts("YES"); 58 else puts("NO"); 59 }
D - Decrementing
题意：
有一个gcd为1的序列，每次可以选择一个数减1，然后将整个序列除以他们的gcd。先把序列变成全1的获胜。问是先手胜还是后手胜。
题解：
博弈题玩弄人性（雾
显然当序列为x,x,x,x+1的形式的时候，游戏才能一步结束。另外一个人肯定要避免让对手到达这样的局面，当且仅当x=1的时候，游戏才会结束。
1、偶数个数为奇数的时候，先手必胜。因为先手可以通过控制奇偶性不变，来到达上述必胜局面。
2、偶数个数为偶数的时候，先手要尝试改变偶数个数的奇偶性。
（1）显然当奇数个数大于1或者奇数中存在1的时候，无法改变，先手必败。因为要改变个数奇偶性的唯一方法是同除一个偶数，而这显然是不可能的……
（2）如果上面都不满足就递归模拟。每次当前操作者肯定不能选偶数进行操作，不然就会把对方送进必胜局面。每次一方随便选择一个奇数，直到出现前面提到的局面的时候就可以停止递归了。最多递归log层。
1 #include<iostream> 2 #include<cstdio> 3 using namespace std; 4 5 int n,a[100009]; 6 int gcd(int a,int b){return b==0?a:gcd(b,a%b);} 7 8 bool Dfs(int x) 9 { 10 int odd=0,even=0,flag=1; 11 for (int i=1; i<=n; ++i) 12 if (a[i]%2) {a[i]--; break;} 13 int g=a[1]; 14 for (int i=2; i<=n; ++i) g=gcd(g,a[i]); 15 for (int i=1; i<=n; ++i) 16 { 17 a[i]/=g; 18 if (!(a[i]%2)) even++; 19 else odd++; 20 if (a[i]==1) flag=0; 21 } 22 if (even%2) return !x; 23 if (odd!=1 || !flag) return x; 24 return Dfs(x^1); 25 } 26 27 int main() 28 { 29 int odd=0,even=0,flag=1; 30 scanf("%d",&n); 31 for (int i=1; i<=n; ++i) 32 { 33 scanf("%d",&a[i]); 34 if (a[i]%2) odd++; 35 else even++; 36 if (a[i]==1) flag=false; 37 } 38 if (even%2) puts("First"); 39 else if (odd!=1 || !flag) puts("Second"); 40 else if (Dfs(1)) puts("First"); 41 else puts("Second"); 42 }
