题意:给出n个模线性方程组 每一个给出一个 xi、yi、ki、ti
(|x[i]−x|+|y[i]−y|)modk[i]=t[i].
求解满足方程的(x,y)的对数,
n,m(1≤n≤10,1≤m≤10^9)
xi,yi,ki,ti(0≤ xi,yi ≤m,2 ≤ ki ≤ 5,0≤ti<ki)
思路:
其中k[i] = 2、3、4、5,其 lcm 为 60,我们可以根据 每一个 x[i] y[i],将整个的平面分为四部分,这样就可以去除绝对值的影响,总共分为n^2块 ,(放一个图),对每一个x,y都进行一次暴力组合,对每一个x,y都进行0~60的延伸即可(红色矩形)。
代码:
#include<bits/stdc++.h>
#define LL long long
using namespace std;
int x[15],y[15],k[15],t[15],xx[15],yy[15];
bool check(int xx,int yy,int n){
for(int i = 0;i<n;i++){
if((abs(x[i] - xx) + abs(y[i] - yy))%k[i] != t[i]) return false;
}
return true;
}
int cal(int l,int r){
r -= l+1;
if(r < 0){
return 0;
}
return r/60+1;
}
int main(){
int T,n,m;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&m);
for(int i = 0;i<n;i++){
scanf("%d%d%d%d",&x[i],&y[i],&k[i],&t[i]);
xx[i+1] = x[i];
yy[i+1] = y[i];
}
yy[n+1] = xx[n+1] = m+1;
sort(xx+1,xx+n+1);
sort(yy+1,yy+n+1);
LL ans = 0;
for(int i = 0;i<=n;i++){
if(xx[i] < xx[i + 1]){
for(int j = 0;j<=n;j++){
if(yy[j] < yy[j + 1]){
for(int dx = 0;dx <60;dx++){
for(int dy = 0;dy <60;dy ++){
if(check(xx[i] + dx,yy[j] + dy,n)){
ans += 1LL*cal(xx[i]+dx,xx[i+1])*cal(yy[j]+dy,yy[j+1]);
}
}
}
}
}
}
}
printf("%lld\n",ans);
}
}