2017年第0届浙江工业大学之江学院程序设计竞赛决赛—L

题目链接:http://115.231.222.240:8081/JudgeOnline/problem.php?cid=1005&pid=11 
—————————————————————————————————————————— 
Problem L: qwb与整数对 
Time Limit: 1 Sec Memory Limit: 128 MB 
Submit: 171 Solved: 18 
[Submit][Status][Web Board] 
Description 
qwb又遇到了一道数学难题，你能帮助他吗？

给出两个整数n和m，请统计满足0＜a＜b＜n并且使得 (a2+b2+m)/(ab) 的结果是整数的整数对(a,b)的个数。

Input 
本题包含多组测试例 。当测试例数据是n=m=0时，表示输入结束。（测试例数量<6000）

每个测试例一行，是两个整数n和m。输入保证0≤n≤1000，-20000< m< 20000。

Output 
对每个测试例，输出测试例的编号（Case X:Y） 以及满足上述要求的整数对(a,b)的个数，输出格式如样例输出所示。

Sample Input 
10 1 
20 3 
30 4 
0 0 
Sample Output 
Case 1: 2 
Case 2: 4 
Case 3: 5 
——————————————————————————————————————————

在线的复杂度实在不科学

看dalao离线做了,我也做了下

离线之后n^2枚举到每个位置计算贡献的时候类似筛法只要找(i*j)的倍数就好了

然后记下这个时候的值询问,满足n的情况就ans++

感觉没有代码直观啊.

附本题代码 
——————————————————————————————————————————

#include <bits/stdc++.h>
typedef long long int LL;
using namespace std;

int n[6666],m[6666],ans[6666];
vector<int>vc[44444];

int main(){
    for(int i=1;~scanf("%d%d",&n[i],&m[i])&&(n[i]||m[i]);i++){
        vc[m[i]+20000].push_back(i);
    }
    for(int i=1;i<=1000;i++){
        for(int j=i+1;j<=1000;j++){
            int x=(i*j)-(i*i+j*j)%(i*j);
            for(int k=x;k<=20000;k+=i*j){
                int vcz=vc[k+20000].size();
                for(int u=0;u<vcz;u++){
                    if(n[vc[k+20000][u]]>j)
                        ans[vc[k+20000][u]]++;
                }
            }
            for(int k=x-i*j;k>=-20000;k-=i*j){
                int vcz=vc[k+20000].size();
                for(int u=0;u<vcz;u++){
                    if(n[vc[k+20000][u]]>j)
                        ans[vc[k+20000][u]]++;
                }
            }
        }
    }

    for(int i=1;(n[i]||m[i]);i++)  printf("Case %d: %d\n",i,ans[i]);

    return 0;
}
   
   
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
   
   
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35