容斥原理学习 M - Calculation 2 HDU - 3501

最新推荐文章于 2017-08-17 11:16:46 发布

小哧

最新推荐文章于 2017-08-17 11:16:46 发布

阅读量464

点赞数 1

分类专栏： acm水题文章标签：容斥原理数论

本文链接：https://blog.csdn.net/JayWRZZ/article/details/76855165

版权

acm水题专栏收录该内容

19 篇文章 0 订阅

订阅专栏

本博客介绍了使用质因数分解和容斥原理解决数论问题的方法，重点讨论了在解决HDU 3501题目时的思路，包括质因数对数的贡献计算以及由欧拉函数推导的解决方案。通过分析，阐述了如何避免计算中的错误并理解数论概念在实际问题中的应用。

摘要由CSDN通过智能技术生成

https://vjudge.net/contest/177165#rank
欢迎来刷我抄的专题.

因为数不大,很快就想到了质因数分解,然后容斥原理做一遍,但是取余的地方没细心.爆long long 却一直没发现.卡了很久.

思路大抵就是 : 质因数分解, 算每一个质因数的贡献,但我们平常所减的是求个数的,这里求完了个数后转变成和其实也蛮简单的, 我们这样考虑 1 - (n-1)中有a个p质因数的倍数, 那么这个质因数的贡献就是p*sum(1+2+..+a), 然后容斥消去重复的就行了.

#include<cstdio>
#include<iostream>
#include<bitset>
#include<cstring>
#define ll long long
const int maxn = 5e5;
using namespace std;
const int mod = 1000000007;
bitset<maxn> judge;
int prime[maxn];

int tot = 0;
ll l,r,n;
ll co[maxn];
void init(){
    judge.reset();
    tot = 0;
    for(int i = 2; i < maxn ; i++){
        if(!judge[i]){
            prime[tot++] = i;
        }
        for(int j = 0; j < tot ;j++){
            if(i*prime[j] >= maxn)break;
            judge[i*prime[j]] = 1;
            if(i%prime[j] == 0){
                break;
            }
        }
    }
}

ll sum(ll n){
    if(n%2==0){
        return ((n/2)%mod*((n+1)%mod))%mod;
    }else{
        return ((n+1)/2%mod*(n%mod))%mod;
    }
}

int cnt = 0;
void slove(ll n){
    cnt = 0;
    for(int i = 0; i < tot && prime[i] <= n;i ++){
        if(n%prime[i] == 0){
            co[cnt++] = prime[i];
            while(n%prime[i] == 0){
                n /= prime[i];
            }
        }
    }
    if(n>1){
        co[cnt++] = n;
    }
}

ll cal(ll a){
    ll ans = 0;
    ll temp = 1;
    int cc = 0;
    for(int i = 1; i < (1<<cnt); i++){
        temp = 1;
        cc = 0;
        for(int j = 0; j <cnt; j++){
            if(i & (1<<j)){
                temp *= co[j];
                cc ++ ;
            }
        }
        if(cc&1)ans += ((sum(a/temp)%mod)*(temp%mod))%mod;
        else ans -= ((sum(a/temp)%mod)*(temp%mod))%mod;
        ans = (ans%mod + mod)%mod;
    }
    return ans;
}

int main(){
    init();
    //FILE *fp;
    //fp = fopen("in.txt", "w+");
    //for(int r = 1; r < mod; r++){
    while(scanf("%lld",&r) && r){
        n = r;
        slove(n);
        printf("%lld\n",cal(r-1)%mod);
    }
}