hdu 6169 - 容斥

最新推荐文章于 2018-09-10 11:05:50 发布

a1214034447

最新推荐文章于 2018-09-10 11:05:50 发布

阅读量652

点赞数

分类专栏：容斥原理

本文链接：https://blog.csdn.net/a1214034447/article/details/77509400

版权

容斥原理专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题目链接：点击打开链接

题解思路:这题就是要求K要是质数，且要小于R，要么肯定是0.那么我们还的知道要是K是sqrt(n)以上的质数那么值肯定是它本身，那么我们只需要先把sqrt(10^11)内的质数求出来就行了。然后去容斥2-K里面的这些质数就可以了，可以用等差公式直接求出质数i的倍数的和，因为数据很大也不适合乌斯和其他，直接用递归容斥更好解决。

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<math.h>
#include<queue>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int mx = 1e6+10,mod = 1e9+7;
ll m,n,q;
int pri[mx];
bool vis[mx];
void init(){
    int top = 0;
    for(int i=2;i<mx;i++){
        if(!vis[i]) pri[top++] = i;
        for(int j=0;j<top&&pri[j]*i<mx;j++){
            vis[pri[j]*i] = 1;
            if(i%pri[j]==0) break;
        }
    }
}
ll dfs(ll val,int ts,ll sum,int pos){
    ll cnt = sum/val,t1,t2;
    if(cnt&1) t1 = cnt%mod,t2 = ((cnt+1)/2)%mod;
    else t1 = (cnt/2)%mod,t2 = (cnt+1)%mod;
    ll ans = ts*((t1*t2%mod)*val%mod);
    for(int i=pos;pri[i]!=q;i++){
        if(val*pri[i]>sum) break;
        ll v = val*pri[i];
        ans = (ans+dfs(v,-ts,sum,i+1)+mod)%mod;
    }
    return ans;
}
bool isprime(ll x){
    int sq = sqrt(1.0*x);
    for(int i=2;i<=sq;i++)
    if(x%i==0) return 0;
    return 1;
}
int main(){
    int t,cas = 1;
    scanf("%d",&t);
    init();
    while(t--){
        scanf("%lld%lld%lld",&m,&n,&q);
        ll sq = sqrt(1.0*n);
        if(!isprime(q)||q>n) printf("Case #%d: 0\n",cas++);
        else if(q>sq){
            if(m>q) printf("Case #%d: 0\n",cas++);
            else printf("Case #%d: %lld\n",cas++,q%mod);
        }else{
            ll cnt1 = dfs(q,1,n,0);
            ll cnt2 = dfs(q,1,m-1,0);
            printf("Case #%d: %lld\n",cas++,(cnt1-cnt2+mod)%mod);
        }
    }
    return 0;
}