专栏首页数据结构与算法BZOJ2820: YY的GCD(反演)

BZOJ2820: YY的GCD(反演)

题解

题意

题目链接

Sol

反演套路题。。

不多说了,就是先枚举一个质数,再枚举一个约数然后反演一下。

最后可以化成这样子

\[\sum_{i = 1}^n \frac{n}{k} \frac{n}{k} \sum_{p \in P, p | k} \mu(\frac{K}{p})\]

然后后面的那一坨可以暴力预处理。。复杂度不清楚,但是显然严格小于调和级数,所以也没啥大问题。

/*

*/
#include<bits/stdc++.h>
#define LL long long 
//#define int long long 
const int MAXN = 1e7 + 10, INF = 1e9 + 7;
using namespace std;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int T, N, M, prime[MAXN], mu[MAXN], tot, vis[MAXN];
LL g[MAXN];
void Get(int N) {
    vis[1] = 1; mu[1] = 1;
    for(int i = 2; i <= N; i++) {
        if(!vis[i]) prime[++tot] = i, mu[i] = -1;
        for(int j = 1; j <= tot && i * prime[j] <= N; j++) {
            vis[i * prime[j]] = 1;
            if(i % prime[j]) mu[i * prime[j]] = -mu[i];
            else {mu[i * prime[j]] = 0; break;}
        }
    }
    for(int i = 1; i <= tot; i++) 
        for(int j = 1; prime[i] * j <= N; j++) g[prime[i] * j] += mu[j];
    for(int i = 1; i <= N; i++) g[i] += g[i - 1];
}
int calc(int K) {
    return g[K];
}
void solve() {
    N = read(); M = read();
    if(N > M) swap(N, M);
    LL ans = 0;
    for(int k = 1, j; k <= N; k = j + 1) {
        j = min(N / (N / k), M / (M / k));
        ans += 1ll * (N / k) * (M / k) * (g[j] - g[k - 1]);
    }
    cout << ans << '\n';
}
signed main() {
    Get(1e7);
    for(int T = read(); T; T--, solve());
    return 0;
}
/*
4
10 10
120 100
123 1234
10000000 10000000
*/

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • 线性筛莫比乌斯函数

    1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<cmath...

    attack
  • BZOJ1101: [POI2007]Zap(莫比乌斯反演)

    Description   FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a ,y<=b,并且gc...

    attack
  • BZOJ1053: [HAOI2007]反素数ant(爆搜)

      对于任何正整数x,其约数的个数记作g(x)。例如g(1)=1、g(6)=4。如果某个正整数x满足:g(x)>g(i) 0<i<x

    attack
  • 线性筛莫比乌斯函数

    1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<cmath...

    attack
  • BZOJ1101: [POI2007]Zap(莫比乌斯反演)

    Description   FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a ,y<=b,并且gc...

    attack
  • 【PAT甲级】1002 A+B for Polynomials (25分)

    This time, you are supposed to find A+B where A and B are two polynomials.

    韩旭051
  • 栈缓冲区溢出

    https://baike.baidu.com/item/%E8%8E%AB%E9%87%8C%E6%96%AF%E8%A0%95%E8%99%AB/90359...

    字节脉搏实验室
  • 洛谷2017 5月月赛R1

    我只想说面对这种难度的题目就是冲着20%的数据暴力。。。 分数:40+20+36.1+38+0+19 T1 签到题 III 题目背景 pj组选手zzq近日学会了...

    attack
  • 面向对象的演进过程

    我们知道 程序 = 数据结构 + 算法,其中数据结构包括数组、栈、队列、链表、树以及图等,而算法是包含顺序、循环、分支三种逻辑结构的代码,为了使算法能够到处复用...

    木可大大
  • 洛谷P2522 [HAOI2011]Problem b(莫比乌斯反演)

    题目描述 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数。 ...

    attack

扫码关注云+社区

领取腾讯云代金券