hdu1098 -数论

最新推荐文章于 2018-09-05 21:37:00 发布

樱木花涛

最新推荐文章于 2018-09-05 21:37:00 发布

阅读量501

点赞数

CC 4.0 BY-SA版权

分类专栏：数学

本文链接：https://blog.youkuaiyun.com/the_ilemon/article/details/40392243

数学专栏收录该内容

4 篇文章

订阅专栏

本文探讨了一个关于整数运算的独特数学问题，目标是在给定的负整数k下，找到使得任意整数x时，5*x^13 + 13*x^5 + k*a*x能被65整除的最小正整数a。通过数学归纳法，证明了当a满足特定条件时，该表达式对于所有整数x都具备被65整除的性质。文章详细解释了如何使用扩展欧几里得算法来求解最小的a值，并提供了实例输入输出，以便读者理解和实践这一过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Ignatius's puzzle

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 6563 Accepted Submission(s): 4544

Problem Description

Ignatius is poor at math,he falls across a puzzle problem,so he has no choice but to appeal to Eddy. this problem describes that:f(x)=5*x^13+13*x^5+k*a*x,input a nonegative integer k(k<10000),to find the minimal nonegative integer a,make the arbitrary integer x ,65|f(x)if
no exists that a,then print "no".

Input

The input contains several test cases. Each test case consists of a nonegative integer k, More details in the Sample Input.

Output

The output contains a string "no",if you can't find a,or you should output a line contains the a.More details in the Sample Output.

Sample Input

Sample Output

Author

eddy

思路：题目的关键是函数式f(x)=5*x^13+13*x^5+k*a*x; 用数学归纳法证明：x取任何值都需要能被65整除..

所以我们只需找到f(1)成立的a,并在假设f(x)成立的基础上,
证明f(x+1)也成立即可。
那么把f(x+1)展开，则f(x+1 ) = f (x) + 5*( (13 1 ) x^12 ...... .....+(13 13) x^0 )+ 13*( (5 1 )x^4+...........+ ( 5 5 )x^0 )+k*a;

很容易证明，除了5*(13 13) x^0 、13*( 5 5 )x^0 和k*a三项以外,其余各项都能被65整除.
那么也只要求出18+k*a能被65整除就可以了.

（参考原文）

即求18+k*a能被65整除的最小正数a，即求解18+k*a=65*n，利用扩展欧几里得算法求解即可

#include<iostream>
#include<cmath>
using namespace std;
void extend_gcd(int x,int y,int &d,int &a,int &b){
	if(y==0){
		d=x;
		a=1;
		b=0;
	}
	else{
		extend_gcd(y,x%y,d,b,a);
		b-=a*(x/y);
//		a-=b*(x/y);
	}
}
int main(){
	int i,j,k,m,n,d,a;
	while(cin>>k){
		 extend_gcd(65,k,d,n,a);
//		 cout<<d<<" "<<a<<" "<<n<<endl;
		 if(18%d)cout<<"no\n";
		 else{
		 	a*=18/d;
		 	if(a>0) a=a-(int)ceil((double)(a)/(65/d))*(65/d);
		 	else a=a+(-a)/(65/d)*(65/d);
		 	cout<<-a<<endl;
		 }
	}
	return 0;
}