利用同余定理,首先计算出 k=n/2016 r=n%2016;
任何一个数都能表示成n=a*2016+i,m=b*2016+j; n*m=a*b*2016^2+a*2016*j+b*2016*i+i*j;
如果i*j是2016的倍数的话,那么n*m一定是2016的倍数, 这样在2016*2016的规模中遍历出符合i*j%2016=0的数,
按照n中有多少个这样的i,m中有多少这样的j,然后相乘,累加起来就得到的区间中所有的数。 这种需要两个for循环。
我们可以用一个就能实现 我们枚举n的余数i 然后2016/__gcd(i,2016) 就是m的所需的余数 然后用m/这个余数 就可以得到m中含有多少个所需余数的倍数 然后n中含有多少个i的余数*m中含有多少个所需的余数 求和就是答案。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
int get(int r, int n) {
return n / r;
}
int main()
{
int n, m;
while(~scanf("%d %d", &n, &m)) {
ll ans = 0;
int r = n % 2016, k = n / 2016;
for(int i = 1; i <= 2016; i++) {
ans += 1LL * (k + (i <