luoguP3518 [POI2011]SEJ-Strongbox

Description

题意就是给出一个密码盒子，这个盒子有密码，密码小于n。
如果x、y是密码，那么(x+y)%n是密码，x可以与y相等
某个人试了k次，前k-1次失败，最后一次成功
问这个盒子最多有多少个密码

sample

42 5
28 31 10 38 24

14

Solution

结论一：如果x为密码，则gcd(x,n)也为密码
因为对于方程x*k-n*c=gcd(x,n)一定有解
结论二：如果x,y为密码，则gcd(x,y)也为密码
因为对于方程ax+by=gcd(x,y)一定有解
ax+by=ax+by+pnx+qny (mod n)
(a+pn)x+(b+qn)y=ax+by (mod n)
可使a+pn>0,b+qn>0

在密码集合中，设x为所有元素的gcd，则x一定为集合中最小的元素
而剩余元素一定为2x,3x,4x……

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int N=2500005;
ll n,i,j,cnt,k;
ll a[N],q[N],f[N];
ll read(){
    ll sum=0;
    char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9'){
        sum=sum*10+c-'0';
        c=getchar();}
    return sum;
}
ll gcd(ll x,ll y){
    ll z;
    while (x%y){
        z=x%y; x=y; y=z;}
    return y;
}
int main(){
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
    n=read(),k=read();
    fo(i,1,k) a[i]=read();
    a[k]=gcd(a[k],n);
    fo(i,1,k-1) a[i]=gcd(a[i],a[k]);
    for(i=1;i*i<=a[k];i++){
        if (a[k]%i==0) {
            q[++cnt]=i;
            if (i*i!=a[k]) q[++cnt]=a[k]/i;}
    }
    sort(q+1,q+cnt+1);
    fo(i,1,k-1)
      f[lower_bound(q+1,q+cnt+1,a[i])-q]=1;
    fo(i,1,cnt)
      if (f[i])
        fo(j,1,i-1)
          if (q[i]%q[j]==0) f[j]=1;
    i=1;
    while (f[i]) i++;
    printf("%lld",n/q[i]);
}