ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

CSUST 简单数学题 题解(质因子分解+并查集)

2021-07-31 20:34:48  阅读:173  来源: 互联网

标签:const int 题解 ll 查集 long 1e6 maxn CSUST


题目链接

题目大意

注意数据范围,因为区间范围最多\(1e6\),所以只要考虑\(1e6\)以内的质因子即可

然后用并查集维护一下集合即可

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
#define fi first
#define se second
#define debug printf("aaaaaaaaaaa\n");
const int maxn=1e6+5,inf=0x3f3f3f3f,mod=1e9+7;
const ll INF=0x3f3f3f3f3f3f3f3f;
const double eps=1e-7;
ll a,b,p;
int fa[maxn];
int prime[maxn],cnt;
bool isprime[maxn];
int tot;
void getprime(int n){
    for(ll i=2;i<=n;i++){//开ll因为后面要计算i*prime[j]
        if(!isprime[i]){
            prime[++cnt]=i;
        }
        for(int j=1;j<=cnt&&i*prime[j]<=n;j++){
            isprime[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
int findd(int x){
    return x==fa[x]?x:fa[x]=findd(fa[x]);
}
int main(){
    getprime(1000000);
    int _; cin>>_;
    while(_--){
        cin>>a>>b>>p;
        for(ll i=a;i<=b;i++){
            fa[i-a+1]=i-a+1;
        }
        for(int i=1;i<=cnt;i++){
            if(prime[i]<p) continue;
            ll beg=(a+prime[i]-1)/prime[i]*prime[i];
            for(ll j=beg;j<=b;j+=prime[i]){
                int x=findd(beg-a+1);
                int y=findd(j-a+1);
                if(x==y) continue;
                fa[y]=x;
            }
        }
        int ans=0;
        for(ll i=a;i<=b;i++){
            if(findd(i-a+1)==i-a+1){
                ans++;
            }
        }
        printf("Case #%d: %d\n",++tot,ans);
    }
    return 0;
}


标签:const,int,题解,ll,查集,long,1e6,maxn,CSUST
来源: https://www.cnblogs.com/hunxuewangzi/p/15085123.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有