代码源每日一题div1 最大公约数

Daimayuan Online Judge

思路：数学枚举+环上问题
如果将环切成$k$段，那么有$x=gcd(d_1,d_2,\cdot \cdot \cdot ,d_k)$
我们可以发现$d_i=x\ast k_i$，因此$sum=\sum d_i=x\ast \sum k_i$，我们发现$x$是数组和的一个因子，由于$sum\le 10^{11}$,计算得到它的因子数大约$4000$，因此我们可以考虑枚举$x$
接下来我们考虑对于每个因子$x$，如果前$i$个数的和$su{m}_i$和前$j$个数的和$su{m}_j$满足：$su{m}_i\%x=su{m}_j\%x\Rightarrow (su{m}_j-su{m}_i)\%x=0$，因此区间$[i+1,j]$的和是$x$的倍数，那么由于是个环，我们可以知道剩余的数的和也是$x$的倍数
那么如果对于某个因子$x$，我们要得到尽可能多的区间个数$tot$，那显然我们需要统计前缀和$su{m}_i\%x$等于哪个数时的个数最多，最多的个数即为我们对于这个因数$x$所能划分的最多区间个数$tot$，显然区间个数小于$tot$时这个因数$x$都是符合的
那我们要求每个$k$对应的最大公约数，只需要求个后缀最大即可

#include <bits/stdc++.h>
using namespace std;
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
#define endl '\n';
typedef long long ll;
typedef pair<int,int> PII;
const int N=1e6+10,mod=1e9+7;
ll n,a[N],sum[N],ans[N];
void solve(ll x){
	ll tt=0;
	map<ll,ll>mp;
	for(int i=1;i<=n;++i){
		mp[sum[i]%x]++;
		tt=max(tt,mp[sum[i]%x]);
	}
	ans[tt]=max(ans[tt],x);
}
int main(){
	IOS;
	cin>>n;
	for(int i=1;i<=n;i++) cin>>a[i];
	for(int i=1;i<=n;i++) sum[i]=sum[i-1]+a[i];
	int len=sqrt(sum[n]);
	for(int i=1;i<=len;i++){
		if(sum[n]%i==0){
			solve(i);
			solve(sum[n]/i);
		}
	}
	for(int i=n-1;i>0;i--) ans[i]=max(ans[i],ans[i+1]);
	for(int i=1;i<=n;++i) cout<<ans[i]<<endl;
	return 0;
}