CSU1803 2016（同余）

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long LL;
/*
同余
题意:求满足a * b 是 2016的倍数的对数，  1<= a <=n, 1<= b <=m;
a[i]表示0-n之间的每个数对2016取余，余数为i的个数
(a*b)%2016==(a%2016)*(b%2016)
a=2016*k+x b=2016*k+y;
a*b%2016==0 => x*y%2016==0
*/
const int maxn=2016+5;
LL a[maxn],b[maxn];
LL n,m;
LL ans;

int main()
{
    while(~scanf("%lld %lld",&n,&m))
    {
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        //预处理余数为i的个数
        for(int i=0;i<2016;i++)
        {
            a[i]=n/2016;
            b[i]=m/2016;
        }
        for(int i=1;i<=n%2016;i++) a[i]++;
        for(int j=1;j<=m%2016;j++) b[j]++;
        ans=0;
         //x*y%2016==0
        for(int i=0;i<2016;i++)
        {
            for(int j=0;j<2016;j++)
            {
                if(i*j%2016==0)
                {
                    ans+=a[i]*b[j];
                }
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}