链接:http://acm.hdu.edu.cn/showproblem.php?pid=1098
给出k值,求最小正整数a使对于任意整数x满足65|f(x),f(x)=5*x^13+13*x^5+k*a*x
介绍下我第一次做的方法,网上大部分都是数学归纳得出的算法,代码看起来要简单些
首先一个多项式是否能被一个整数整除,只要将[1,最高次幂+1]依次代入x检测是否整除,若全部整除,则x为任意整数都能整除,否则不能整除(证明见刘汝佳《算法竞赛》数论例8)
于是依次可以得出等式:
(5*x^13+13*x^5+k*a*x) mod 65=0
(k*a*x) mod 65=65-(5*x^13+13*x^5)mod 65
因为上式a才是要求的未知数,为了下面表示方便以及和代码对应令 a=kx, xx=a
得:(a*xx) mod 65=65 -(5*x^13+13*x^5)mod 65
令 c = 65 - (5*x^13+13*x^5)mod 65
则得同余方程 a*xx = c(mod 65)
引入变量y得二元不定式 : a*xx+65*y=c
再将[1,14]依次代入x,求得c,a,解该方程最小正整数解xx,取所有解中的最大值便是本题的答案
注意计算c的时候x^13比较大,用int不行,但用long或long long hdu又过不了,可以只计算(5*x+13*x)mod 65(原理还在研究中)
#include<stdio.h>
#include<string.h>
void ex_gcd(int a,int b,int &d,int &x,int &y)
{
if(b==0) {d=a,x=1,y=0;}
else {ex_gcd(b,a%b,d,y,x);y-=x*(a/b);}
}
int main()
{
int a,c,i;
int x,y,d;
int k;
while(scanf("%d",&k)!=EOF)
{
int flag=0;
int ans=0;
for(i=1;i<=14;i++)
{
a=(int)k*i;
c=65-((5*i+13*i)%65);
ex_gcd(a,(int)65,d,x,y);
if(c%d)
{
flag=1;
printf("no\n");
break;
}
x=x*(c/d);
int t=65/d;
x=(x%t+t)%t;
if(x>ans) ans=(int)x;
}
if(!flag) printf("%d\n",ans);
}
return 0;
}