[Noi2002]Savage 题解

[Noi2002]Savage

时间限制: 5 Sec  内存限制: 64 MB

题目描述

输入

第1行为一个整数N(1<=N<=15)，即野人的数目。

第2行到第N+1每行为三个整数Ci, Pi, Li表示每个野人所住的初始洞穴编号，每年走过的洞穴数及寿命值。

(1<=Ci,Pi<=100, 0<=Li<=10^6 )

输出

仅包含一个数M，即最少可能的山洞数。输入数据保证有解，且M不大于10^6。

样例输入

3
1     3      4
2     7      3
3     2      1

样例输出

6
//该样例对应于题目描述中的例子。
　　这道题明眼人都能看得出来是扩展欧几里得，然而怎么搞就是一个问题了，如果对扩展欧几里得不是太熟可以先做一下 青蛙的约会 裸题，但要注意一下你的模板必须正确，否则像本博主这样的蒟蒻调半天才发现模板有错就崩了。（本博主扩展欧几里得的板子来自Q某犇，据Q某犇说他也是找了好久才发现神利.代目学长写正确模板，其余好多都是错误的）
　　首先我们完全可以把青蛙那道题的主要代码都搬过来，由于这道题并不符合单调性，想二分的同学就扑街了，我也是其中之一。
　　由于n很小m最大不过1000,000我们从小到大挨个枚举就好了，时间复杂度最坏就是n^2*log（max(c[i])）*1000000，反正时间有5秒，而且这只是最坏复杂度，所以貌似是可以的。然后就简单了，两个野人相遇的条件是他们能相遇且两人都存活，因此我们只要算出他们两个相遇的时间是否比他们中寿命最短的那个长就行了。
　　

 #include<iostream>
 #include<cstdlib>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<map>
 #include<queue>
 #include<string>
 #include<cmath>
 using namespace std;
 int n,c[],p[],l[];
 int exgcd(int a,int b,int &x,int &y){
     if(b==)
     {
         x=;
         y=;
         return a;
     }
     int t=exgcd(b,a%b,x,y);
     int o=x;
     x=y;
     y=o-a/b*y;
     return t;
 }
 bool check2(int a,int b,int L){
     int x,y;
     int gcd=exgcd(p[a]-p[b],L,x,y);
     if(((c[b]-c[a])%gcd)!=)
         return ;
     int aa=p[a]-p[b],bb=L;
     aa/=gcd,bb/=gcd;
     exgcd(aa,bb,x,y);
     bb=abs(bb);
     x=x*(c[b]-c[a])/gcd;
     x%=bb;
     if(x<) x+=bb;
     if(x>min(l[a],l[b])) return ;
     return ;
 }
 bool check1(int L){
     for(int i=;i<=n;i++)
     {
         for(int j=i+;j<=n;j++)
         {
             if(!check2(i,j,L))
             {
                 return ;
             }
         }
     }
     return ;
 }
 int main(){
 //    freopen("savage.in","r",stdin);
 //    freopen("savage.out","w",stdout);
     scanf("%d",&n);
     int st=;
     for(int i=;i<=n;i++)
     {
         scanf("%d%d%d",&c[i],&p[i],&l[i]);
         st=max(st,c[i]);
     }
     int ans=;
     for(int i=st;;i++)
     {
         if(check1(i))
         {
             ans=i;
             break;
         }
     }
     printf("%d\n",ans);
 //    while(1);
     return ;
 }