NOI 2002 荒岛野人

NOI 2002 荒岛野人

人生第一次做NOI的题祭！！！

大概是NOI最简单的一道题

克里特岛以野人群居而著称。岛上有排列成环行的M个山洞。这些山洞顺时针编号为1,2,…,M。岛上住着N个野人，一开始依次住在山洞C1,C2,…,CN中，以后每年，第i个野人会沿顺时针向前走Pi个洞住下来。

每个野人i有一个寿命值Li，即生存的年数。

下面四幅图描述了一个有6个山洞，住有三个野人的岛上前四年的情况。三个野人初始的洞穴编号依次为1，2，3；每年要走过的洞穴数依次为3，7，2；寿命值依次为4，3，1。

（图片就不考了）

 

我们发现这道题奇怪的给了答案的数据规模，这就说明我们可以枚举答案。

然后我们发现n的数据规模很小，n^2*m似乎可以过。

我们对于每一个答案，枚举任意两个人是否不会相遇。

对于每两个野人，我们可以根据他们的C和P列出同余方程：

Ci+Pi*x≡Cj+Pj*x(mod b)  其中x是年数，b是枚举的答案。

然后，将这个同余方程转化为不定方程，然后求解。

Ci+Pi*x=Cj+Pj*x+b*-y（-y是一个构造的新的未知数）

(Pi-Pj)*x+b*y=Cj-Ci

然后用拓展欧几里得求解（这是解析）求解

最后一定要注意！！！（因为在这里卡了几个小时）求完特解以后，要保证最终解为正。

通常情况，只需要(x*(Cj-Ci)/d(最小公约数）+b/d)%b/d就可以了。

但是注意Cj-Ci可能为负，而且c++的负数取模的机制是取决于被模数，被模数为正则余数为正，为负则反之。

而因为Cj-Ci可能为负，可能导致最终答案依然为负。

故此，我们将x*(Cj-Ci)/d先模b/d，保证其绝对值是小于b/d的，但是，b/d也可能为负（因为d可能为负），所以即使加上一个b/d也不能保证结果为正。

因此，我们把最后加上的b/d变成一个它的倍数且保证它为正的，可以替换为b或者b/d的绝对值。

代码如下：

 

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#define MAXN 10010
#define in(a) a=read()
#define REP(i,k,n)  for(int i=k;i<=n;i++)
using namespace std;
inline int read(){
    int x=0,f=1;
    char ch=getchar();
    for(;!isdigit(ch);ch=getchar())
        if(ch=='-')
            f=-1;
    for(;isdigit(ch);ch=getchar())
        x=x*10+ch-'0';
    return x*f;
}
int n,m;
int c[100010],p[100010],l[100010];
int t;
inline void exgcd(int a,int b,int &d,int &x,int &y){
    if(b==0)  x=1,y=0,d=a;
    else  exgcd(b,a%b,d,x,y),t=x,x=y,y=t-a/b*y;
}
inline bool check(int k){
    REP(i,1,n)
        REP(j,i+1,n){
            int a=p[i]-p[j],b=k,x,y,d;
            exgcd(a,b,d,x,y);
            if((c[j]-c[i])%d!=0)  continue;
            x=((x*(c[j]-c[i])/d)%(b/d)+(abs(b/d)))%(b/d);//保证x为正
            if(x<=min(l[i],l[j]))  return 0;
        }
    return 1;
}
int main(){
    in(n);
    REP(i,1,n)  in(c[i]),in(p[i]),in(l[i]),m=max(m,c[i]);
    REP(i,m,1000000) if(check(i)){ cout<<i<<endl;  return 0;}
    return 0;
}

 

posted @ 2019-03-25 15:31 Dijkstra·Liu 阅读( ...) 评论( ...) 编辑 收藏