HDU 4596 J - Yet another end of the world（数论）

题意：

给你n个 $x_i$，$y_i$，$z_i$，让你判断是否存在一个数num， 满足形如下式的条件，满足的话则代表此飞船在这个虫洞的作用力下，那么这个num号飞船不能再同时再处于另一个冲动的引力下，否则会被潮汐力撕碎。。。现在问你是否有被撕碎的危险。

$$y <= num(mod x_i) <= z$$

思路：

遍历两两组合的全部。对于任意的两个有如下的特点。

$$num \mod x_1 = r_1 \tag1$$
$$num \mod x_2 = r_2 \tag2$$
即：
$$num = k_1*x_1 + r_1 = k_2*x_2 + r_2$$
经过化简：
$$k_1*x_1 - k_2*x_2 = r_2 - r_1 \tag3$$
(3)式是典型的二元一次不定方程，有解的充要条件是 $gcd(k_1,k_2) | (r_2-r_1)$ 即：( $r_2$ - $r_1$ ) % gcd( $k_1$, $k_2$) = 0

对于这两个区间来说，如果相交，那么必可选出相等的点使得 $r_2$ - $r_1$=0，从而有解。但不相交的情况也有可能有解，我们假设在左边的那个区间为y1，z1，右边的为y2，z2，这时 $r_2$ - $r_1$ 的取值范围就在 $[y_2-z_1,z_2-y_1]$ 间了，我们这时只需判断这个区间里是否有模g = gcd($k_1$,$k_2$)等于0的数就好了。

判断的方法：如果 $y_2-z_1$ % g 不为零并且 $z_2-y_1$ 并没有超过g的循环节。具体代码如下：

代码：

#include <iostream>
#include <cstdio>
#include <string.h>

typedef long long int lli;

using namespace std;
lli flag = 0;
lli gcd(lli a,lli b){
    return b == 0 ? a:gcd(b,a%b);
}
struct node{
    lli x,y,z;
}a[1050];
bool f(int ii,int jj){
    if( (a[ii].z >= a[jj].y && a[ii].y < a[jj].z)
        || (a[jj].z >= a[ii].y && a[jj].y < a[ii].z)){
        return true;
    }
    if(a[ii].y > a[jj].z){
        swap(ii,jj);
    }
    lli from = a[jj].y-a[ii].z;
    lli to = a[jj].z-a[ii].y;
    lli gc = gcd(a[ii].x,a[jj].x);
    if(from % gc != 0 && to-from <= gc-from%gc-1){
        return false;
    }
    return true;;
}
int main(){
    int n;
    while(~scanf("%d",&n)){
        flag = 0;
        for(int i = 1;i <= n;i++){
            scanf("%lld%lld%lld",&a[i].x,&a[i].y,&a[i].z);
        }
        for(int i = 1; i< n;i++){
            for(int j = i+1;j <= n;j++){
                if(f(i,j)){
                    flag = 1;
                    break;
                }
            }
            if(flag == 1) break;
        }
        if(flag == 0){
            puts("Can Take off");
        }
        else
            puts("Cannot Take off");
    }
}