NOIP 2012 提高组 复赛 day2 mod 同余方程
1.一看题目,还好数论里瞄过,没被数学符号给吓到,要不要将数论书再翻一翻,还是将该题AC了再说。
2.数学有储备,见了该题就不慌了。
3.看了该题,心里明白,一不小心就容易数据溢出,为了争取尽可能多的分,long long是免不了了。心里有个小九九,用C来写,写好后,正好去linux下用Arbiter进行测试,看看对%lld的支持。
4.简单做了运算int 最大值2*10^9 long long 最大值4*10^9*10^9,运算过程中随时注意取模,应该不会溢出。
5.编写代码,提交50分,剩下5个数据,全报TLE,超时了。
6.突然间想到一个办法,思路来自(10*2+1)/3=7,提交60分,剩下4个数据,全报TLE,超时了。
7.网络中搜索了一下,说是扩展欧几里德算法,好吧,开始学习。
8.翻书的过程中发现,关于信息学的书涉及的数学,比纯数学的书涉及的数学,讲解得更具体,更接近计算机的具体实现。故,翻书的顺序,先翻信息学的书,要扩充,完毕知识时,再翻纯数学的书。
9.学习一天未果,好吧,上代码学习,并进行相应跟踪。
10.稀里糊涂的编码,用的是扩展欧几里德算法,测试样例发现是-3,果断加上10,简单修改代码,提交AC。
11.还是要弄明白原理。再学习学习。
12.关于扩展欧几里德算法,这两篇文章介绍得不错,建议先阅读第一篇,之后第二篇。
http://www.cnblogs.com/ka200812/archive/2011/09/02/2164404.html
http://www.cnblogs.com/comeon4mydream/archive/2011/07/18/2109060.html
13.2011年时POJ1061青蛙的约会,比较流行,也即扩展欧几里德算法,2012年进行改编成为NOIP提高组复赛day2第一题。
附上AC代码,编译环境Dev-C++4.9.9.2
//2012 mod3
#include <stdio.h>
void gcd(int a,int b,int &d,int &x,int &y){
if(!b){
d=a;
x=1;
y=0;
}else{
gcd(b,a%b,d,y,x);
y-=x*(a/b);
}
}
int main(){
int a,b,d,x,y;
scanf("%d%d",&a,&b);
gcd(a,b,d,x,y);
while(x<0){
x+=b;
}
printf("%d\n",x);
return 0;
}
附上60分代码,编译环境Dev-C++4.9.9.2
//2012 同余方程 mod2
#include <stdio.h>
int main(){
long long a,b,i,j;
scanf("%lld%lld",&a,&b);
i=1;
while(1){
if((b*i+1)%a==0)
break;
i++;
}
printf("%lld\n",(b*i+1)/a);
return 0;
}
附上50分代码,编译环境Dev-C++4.9.9.2
//2012 同余方程
#include <stdio.h>
int main(){
long long a,b,i,t,j;
scanf("%lld%lld",&a,&b);
i=1;
t=a;
t%=b;
while(1){
j=i;
j%=b;
if((t*j)%b==1)
break;
i++;
}
printf("%lld\n",i);
return 0;
}
附上poj1061青蛙的约会AC代码,编译环境Dev-C++4.9.9.2(很多小的细节)
//poj1061 青蛙的约会
#include <stdio.h>
void gcd(long long a,long long b,long long &d,long long &x, long long &y){
if(!b){
d=a;
x=1;
y=0;
}else{
gcd(b,a%b,d,y,x);
y-=a/b*x;
}
}
int main(){
long long a,b,d,x,y,b2;
long long t1,t2,t3,t4,t5;
scanf("%lld%lld%lld%lld%lld",&t1,&t2,&t3,&t4,&t5);
a=t3-t4;
b=t5*-1;
gcd(a,b,d,x,y);
if((t2-t1)%d!=0){
printf("Impossible\n");
}else{
b2=b/d;
x=((t2-t1)/d*x%b2+b2)%b2;
printf("%lld\n",x);
}
return 0;
}