记录编号	357750	评测结果	AAAAAAAAAA
题目名称	[BZOJ 2820] YY的GCD	最终得分	100
用户昵称	sxysxy	是否通过	通过
代码语言	C++	运行时间	12.820 s
提交时间	2016-12-12 15:10:38	内存使用	162.49 MiB

//wc2016教师授课讲义
#include <cstdio>
#include <cstring>
#include <cstdarg>
#include <queue>
#include <vector>
#include <list>
#include <algorithm>
using namespace std;
#define MAXN 10000001
bool vis[MAXN];
int primes[MAXN];
int miu[MAXN];
int tot;
typedef long long LL;
LL f[MAXN];  
#define MAXT 10002
struct que
{
    int n, m;
}qs[MAXT];
void pre(int n)
{
    miu[1] = 1;
    for(int i = 2; i <= n; i++)
    {
        if(!vis[i])
        {
            primes[tot++] = i;
            miu[i] = -1;
        }
        for(int j = 0; j < tot; j++)
        {
            int k = primes[j]*i;
            if(k > n)break;
            vis[k] = true;
            if(i%primes[j])miu[k] = -miu[i];
            else break;
        }
    }
    for(int i = 0; i < tot; i++)
        for(int j = 1; j*primes[i] <= n; j++)
            f[j*primes[i]] += miu[j];
    for(int i = 1; i <= n; i++)f[i] += f[i-1];
}
LL calc(int n, int m)
{
    LL ans = 0;
    if(n > m)swap(n, m);
    for(int i = 1, j; i <= n; i = j+1)
    {
        j = min(n/(n/i), m/(m/i));
        ans += (f[j]-f[i-1])*(n/i)*(m/i);
    }
    return ans;
}
int main()
{
    freopen("YYnoGCD.in", "r", stdin);
    freopen("YYnoGCD.out", "w", stdout);
    int T;
    scanf("%d", &T);
    int n = 0;
    for(int i = 0; i < T; i++)
    {
        scanf("%d %d", &qs[i].n, &qs[i].m);
        n = max(n, max(qs[i].n, qs[i].m));
    }
    pre(n);
    for(int i = 0; i < T; i++)
        printf("%lld\n", calc(qs[i].n, qs[i].m));
    return 0;
}