zoukankan      html  css  js  c++  java
  • P1621 集合

    传送门

    解题思路

    因为是在[a,b]范围内将两个拥有大于等于P的公共质因数的整数进行合并,关于质数的操作我们很容易会联想到素数筛,实际上就是这样,我们先通过欧拉筛,然后我们找到第一个大于等于a的质数,然后将[a,b]范围内的该质数的倍数进行合并操作,最后从a到b数一下有多少个分类即可

    Code

    #include<bits/stdc++.h>
    using namespace std;
    #define endl "
    "
    
    const int N = 100005;
    
    int primes[N];
    bool vis[N];
    int b,a,p;
    int fa[N];
    
    
    void init() {
    	for(int i = 1;i < N; ++i) fa[i] = i;
    }
    
    void get_prime(int n) {
    	vis[0] = vis[1] = true;
    	for(int i = 2;i <= n; ++i) {
    		if(!vis[i]) 
    			primes[++primes[0]] = i;
    		for(int j = 1;j <= primes[0] && primes[j] * i <= n; ++j) {
    			vis[primes[j] * i] = true;
    			if(i % primes[j] == 0) break;
    		}
    	}
    }
    
    
    int find(int x) {
    	int t = x;
    	while(t != fa[t]) {
    		t = fa[t];
    	}
    	while(x != fa[x]) {
    		int temp = fa[x];
    		fa[x] = t;
    		x = temp;
    	}
    	return x;
    }
    
    void merge(int l,int r) {
    	l = find(l);
    	r = find(r);
    	if(l != r) {
    		fa[r] = l;
    	}
    }
    
    int main()
    {
    	scanf("%d%d%d",&a,&b,&p);
    	init();
    	get_prime(b);
    	int i = lower_bound(primes+1,primes+1+primes[0],p) - primes;
    //	while(primes[i] < a) i++;
    	for(;i < primes[0] && primes[i] <= b; ++i) {
    		int k = primes[i];
    		for(int j = 1;j <= b && j * k <= b; ++j) {
    				merge(k,j * k);
    		}
    	}
    	set<int> ans;
    	for(i = a; i <= b; ++i) {
    //		printf("i = %d  find() = %d
    ",i,find(i));
    		ans.insert(find(i));
    	}
    	printf("%d
    ",ans.size());
    	return 0;
    }
    /*
    10 20 3
    */ 
    
  • 相关阅读:
    js正则匹配
    包含HTML的字符串去掉HTML标签
    smart-table 服务端请求真分
    禁用H5 表单验证novalidate
    webpack
    linux 进程查看及杀死进程
    配置ca服务器和http,mail加密
    mysql权限
    mysql查询
    mysql储存引擎
  • 原文地址:https://www.cnblogs.com/Mangata/p/14729078.html
Copyright © 2011-2022 走看看