csu1803——2016（数学）

Description

给出正整数 n 和 m，统计满足以下条件的正整数对 (a,b) 的数量：
1. 1≤a≤n,1≤b≤m;
2. a×b 是 2016 的倍数。
Input

输入包含不超过 30 组数据。
每组数据包含两个整数 n,m (1≤n,m≤109).
Output
对于每组数据，输出一个整数表示满足条件的数量。
Sample Input
32 63
2016 2016
1000000000 1000000000
Sample Output
1
30576
7523146895502644

硬算肯定不行。注意到1~2016和2017~4032乘其他任何区间的结果是一样的，因为是两倍的关系，所以1~4032乘任何区间的结果就是上述两个区间的两倍，所以只要统计1~2016的倍数就行，多出来的部分再加到前面，最后只要统计1~2016区间的就行

#include <iostream>
#include <cstring>
#include <string>
#include <vector>
#include <queue>
#include <cstdio>
#include <set>
#include <math.h>
#include <algorithm>
#include <queue>
#include <iomanip>
#include <ctime>
#define INF 0x3f3f3f3f
#define MAXN 100005
#define Mod 1000000007
using namespace std;
long long a[3000],b[3000];
int main()
{
    long long n,m;
    while(~scanf("%lld%lld",&n,&m))
    {
        int k=n/2016;
        int t=n%2016;
        for(int i=1;i<=t;++i)
            a[i]=k+1;
        for(int j=t+1;j<=2016;++j)
            a[j]=k;
        k=m/2016;
        t=m%2016;
        for(int i=1;i<=t;++i)
            b[i]=k+1;
        for(int j=t+1;j<=2016;++j)
            b[j]=k;
        long long ans=0;
        for(int i=1;i<=2016;++i)
        {
            for(int j=1;j<=2016;++j)
            {
                if((i*j)%2016==0)
                    ans=ans+a[i]*b[j];
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}