[HDU1695]GCD（莫比乌斯反演）

最新推荐文章于 2017-04-12 07:43:21 发布

Clove_unique

最新推荐文章于 2017-04-12 07:43:21 发布

阅读量668

点赞数

分类专栏：题解莫比乌斯反演文章标签： HDU 数论

本文链接：https://blog.csdn.net/clove_unique/article/details/51039266

版权

题解同时被 2 个专栏收录

1097 篇文章 4 订阅

订阅专栏

莫比乌斯反演

22 篇文章 0 订阅

订阅专栏

题目描述

传送门

题解

又是根据题目化出来一坨式子。
$\sum\limits_{i=1}^b\sum\limits_{j=1}^d[(i,j)=k]$
$=\sum\limits_{i=1}^b\sum\limits_{j=1}^d[({i\over k},{j\over k})=1]$
$=\sum\limits_{i=1}^b\sum\limits_{j=1}^d\sum\limits_{t|({i\over k},{j\over k})}\mu (t)$
$=\sum\limits_{i=1}^b\sum\limits_{j=1}^d\sum\limits_{t=1}^{min({i \over k},{j\over k})}[t|{i\over k}][t|{j\over k}] \mu (t)$
令i=mk，j=nk
$=\sum\limits_{t=1}^{min(m,n)}\sum\limits_{m=1}^{b\over k}[t|m]\sum\limits_{n=1}^{d\over k}[t|n] \mu (t)$
$=\sum\limits_{t=1}^{min(m,n)}\lfloor{{b\over k}\over t}\rfloor\lfloor{{d\over k}\over t}\rfloor \mu (t)$
这样看起来就美观一点了。
听说可以分块求？
反正我打的暴力。
去重的话可以把重复的那一段单独求，然后/2，得到了重复的数，用原答案减去这个数，就为最后的答案。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
#define LL long long

int T,a,b,c,d,k;
int n,m,mu[100005],p[100005],prime[100005];
LL ans1,ans2; 

inline void get_mu(int n){
    mu[1]=1;
    for (int i=2;i<=n;++i){
        if (!p[i]){
            prime[++prime[0]]=i;
            mu[i]=-1;
        }
        for (int j=1;j<=prime[0]&&i*prime[j]<=n;++j){
            p[i*prime[j]]=1;
            if (i%prime[j]==0){
                mu[i*prime[j]]=0;
                break;
            }
            else mu[i*prime[j]]=-mu[i];
        }
    }
} 

int main(){    
    get_mu(100000);
    scanf("%d",&T);
    for (int t=1;t<=T;++t){
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        if (!k||!b||!d){
            printf("Case %d: 0\n",t);
            continue;
        }
        if (b>d) swap(b,d);
        n=(int)b/k; m=(int)d/k; ans1=ans2=0;
        for (int i=1;i<=n;++i)
          ans1+=(LL)(n/i)*(m/i)*mu[i];
        for (int i=1;i<=n;++i)
          ans2+=(LL)(n/i)*(n/i)*mu[i];
        printf("Case %d: %I64d\n",t,ans1-ans2/2);
    }
}