HDU 6639 Faraway【数学】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6639
题目大意：给定$n,m$以及$n$组$x_i,y_i,k_i,t_i$，且条件为$(|x_i-x_e|+|y_i-y_e|)\%k_i=t_i$，计算满足条件的所有的$x_e,y_e$对数。

思路：用$x_i,y_i$将所有区域分割成$(n+1{)}^2$个区域，然后枚举每一个区域，因为模数$k_i$在$[2,5]$内，$lcm(2,3,4,5)=60$，那么对于每个区域中，满足条件的点是呈最大步长为60的循环的。那么就可以暴力枚举$60\ast 60$的小矩阵中的所有整数点，对符合条件的点计数，然后可以$O(1)$算出当前区域内所有满足条件的点数，所有区域的点数和就是答案。

时间复杂度：$O({60}^2\ast n^3)$

AC代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int n,m;
ll xi[15],yi[15],ki[15],ti[15],tx[15],ty[15];
inline int abs(int x){return x>0?x:-x;}
bool check(int x,int y){
    for(int i=0;i<n;i++){
        if((abs(xi[i]-x)+abs(yi[i]-y))%ki[i]!=ti[i]){
            return false;
        }
    }
    return true;
}
int cnt(int l,int r){//区间长度为[l,r]的有多少个数满足条件
    r-=l+1;
    if(r<0){
        return 0;
    }
    return r/60+1; 
}
int main()
{
    int t;
    cin>>t;
    while(t--){
        cin>>n>>m;
        int nn=1;
        tx[0]=0;//下边界
        ty[0]=0;
        tx[1]=m+1;//上边界
        ty[1]=m+1;
        for(int i=0;i<n;i++){
            cin>>xi[i]>>yi[i]>>ki[i]>>ti[i];
            tx[++nn]=xi[i];
            ty[nn]=yi[i];
        }
        sort(tx+1,tx+nn+1);
        sort(ty+1,ty+nn+1);
        ll ans=0;
        for(int i=0;i<nn;i++){
            if(tx[i]<tx[i+1]){//分割的矩阵区域不是一条线(合法)
                for(int j=0;j<nn;j++){
                    if(ty[j]<ty[j+1]){
                        for(int x=0;x<60;x++){//暴力从0~59枚举x,y
                            for(int y=0;y<60;y++){
                                if(check(tx[i]+x,ty[j]+y)){
                                    ans=ans+(ll)cnt(tx[i]+x,tx[i+1])*cnt(ty[j]+y,ty[j+1]);
                                }
                            }
                        }
                    }
                }
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}