CSU 1803 2016(数论)

题解：

可以求出在 [1,N] 中，模为 [0,2016] 的数的个数; 以及在 [1,M] 中，模为 [0,2016] 的数的个数。复杂度 O(2016)
因为 x∗y%2016=x%2016∗y%2016    ，那么，接下来就只需要统计求和就OK了。

Description

 给出正整数 n 和 m，统计满足以下条件的正整数对 (a,b) 的数量：

1. 1≤a≤n,1≤b≤m;

2. a×b 是 2016 的倍数。

Input

输入包含不超过 30 组数据。

每组数据包含两个整数 n,m (1≤n,m≤10 ⁹).

Output

对于每组数据，输出一个整数表示满足条件的数量。

Sample Input

32 63
2016 2016
1000000000 1000000000

Sample Output

1
30576
7523146895502644

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<stack>
#include<string>
#include<queue>
#include<algorithm>
#include<vector>
#define LL long long
#define INF 99999999999999999;
#define maxn 200005
using namespace std;
int main()
{
	
	LL n, m;
	while (scanf("%lld%lld", &n, &m) != EOF)
	{
		LL a[2020] = { 0 }, b[2020] = {0};
		LL aa = n / 2016, bb = m / 2016;
		n = n % 2016, m = m % 2016;
		for (int i = 1; i <= n; i++)
			a[i]++;
		for (int i = 1; i <= m; i++)
			b[i]++;
		for (int i = 1; i <= 2016; i++)
		{
			a[i] += aa;
			b[i] += bb;
		}
		LL sum = 0;
		for (int i = 1; i <= 2016; i++)
			for (int j = 1; j <= 2016; j++)
			{
				if (i*j % 2016 == 0)
					sum += a[i] * b[j];
			}
		printf("%lld\n", sum);
	}
}