[poi2011]bzoj 2277 —— strongbox·[洛谷3518]

·问题描述·

　　有一个密码箱，0到n-1中的某些数是它的密码。且满足：如果a和b都是它的密码，那么（a+b）%n也是它的密码。某人试了k次密码，前k-1次都失败了，最后一次成功。

　　问：该密码箱最多有多少个密码？

·输入格式·

　　输入第一行两个整数分别表示n，k。

　　第二行为k个用空格隔开的非负整数，表示每次试的密码。（数据保证存在合法解）

·输出格式·

　　输出一行一个数，表示结果。

·输入样例·

42 5

28 31 10 38 24

·输出样例·

14

·数据范围·

对于100%的数据：1<=k<=250000，k<=n<=10^14。

 

Solution：

　　本题考察数学。由题意可知，若x为密码则（x+x）%n为密码，则p*x%n（0

　　由引理：a*x+b*y=c（均为整数），有整数解的充要条件是gcd(a,b)|c。证明很容易：不妨设a=p*gcd(a,b),b=q*gcd(a,b) --> a*x+b*y=(p+q)*gcd(a,b)=c，显然要有整数解，则gcd(a,b)|c。

　　回到本题的条件：p*x-n*c=q。有解则必定满足gcd(x,n)|q，所以必定有p*x-n*c=gcd(x,n)成立，等价于p*x%n=gcd(x,n)，则gcd(x,n)一定为一个密码。类似的，对于不同的密码x和y，存在（p*x+q*y）%n为密码，由引理必定存在p*x+q*y=gcd(x,y)，与单个x是密码同理gcd(x,y)一定是密码。

　　而要使得密码最多，由x是密码则p*x%n（0为密码可知，当x为最小时，密码最多有n/x个。

　　具体实现时，我们先求出a[k]=gcd(a[k],n)，再使a[i]=gcd(a[i],a[k])，然后从新的a[k]中删去所有是a[i]因子的因子，最后输出答案就是n除以没被删的最小的因子。

代码：

 

 1 /*数学一本通上的例题——by 520*/
 2 #include
 3 #define il inline
 4 #define ll long long
 5 using namespace std;
 6 ll n,k,tot,a[250005],q[250000],p[250000],cnt=1;
 7 il ll gi(){
 8     ll a=0;char x=getchar();bool f=0;
 9     while((x<'0'||x>'9')&&x!='-')x=getchar();
10     if(x=='-')x=getchar(),f=1;
11     while(x>='0'&&x<='9')a=a*10+x-48,x=getchar();
12     return f?-a:a;
13 }
14 il ll gcd(ll a,ll b){
     return b?gcd(b,a%b):a;}
15 int main()
16 {
17     freopen("strongbox.in","r",stdin);
18     freopen("strongbox.out","w",stdout);
19     n=gi(),k=gi();
20     for(int i=1;i<=k;i++)a[i]=gi();
21     a[k]=gcd(a[k],n);
22     for(int i=1;igcd(a[k],a[i]);
23     for(ll i=1;i*i<=a[k];i++)
24         if(a[k]%i==0){
25             q[++tot]=i;
26             if(i*i!=a[k])q[++tot]=a[k]/i;
27         }
28     sort(q+1,q+tot+1);
29     for(int i=1;i1,q+tot+1,a[i])-q]=1;
30     for(int i=1;i<=tot;i++)
31         if(p[i])
32             for(int j=1;j)
33                 if(q[i]%q[j]==0)p[j]=1;
34     while(p[cnt])cnt++;
35     cout<q[cnt];
36     return 0;
37 }

 

转载于:https://www.cnblogs.com/five20/p/8414577.html