2016 (同余定理）

题目：

给出正整数 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

同余定理介绍：

同余的定义是这样的：  两个整数，a,b,如果他们同时除以一个自然数m，所得的余数相同，则称a,b对于模m同余。。记作a≡b（mod.m）。读作：a同余于b模m。  同余的性质也比较多，主要有以下一些：

 1.对于同一个除数，两个数之和（或差）与它们的余数之和（或差）同余。 

2.对于同一个除数，两个数的乘积与它们余数的乘积同余。 

3.对于同一个除数，如果有两个整数同余，那么它们的差就一  定能被这个除数整除。 

4.对于同一个除数，如果两个数同余，那么他们的乘方仍然同余。

我们这里用的是第二条：对于同一个除数，两个数的乘积与它们余数的乘积同余。

参考代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

int main()
{
	int k,t;
	long long ans;
	long long m,n;
	long long a[2016],b[2016];
	
	while(scanf("%lld %lld",&m,&n) == 2)
	{
		memset(a,0,sizeof(a));
		memset(b,0,sizeof(b));
		ans = 0;
		
		k = m/2016;
		t = m%2016;
		for(int i = 1;i <= t;i ++)
		a[i] = k+1;
		for(int i = t+1;i <= 2016;i ++)
		a[i] = k;
		
		k = n/2016;
		t = n%2016;
		for(int i = 1;i <= t;i ++)
		b[i] = k+1;
		for(int i = t+1;i <= 2016;i ++)
		b[i] = k;
		
		for(int i = 1;i <= 2016;i ++)
		{
			for(int j = 1;j <= 2016;j ++)
			{
				if(i*j % 2016 == 0)
				ans += a[i]*b[j];
			}
		}
		printf("%lld\n",ans);
	}
	return 0;
}