BZOJ2820：YY的GCD——题解

http://www.lydsy.com/JudgeOnline/problem.php?id=2820

Description

神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种

傻×必然不会了，于是向你来请教……多组输入

Input

第一行一个整数T 表述数据组数接下来T行，每行两个正整数，表示N, M

Output

T行，每行一个整数表示第i组数据的结果

Sample Input

2
10 10
100 100

Sample Output

30
2791

————————————————————————

看hzw的博客吧，他讲的蛮清楚的……

我主要是没有可以写数学公式的东西……

http://hzwer.com/6142.html

#include<cstdio>
#include<queue>
#include<cctype>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=;
int miu[N],su[N],sum[N];
ll f[N];
bool he[N];
void Euler(int n){
    int tot=;
    miu[]=;
    for(int i=;i<=n;i++){
    if(!he[i]){
        su[++tot]=i;
        miu[i]=-;
    }
    for(int j=;j<=tot;j++){
        if(i*su[j]>=n)break;
        he[i*su[j]]=;
        if(i%su[j]==){
        miu[i*su[j]]=;break;
        }
        else miu[i*su[j]]=-miu[i];
    }
    }
    for(int i=;i<=tot;i++){
    int p=su[i];
    for(int j=;j*p<=n;j++)f[j*p]+=miu[j];
    }
    for(int i=;i<=n;i++)f[i]+=f[i-];
    return;
}
int main(){
    Euler();
    int t;
    scanf("%d",&t);
    while(t--){
    int a,b;ll ans=;
    scanf("%d%d",&a,&b);
    if(a>b)swap(a,b);
    for(int i=,j;i<=a;i=j+){
        j=min(a/(a/i),b/(b/i));
        ans+=(f[j]-f[i-])*(a/i)*(b/i);
    }
    printf("%lld\n",ans);
    }
    return ;
}