能量采集 BZOJ - 2005

最新推荐文章于 2019-11-01 19:12:22 发布

GoneWithTheWind_yin

最新推荐文章于 2019-11-01 19:12:22 发布

阅读量349

点赞数

分类专栏： BZOJ 数学莫比乌斯容斥欧拉函数

本文链接：https://blog.csdn.net/GoneWithTheWind_yin/article/details/78251128

版权

数学同时被 3 个专栏收录

17 篇文章 0 订阅

订阅专栏

BZOJ

13 篇文章 0 订阅

订阅专栏

莫比乌斯

2 篇文章 0 订阅

订阅专栏

题目传送门

思路：这个题目可以用容斥来做
我们令f[x]为gcd(i, j) = x的个数
令g[x]为i, j因子为x的个数(n / x) * (m / x)
我们再把其中最大公约数不是x的减掉就是要求的答案了

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <list>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>

#define MAXN 100010
#define MAXE 5
#define INF 0x7ffffff
#define MOD 100003
#define LL long long
#define ULL unsigned long long
#define pi 3.14159

using namespace std;

LL f[MAXN];

int main() {
    std::ios::sync_with_stdio(false);
    LL n, m;
    cin >> n >> m;
    LL sum = 0;
    for (LL i = min(n, m); i >= 1; --i) {
        f[i] = (n / i) * (m / i);
        for (LL j = 2 * i; j <= min(n, m); j += i) {
            f[i] -= f[j];
        }
        sum += f[i] * (2 * i - 1);
    }
    cout << sum << endl;
    return 0;
}