bzoj 1407: [Noi2002]Savage【扩展欧几里得+中国剩余定理】

首先答案不会很大,所以枚举答案m,于是把问题转为了判定;

关于如何判定:

首先题目中虽然没说但是数据是按照初始洞穴编号排的序,所以并不用自己重新再排

假设当前答案为m,相遇时间为x,野人i和j,那么可以列出同余式;

\[x(p[i]-p[j])\equiv c[j]-c[i](mod\ m)
\]

\[x(p[i]-p[j])+ym=c[j]-c[i]
\]

于是可解exgcd。由于并不是互质的,所以最后的算天数需要m/d

#include<cstdio>
using namespace std;
const int N=20;
int n,mx,a[N],b[N],c[N],d,x,y;
void exgcd(int a,int b,int &d,int &x,int &y)
{
if(!b)
{
d=a,x=1,y=0;
return;
}
exgcd(b,a%b,d,y,x);
y-=a/b*x;//!!
}
bool ok(int m)
{
for(int i=1;i<n;i++)
for(int j=i+1;j<=n;j++)
{
int aa=((a[j]-a[i])%m+m)%m,bb=((b[i]-b[j])%m+m)%m;
exgcd(bb,m,d,x,y);
if(aa%d)
continue;
aa/=d;
int mm=m/d,k=(aa*x%mm+mm)%mm;
if(k<=c[i]&&k<=c[j])
return 0;
}
return 1;
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d%d%d",&a[i],&b[i],&c[i]);
mx=max(a[i],mx);
a[i]--;
}
for(int i=mx;;i++)
if(ok(i))
{
printf("%d\n",i);
break;
}
return 0;
}
上一篇:Hadoop开发第4期---分布式安装


下一篇:Hadoop简介与分布式安装