洛谷【P2257】 YY的GCD

出处：http://www.cnblogs.com/peng-ym/p/8652288.html   （  直接去出处那看就好了 ）

题目描述

• 神犇YY虐完数论后给傻×kAc出了一题
• 给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对
  kAc这种傻×必然不会了，于是向你来请教……
• 多组输入

#include<bits/stdc++.h>
#define LL long long
#define ULL unsigned long long
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define dep(i,j,k) for(int i=k;i>=j;i--)
#define INF 0x3f3f3f3f
#define mem(i,j) memset(i,j,sizeof(i))
#define make(i,j) make_pair(i,j)
#define pb push_back
using namespace std;
const int N=1e7+;
bool vis[N];
int pre[N],sum[N],mu[N],tot;
void init() {
    mu[]=;
    rep(i,,) {
        if(!vis[i]) { pre[++tot]=i; mu[i]=- ; }
        rep(j,,tot) {
             if(i*pre[j]>N-) break;
             vis[i*pre[j]]=;
             if(i%pre[j]==) break;
             mu[i*pre[j]]=-mu[i];
        }
    }
    rep(j,,tot) {
        for(int i=;i*pre[j]<=N-;i++) {
            sum[i*pre[j]]+=mu[i];
        }
    }
    rep(i,,N-) sum[i]+=sum[i-];
}
int main() {
    init();
    int t,n,m;
    scanf("%d",&t);
    while(t--) {
        scanf("%d %d",&n,&m); if(n>m) swap(n,m);
        LL ans=;
        for(int l=,r;l<=n;l=r+) {
            r=min(n/(n/l),m/(m/l));
            ans+=1LL*(n/l)*(m/l)*(sum[r]-sum[l-]);
        }
        printf("%lld\n",ans);
    }
    return ;
}