千鸽笼 题解

千鸽笼

时间限制: 1000ms

空间限制: 53687091kB

[题目描述]

你有 $(B-A+1)$ 只鸽子，它们的编号为从 $A$ 到 $B$ 的整数。一开始每只鸽子都被单独关在一个笼子里。

现在为了节省空间，你需要将一些鸽子赶到另外的笼子里。你每次可以进行以下操作：

1.选择两只在不同笼子里的鸽子，使得它们的编号拥有至少一个大于等于P的公共质因数，然后钦定它们为它们当前所在笼子的首领。（否则这两群鸽子就会打起来）

2.将这两个笼子里的所有鸽子赶到同一个笼子里。

反复如上操作，直到剩下的鸽笼都无法合并为止。现在你需要求出最终还有几个鸽笼里有鸽子。

[输入格式]

一行三个正整数 $A,B,P$ 。

[输出格式]

一行一个正整数，表示答案。

[输入样例]

10 20 3

[输出样例]

7

[数据范围]

对于 $100$% 的数据, $1<=A<=B<=100000,2<=P<=B$ 。

---

题意：

有 $A$ ~ $B$ 这 $(B-A+1)$ 个数，现在要将这几个数合并几组，要求组的所有数的公共质因数不小于 $p$ ，求最少能分几组

一看到“合并”就想到了要用并查集

(并查集模板：)

void init () {
    for (int i = 1 ; i <= B ; i ++)  f [i] = i ;
}
int find (int x) {
    return f [x] == x ? x : f [x] = find (f [x]) ;
}
void merge (int x , int y) {
    int fx , fy ;
    fx = find (x) ; fy = find (y) ;
    if (fx != fy)  f [fx] = fy ;
}

那么想想：只要选定一个公共质因数 x ，保证 x >= p ， 那么将它的倍数都合并，不就行了？

for (int i = 2 ; i <= B ; i ++)
        if (! prime [i] && i >= k) {
            sum = 0 ; 
            for (int j = 1 ; j <= B / i ; j ++)
                if (i * j >= A) {
                    sum ++ ; 
                    if (sum == 1)  p = i * j ;
                    else  merge (p , i * j) ;  //merge
                }
        }

最后在计数一下(我喜欢用map，因为map能放一存一，唯一缺点速度较慢)
输出答案

AC代码：

#include <bits/stdc++.h>
#define int long long
#define mod (1e9+7)
#pragma G++ optimize (2)
#pragma G++ optimize (3)
using namespace std ;
int  A , B , k , f [100010] , ans , sum , p , prime [100010] ;
map <int , int> js ;
void init () {
	for (int i = 1 ; i <= B ; i ++)  f [i] = i ;
}
int find (int x) {
	return f [x] == x ? x : f [x] = find (f [x]) ;
}
void merge (int x , int y) {
	int fx , fy ;
	fx = find (x) ; fy = find (y) ;
	if (fx != fy)  f [fx] = fy ;
}
main () {
	ios::sync_with_stdio (false) ;
	cin.tie (NULL) ; cout.tie (NULL) ;
	cin >> A >> B >> k ;
	init () ;
	for (int i = 2 ; i <= B ; i ++)
		if (! prime [i])
			for (int j = i * i ; j <= B ; j += i)
				prime [j] = 1 ;
	for (int i = 2 ; i <= B ; i ++)
		if (! prime [i] && i >= k) {
			sum = 0 ; 
			for (int j = 1 ; j <= B / i ; j ++)
				if (i * j >= A) {
					sum ++ ; 
					if (sum == 1)  p = i * j ;
					else  merge (p , i * j) ;
				}
		}
	for (int i = A ; i <= B ; i ++)
		if (js [find (i)] == 0) {
			ans ++ ; 
			js [find (i)] ++ ;
		}
	cout << ans ;
	return 0 ;
}