bzoj1407: [Noi2002]Savage

这个题，把式子列出来，就不难发现是个拓展欧几里得求线性解，然后分类讨论一下解得情况就可以了。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
int gcd(int a,int b)
{
if(a<0) a=-a;
if(b<0) b=-b;
if(b==0)
{
return a;
}
return gcd(b,a%b);
}
int exgcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;y=0;
return a;
}
int tmp=exgcd(b,a%b,x,y);
int temp=y;y=x-a/b*y;x=temp;
return tmp;
}
int c[20],p[20],cnt,l[20],n,m,pp[20][20],cc[20][20];
bool solve(int a,int b)
{
int d,x,y;
d=exgcd(pp[a][b],m,x,y);
if(cc[a][b]%d!=0)
{
return true;
}
int t1=cc[a][b]/d;
int t2=m/d;
int ans=t1*x;
ans=(ans%t2+t2)%t2;
if(ans<=min(l[a],l[b])) return false;
return true;
}
bool check()
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<i;j++)
{
if(!solve(i,j))
{
return false;
}
}
}
return true;
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d%d%d",&c[i],&p[i],&l[i]);
c[i]--;
cnt=max(cnt,c[i]);
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<i;j++)
{
pp[i][j]=p[i]-p[j];
cc[i][j]=c[j]-c[i];
if(pp[i][j]<0)
{
pp[i][j]*=-1,cc[i][j]*=-1;
}
}
}
for(m=cnt+1;;m++)
{
if(check())
{
printf("%d\n",m);
return 0;
}
}
return 0;
}