首页 > 代码库 > Codeforces 354C Vasya and Beautiful Arrays[dp+暴力]
Codeforces 354C Vasya and Beautiful Arrays[dp+暴力]
题意:
给出n个整数,对每个整数可以减去0-k的任意一个数
求这样操作后,n个数的最大GCD是多少
分析:
我们首先可以知道n个整数中最小的数是多少
而且,最终的答案肯定不大于这个数
这个n个整数中最小的数是答案的上限
然后对于答案的下限
可以肯定的是
1肯定是答案的下限
2呢?3呢?为什么1一定是
其实,0-k+1,都可以作为答案
为什么?
可以把k想象成一个剪刀
对k+1来说,任何数都可以剪掉0-k变成k+1的倍数(任何数模k+1的结果都是0-k)
所以0-k也可以,综上0-k+1都可以,所以答案的下限是k+1
答案处于[k+1,minn]
如果minn<k+1,那么就是minn
如果minn>k+1,则我们枚举这其中所有的数
看最大的能使得所有数列中的数经过操作之后GCD(all)确实能等于这个数的数是多少
也就是验证这个数的可行性,找出可行的数中最大的那个
怎么验证呢?给出一个数m,怎么才能知道可行不可行呢?
一步步来推理
每个数模m,结果都是0--m-1
对于0--m-1,我们的剪刀最多只能剪掉k
所以如果数列中某个数模m的结果大于k,则m不可行
那么我们是不是能这样做
对数列中的每一个数都模m呢
估计一下复杂度,枚举是100W,最多有30万个数,则100w*30w,T!
转化一下这个问题
我们是对每个数尝试一下它是不是处于可行区间
反过来,如果可行区间有这个数的话,是一样的
如果所有可行区间中存在的数的和是n,则这个m是可行的
可行区间是什么?
一个数模m即是x =m*d+x%m
x%m小于等于k,也就是x处于[m*d,m*d+k]
那么所有的可行区间就是
[m*1,m*1+k],[m*2,m*2+k],[m*3,m*3+k],[m*4,m*4+k]...
有logn个这样的区间
算下复杂度,n*logn
这些区间内出现的数的总和我们可以累加前缀和之差得到
#include <cstdio> #include <cstdlib> #include <cstring> #include <cmath> #include <algorithm> #include <iostream> using namespace std; const int NN=555555; int vistmp[NN*2]; int vis[NN*2]; int f[NN]; int n,k; int maxn; bool check(int s){ int sum=0; for(int i=1;i*s<=maxn;i++){ sum += vis[min(maxn,i*s+k)] - vis[i*s-1];//注意细节 if(s==7){ } } return sum==n; } int main(){ #ifndef ONLINE_JUDGE freopen("G:/in.txt","r",stdin); //freopen("G:/myout.txt","w",stdout); #endif cin>>n>>k; int minn=1<<30; for(int i=1;i<=n;i++){ cin>>f[i]; minn=min(minn,f[i]); maxn=max(maxn,f[i]); vistmp[f[i]]++; } vis[minn]=vistmp[minn]; for(int i=minn+1;i<=maxn;i++) vis[i]=vis[i-1]+vistmp[i]; if(minn<=k+1){ cout<<minn<<endl; return 0; } for(int i=minn;i>=k+1;i--){ if(check(i)){ cout<<i<<endl; return 0; } } }
Codeforces 354C Vasya and Beautiful Arrays[dp+暴力]