数学专项number_theory：UVa 718

很容易判断哪些电梯可以到达A，哪些电梯可以到到B。然后两两判断电梯是否可以在某一层相连，即建个图。这里就需要解一个具有两个方程的模方程组了。

这里，我分了两种情况。首先是当两者的X互质，即模互质时可以直接用中国剩余定理。而当两者不互质时，则采用枚举的方法。

最后再以可以到达A的电梯进行dfs就行了。只不过还是很慢，近1s，估计是想复杂了……

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long LL;
LL f,a,b;
int e;
LL x[110],y[110];
vector<int> start;
int end[110],vis[110];
vector<int> g[110];
LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}
void exgcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b) {d=a;x=1;y=0;}
    else
    {
        exgcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}
int china(int p1,int p2)
{
    LL M=x[p1]*x[p2];
    LL ans=0,d,q,w;
    w=M/x[p1];
    exgcd(x[p1],w,d,d,q);
    ans=(ans+w*q*(y[p1]%x[p1]))%M;
    w=M/x[p2];
    exgcd(x[p2],w,d,d,q);
    ans=(ans+w*q*(y[p2]%x[p2]))%M;
    while(ans<f)
    {
        if(ans>=y[p1]&&ans>=y[p2]) return 1;
        ans+=M;
    }
    return 0;
}
int judge(int p1,int p2)
{
    if(y[p1]<y[p2]) swap(p1,p2);
    LL g=x[p2]*gcd(x[p1],x[p2]);
    LL i=0;
    while(y[p1]+i*x[p1]<f&&i<g)
    {
        if((y[p1]-y[p2]+i*x[p1])%x[p2]==0) return 1;
        i++;
    }
    return 0;
}
int dfs(int u)
{
    vis[u]=1;
    if(end[u]) return 1;
    int ans=0;
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if(!vis[v]) ans|=dfs(v);
    }
    return ans;
}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        start.clear();
        memset(end,0,sizeof(end));
        cin>>f>>e>>a>>b;
        for(int i=0;i<e;i++) g[i].clear();
        for(int i=0;i<e;i++) cin>>x[i]>>y[i];
        for(int i=0;i<e;i++)
        {
            if(a>=y[i]&&(a-y[i])%x[i]==0) start.push_back(i);
            if(b>=y[i]&&(b-y[i])%x[i]==0) end[i]=1;
        }
        for(int i=0;i<e;i++)
        {
           for(int j=0;j<e;j++)
           {
               if(i==j) continue;
               if(gcd(x[i],x[j])==1)
               {
                   if(china(i,j))
                   {
                       g[i].push_back(j);
                       g[j].push_back(i);
                   }
               }
               else
               {
                   if(judge(i,j))
                   {
                       g[i].push_back(j);
                       g[j].push_back(i);
                   }
               }
           }
        }
        int flag=0;
        for(int i=0;i<start.size();i++)
        {
            memset(vis,0,sizeof(vis));
            if(dfs(start[i])) {flag=1;break;}
        }
        if(flag) cout<<"It is possible to move the furniture."<<endl;
        else cout<<"The furniture cannot be moved."<<endl;
    }
    return 0;
}