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;
}