HDU 5212 Code【莫比乌斯反演】

原创 2016年05月30日 23:05:02

题目链接:

http://acm.hdu.edu.cn/showproblem.php?pid=5212

题意:

给定序列,1i,jn,求gcd(a[i],a[j])(gcd(a[i],a[j])1)之和。

分析:

同样我们设
f(d):满足gcd(x,y)=dx,y均在给定范围内的(x,y)的对数。
F(d):满足d|gcd(x,y)x,y均在给定范围内的(x,y)的对数。
反演后我们得到

f(x)=x|dμ(d/x)F(d)

由于序列给定,这里的F(d)我们可以通过枚举d,来找d的倍数的个数,那么F(d)=cnt[d]cnt[d],枚举最大公约数求出f(d),那么答案即为f(d)d(d1)的和。时间复杂度O(nlogn)

代码:

/*
-- Hdu 5212
-- mobius
-- Create by jiangyuzhu
-- 2016/5/30
*/
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <cmath>
#include <stack>
using namespace std;
typedef long long ll;
#define sa(n) scanf("%d", &(n))
#define sal(n) scanf("%I64d", &(n))
#define pl(x) cout << #x << " " << x << endl
#define mdzz cout<<"mdzz"<<endl;
const int maxn = 1e4+ 5 , mod = 1e4 + 7;
int tot = 0;
int miu[maxn], prime[maxn], a[maxn];
int cnt[maxn], F[maxn];
bool flag[maxn];
void mobius()
{
    miu[1] = 1;
    tot = 0;
    for(int i = 2; i < maxn; i++){
        if(!flag[i]){
            prime[tot++] = i;
            miu[i] = -1;
            cnt[i] = 1;
        }
        for(int j = 0; j < tot && i * prime[j] < maxn; j++){
            flag[i * prime[j]] = true;
            cnt[i * prime[j]] = cnt[i] + 1;
            if(i % prime[j]){
                miu[i * prime[j]] = -miu[i];
            }
            else{
                miu[i * prime[j]] = 0;
                break;
            }
        }
    }
}
int main (void)
{
    mobius();
    int n;
    while(~sa(n)){
        int maxa = 0;
        memset(cnt, 0, sizeof(cnt));
        memset(F, 0, sizeof(F));
        for(int i = 0; i < n; i++) {
            sa(a[i]);
            cnt[a[i]]++;
            maxa = max(maxa, a[i]);
        }
        for(int i = 1; i <= maxa; i++){
            for(int j = i; j <= maxa; j += i){
                F[i] += cnt[j];
            }
        }
        ll ans = 0;
        ll tmp = 0;
        for(int i = 1; i <= maxa; i++){
            tmp = 0;
            for(int j = i; j <= maxa; j += i){
                 tmp += miu[j/ i]  *  F[j] * 1ll * F[j] % mod;
            }
            ans =( ans + tmp * 1ll * i % mod * (i  - 1)% mod) % mod;
        }
        printf("%I64d\n", ans);
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

hdu 5212 Code (容斥)

题意: int calc {      int res=0;      for(int i=1;i          for(int j=1;j          {              res...
  • My_ACM_Dream
  • My_ACM_Dream
  • 2015年04月26日 21:52
  • 538

HDU 5212

来自我的新博客HDU 5212Description:      ~~~~~~ 给定一个数列{aN}\{a_N\} ,然后求 Ans=∑i=1N∑j=1Ngcd(ai,aj)∗(gcd(ai,aj...
  • qq_21995319
  • qq_21995319
  • 2015年09月15日 21:39
  • 413

HDU5212 Code

题意:输出那个程序的结果 这题其实就是求任意一个数与每一个数的gcd×(gcd-1),然后全部加起来。。 把每个数值对应的个数记下来,然后从最大的数开始,找这个数的倍数,如果gcd是x那个构成这个...
  • eeeaaaaa
  • eeeaaaaa
  • 2015年04月27日 12:43
  • 238

hdu 5212(容斥原理)

Code Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submi...
  • u013268685
  • u013268685
  • 2015年04月27日 01:43
  • 940

hdu5212 Code 莫比乌斯

Problem Description WLD likes playing with codes.One day he is writing a function.Howerver,his co...
  • u013044116
  • u013044116
  • 2015年10月22日 21:53
  • 270

HDU 5212 Code(容斥 或 莫比乌斯反演)

Code Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) ...
  • LYHVOYAGE
  • LYHVOYAGE
  • 2015年04月27日 09:49
  • 2215

HDU 5212 Code(容斥)

题意: 输入n个数, 求这n个数的n^2个两两组合的最大公约数x,x(x-1)之和 。 换个思路, 我们如果求出了最大公约数为x的对数cnt, 那么答案就增加x(x-1)cnt^2 那么怎么求以x为最...
  • weizhuwyzc000
  • weizhuwyzc000
  • 2015年12月12日 15:09
  • 665

hdu 5212(容斥原理)

Code Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submi...
  • u013268685
  • u013268685
  • 2015年04月27日 01:43
  • 940

Code HDU 5212(容斥原理)

Code Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submissi...
  • u013167299
  • u013167299
  • 2015年04月27日 22:05
  • 464

hdu 5212(莫比乌斯反演)

Code Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Pro...
  • hexianhao
  • hexianhao
  • 2016年07月12日 23:05
  • 169
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 5212 Code【莫比乌斯反演】
举报原因:
原因补充:

(最多只允许输入30个字)