2020-11-15

最新推荐文章于 2022-11-20 15:26:36 发布

luo_wh

最新推荐文章于 2022-11-20 15:26:36 发布

阅读量238

点赞数

分类专栏：题解文章标签：数学

本文链接：https://blog.csdn.net/qq_44599840/article/details/109710346

版权

题解专栏收录该内容

2 篇文章 0 订阅

订阅专栏

ICPC江西省赛A题个人题解

暴力的做法

void solve(){
    int n;
    scanf("%d",&n);
    long long res = 0;
    for(int i = 1;i <= n;i++){
        for(int j = 1;j <= i;j++){
            if(gcd(i,j) == 1){
                //_val[j]为j的F值
                res += _val[j];
            }
        }
    }
    printf("%lld\n",res);
}

也可以这样理解：

void solve(){
    int n;
    scanf("%d",&n);
    long long res = 0;
    for(int i = 2;i <= n;i++){
        for(int j = 1;j <= i;j++){
            if(gcd(i,j) != 1){
                //_val[j]为j的F值
                sum[i] -= _val[j];
            }
        }
    }
    for(int i = 1;i <= n;i++)
        res += sum[i];
    printf("%lld\n",res);
}

其中 _val[MAX_N]来存正整数的 f值，long long sum[MAX_N] 来存_val 的前缀和。

int get_val(int num){
    int res = 0;
    while (num){
        res += num % 10;
        num /= 10;
    }
    return res;
}

_val[0] = 0;
    for(int i = 1;i <= n;i++){
        _val[i] = get_val(i);
        sum[i] = sum[i-1] + _val[i];
    }

思路

对于每一个质数 i 在这里插入图片描述
我们在excel中打表，棕色的是gcd为1的组合，蓝色的为不互质的组合；
我们可以看到，如果n能被2整除，所有小于n 的2的倍数都不与n互质。
我们可以发现，对每一个质数 x ，都与x的倍数不互质(这不废话吗。。。)
但对每一个非质数 x，都与他的质因数的倍数不互质。
例如
4：4与4，6，8都不互质；
6：6 与 8，9，10，12都不互质。

即我们可以这样解：
对一个质数p，所有p的倍数x 都要减掉 _val[p * k] (k为小于 x / p的所有正整数)

last = 0;
for(int j = 1;j <= n / p;j++){
        last += _val[p * j];
        sum[p * j] -= last;
}

但是又发现了一个问题，比如6，12，30这样，能被多个质数整除的数。多减了很多次。所以我们就要加回来。加的次数过多，又要减回去。
例如被2和3两个质数整除的数，都能被6整除，我们可以用6的倍数加回来。
又如能被2，3，5三个质数整除的数，都能被30整除，我们可以用30的倍数减回去。
在草稿纸上算发现，能被奇数个偶数整除的数要加,否则就减。
这个能被多个质数整除的最小的数就是 $p 1 * p 2 * p 3 . . * p k$
$1 - C_k^0 + C_k^1 - C_k^2 + C_k^3 ... - (-1)^kC_k^k = 1 + (1 - 1)^k = 1$

		last = 0;
		if(k % 2){
            for(int j = 1;j <= n / i;j++){
                last += _val[i * j];
                sum[i * j] -= last;
            }
        }else{
            for(int j = 1;j <= n / i;j++){
                last += _val[i * j];
                sum[i * j] += last;
            }
        }

如何找到这个 $p 1 * p 2 * p 3 . . . * p k$

质因数分解一下就好。如果i的质因数的次数都小于等于1，那他就是天选之子。否则就会重复计算。

int divide(int num){
    int t = num;
    int res = 0;
    for(int i = 1;i <= num / i;i++){
        if(t % prim[i] == 0){
            res++;
            t /= prim[i];
            if(t % prim[i] == 0)
                return 0;
        }
    }
    if(t > 1)
        res++;
    return res;
}

	for(int i = 2;i <= n;i++){
        ll last = 0;
        int k = divide(i);
        if(!k)
            continue;
        if(k % 2){
            for(int j = 1;j <= n / i;j++){
                last += _val[i * j];
                sum[i * j] -= last;
            }
        }else{
            for(int j = 1;j <= n / i;j++){
                last += _val[i * j];
                sum[i * j] += last;
            }
        }
    }

AC代码

#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <vector>

using namespace std;

const int MAX_N = 1e5 + 5;
typedef long long ll;

ll _val[MAX_N],sum[MAX_N];
int prim[MAX_N],tot,v[MAX_N];
bool vis[MAX_N];

void get_prim(){
    int N = MAX_N - 5;
    tot = 0;
    memset(v,0,sizeof v);
    for(int i = 2;i <= N;i++){
        if(!v[i]){
            v[i] = i;
            prim[++tot] = i;
        }
        for(int j = 1;j <= tot;j++){
            if(prim[j] > v[i] || prim[j] > N/i)
                break;
            v[i * prim[j]] = prim[j];
        }
    }
}

int divide(int num){
    int t = num;
    int res = 0;
    for(int i = 1;i <= num / i;i++){
        if(t % prim[i] == 0){
            res++;
            t /= prim[i];
            if(t % prim[i] == 0)
                return 0;
        }
    }
    if(t > 1)
        res++;
    return res;
}


int get_val(int num){
    int res = 0;
    while (num){
        res += num % 10;
        num /= 10;
    }
    return res;
}

void solve(){
    int n;
    scanf("%d",&n);
    memset(vis,0,sizeof vis);
    _val[0] = 0;
    for(int i = 1;i <= n;i++){
        _val[i] = get_val(i);
        sum[i] = sum[i-1] + _val[i];
    }
    ll res = 0;
    for(int i = 2;i <= n;i++){
        ll last = 0;
        int k = divide(i);
        if(!k)
            continue;
        if(k % 2){
            for(int j = 1;j <= n / i;j++){
                last += _val[i * j];
                sum[i * j] -= last;
            }
        }else{
            for(int j = 1;j <= n / i;j++){
                last += _val[i * j];
                sum[i * j] += last;
            }
        }
    }
    for(int i = 1;i <= n;i++){
        res += (ll)sum[i];
    }
    printf("%lld\n",res);
}

int main(){
    get_prim();
    solve();
    return 0;
}