数学/扩展欧几里得/NOI2002 荒岛野人 savage

题目：

　　岛上有排列成环行的M个山洞。这些山洞顺时针编号为1,2,…,M。岛上住着N个野人，一开始依次住在山洞C1,C2,…,CN中，以后每年，第i个野人会沿顺时针向前走Pi个洞住下来。每个野人i有一个寿命值Li，即生存的年数。奇怪的是，虽然野人有很多，但没有任何两个野人在有生之年处在同一个山洞中，使得小岛一直保持和平与宁静，这让科学家们很是惊奇。他们想知道，至少有多少个山洞，才能维持岛上的和平。

分析：

　　首先分析两个野人i,j的情况。如果i,j在x年相遇并且在第x年时俩人都活着，则可以得到一个同余方程：ci+x*pi≡cj+x*pj (mod M).现在我们要求不能使野人们相遇，则需要让这个同余方程无解，或解出的最小的x比两个人中任何一人的寿命长。

上述是两个野人之间的分析，但对于n个野人呢，由于n<=15，我们完全可以双重循环两两判断比较。M的初始值为max{c[i]}（易证，若山洞个数比初始位置还少,是不可行的），每次两两判断时，若有一对能相遇，则这个M不可行，需要M++,直到找到一个M，使每两个野人都不能相遇，则这个M为最终答案。

 

具体实现：

同余方程：ci+x*pi≡cj+x*pj (mod M)

移项可得：(pi-pj)*x≡cj-ci (mod M)

即求满足(pi-pj)*x-M*y=cj-ci 的最小的正整数x

其中在每一步,pi-pj已知，M已知，cj-ci已知，方程即转化为ax+by=c，很眼熟吧？用扩展欧几里得就能求出一组可行的(x,y)的解，但是此时的x并不是最小的正整数解，还需借助通解式x=x0+kb来求出最小正整数解，具体用扩展欧几里得求解二元一次不定方程的解法在我另一随笔上有详细介绍，请移步http://www.cnblogs.com/Rinyo/archive/2012/11/25/2787419.html

 1 /*

 2     PROBLEM:NOI2002 savage

 3     AUTHER:Rinyo

 4     MEMO:扩展欧几里得

 5 */

 6 #include<cstdio>

 7 #include<cstdlib>

 8 int n,m;

 9 int c[20],p[20],l[20];

10 inline int cmin(const int &x,const int &y) {return x<y?x:y;}

11 inline int cmax(const int &x,const int &y) {return x>y?x:y;}

12 int exgcd(int a,int b,int& x,int& y)

13 {

14     if (b==0) {x=1;y=0;return a;}

15     int d=exgcd(b,a%b,x,y);

16     int t=x;x=y;y=t-(a/b)*y;

17     return d;

18 }

19 int main()

20 {

21     scanf("%d",&n);

22     for (int i=1;i<=n;i++) {scanf("%d%d%d",&c[i],&p[i],&l[i]);m=cmax(m,c[i]);}

23     for (int k=m;;k++)

24     {

25         bool flag=true;

26         for (int i=1;i<=n-1;i++)

27         {

28             for (int j=i+1;j<=n;j++)

29             {

30                 int x,y;

31                 int pp=p[i]-p[j],cc=c[j]-c[i];

32                 int d=exgcd(pp,k,x,y);

33                 if (cc%d) continue;

34                 int now=x*cc/d;

35                 now=now % (k/d);

36                 if (now<0) now+=abs(k/d);

37                 if (now<=cmin(l[i],l[j])) {flag=false;break;}

38             }

39             if (!flag) break;

40         }

41         if (flag) {printf("%d",k);break;}

42     }

43     return 0;

44 }