[2017浙工大之江学院决赛 L] qwb与整数对（离线，筛）

题目链接：http://115.231.222.240:8081/JudgeOnline/problem.php?cid=1005&pid=11

这题不会，看了柠檬巨的题解才知道可以筛出来。

枚举a、b，然后求下m最小能是多少，可以满足原式为整数。可以用(a*b)-(a*a+b*b)%(a*b)求。然后在m的基础上算上a*b的倍数再更新到对应case上即可。

#include <bits/stdc++.h>
using namespace std;

const int maxn = 6666;
const int maxm = 44444;
const int Q =  20000;
int n[maxn], m[maxn], ret[maxn];
vector<int> id[maxm];

int main() {
    // freopen("in", "r", stdin);
    for(int i = 1; ~scanf("%d%d",&n[i],&m[i]) && (n[i]||m[i]); i++) {
        id[m[i]+Q].push_back(i);
    }
    for(int a = 1; a <= 1000; a++) {
        for(int b = a+1; b <= 1000; b++) {
      int M = (a*b)-(a*a+b*b)%(a*b);
            for(int t = 0; M+t*(a*b) <= Q; t++) {
                for(int i = 0; i < id[M+t*(a*b)+Q].size(); i++) {
                    if(b < n[id[M+t*(a*b)+Q][i]]) ret[id[M+t*(a*b)+Q][i]]++;
                }
            }
            for(int t = 1; M-t*(a*b) >= -Q; t++) {
                for(int i = 0; i < id[M-t*(a*b)+Q].size(); i++) {
                    if(b < n[id[M-t*(a*b)+Q][i]]) ret[id[M-t*(a*b)+Q][i]]++;
                }
            }
        }
    }
    for(int i = 1; n[i]||m[i]; i++) {
        printf("Case %d: %d
", i, ret[i]);
    }
    return 0;
}