ACM今日学习总结(2015.07.28)

一、欧几里得定理相关应用
int gcd(int a,int b){return !b ? a : gcd(b,a % b);}
定理1 gcd(a,b)是ax+by的线性组合的最小正整数，x, y∈z;
定理2 如果ax+by=c，x, y∈z；则c% gcd(a,b) ==0;
定理3
如果a，b是互质的正整数，c是整数，且方程 ax+by=c有一组整数解x0，y0则此方程的一切整数解可以表示为 x=x0+bt;y=y0-at;t∈z;
扩展欧几里得定理
对于两个不全为0的整数a、b，必存在一组解x,y，
使得ax+by==gcd(a,b)
由定理：ax+by==gcd(a,b)，若b==0，则gcd(a,0)==a。
原式变为ax+by==a –> x==1,y==0。
若x,y表示第一次递归时的值，x1,y1表示第二次递归时的值。
则gcd(a,b)==gcd(b,a%b)，有ax+by==b*x1+(a%b)*y1。
b*x1+(a%b)y1==b*x1+(a-(a/b)*b)*y1==a*y1+b(x1-(a/b)*y1)，
最终得到ax+by==a*y1+b*(x1-(a/b)*y1)
即上一深度的x等于下一深度的y1
上一深度的y等于下一深度的x1-(a/b)*y1。
对于一般的ax+by==c，x1=x*(c/gcd(a,b)),y1=y*(c/gcd(a,b))
需要注意，上面推导时用的除法都是整型除法

int gcd(int a,int b)
{
      int t,d;
      if(b==0)
      {
           x=1;
           y=0; 
           return a;
       }
      d=gcd(b,a%b);
      t=x;
      x=y;
      y=t-(a/b)*y; 
      return d;
}  

上面的程序中,x和y我是用全局变量保存的
再深入一点，就解出这么一组解其实一般来说是解决不了什么问题的。没有哪个ACM的题这么简单吧。
比如我们现在要得到所有的解，那么这所有的解究竟是什么呢？
直接说吧,假设d=gcd(a,b). 那么x=x0+b/d*t; y=y0-a/d*t;其中t为任意常整数。

二、初学者常见问题

1）Main函数必须返回int类型(正式比赛)
2）不要在for语句中定义类型
3）__int64不支持，可以用long long代替
4）使用了汉语的标点符号
5）itoa不是ansi函数
能将整数转换为字符串
而且与ANSI标准兼容的方法是使用sprintf()函数
int num = 100;
char str[25];
sprintf(str, ” %d” , num);
6）另外，拷贝程序容易产生错误