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;
}