poj 2773 容斥原理+二分(求与n的第k互质数)

点击打开链接

#include <iostream>
#include <cstring>
#include <vector>
using namespace std;
typedef long long ll;
const int N=1e6+20;
const ll inf=1e12;
ll n,k,num,ans;
int fac[N];
void Get(ll n)
{
	for(int i=2;i*i<=n;i++)
	{
		if(n%i==0)
		{
			fac[num++]=i;
			while(n%i==0)
			n/=i;
		}
	}
	if(n>1)
	fac[num++]=n;
}
ll dfs(int id,ll m)//容斥原理算1~m与n不互质的个数 
{
	ll cnt=0;
	for(int i=id;i<num;i++)
	{
		cnt+=m/fac[i]-dfs(i+1,m/fac[i]);
	}
	return cnt;
}
int main()
{
	while(cin>>n>>k)
	{
		num=0;
		Get(n);
		ll l=1,r=inf;
		while(l<=r)//1~x1中与n互质个数小于1~x2中与n互质个数 x1<x2 满足单调性用二分求第k个 
		{
			
			ll m=(l+r)/2;
			ll cnt=m-dfs(0,m);//求1~m中有多少个数与n互质 
			if(cnt>=k)
			{
				r=m-1;
				if(cnt==k)
				ans=m;//记录最左边成立的m 
			}
			else 
			{
				l=m+1;
			}
			
		}
		cout<<ans<<endl;
	}
	return 0;
}