codeforces 920g List Of Integers

      题意：求大于x的与p互质的第k个数。

      首先我们可以二分将题目转换为判定性问题，即二分一个mid，判断大于x小于等于mid与p互质有几个数字，那么我们可以将问题转化为，小于等于mid的与p互质的数字数量减小于等于x的与p互质的数字数量。

      那么如何求解小于等于x的与p互质数字数量呢，由于x小于等于10的6次方，所以质因子个数很少，所以我们可以听过容斥原理n-有一个质因子的数的个数+有2个质因子数的个数-...，而求x以内能整除y的数的个数就是x/y了。

      那么我们就通过了二分+容斥解决了这个问题。

      下附AC代码。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define maxn 1000005
using namespace std;
typedef long long ll;
int tot;
ll a,p,k;
int pri[maxn],vis[maxn];
void shai()
{
	for(int i=2;i<=1000000;i++)
	{
		if(!vis[i])
			pri[++tot]=i;
		for(int j=1;j<=tot && i*pri[j]<=1000000;j++)
		{
			vis[i*pri[j]]=1;
			if(i%pri[j]==0)
			break;
		}
	}
}
int cnt=0;
ll mypri[maxn];
void getpri(ll now)
{
	for(int i=1;i<=tot && now!=1;i++)
	{
		if(now%pri[i]==0)
		{
			mypri[++cnt]=pri[i];
			while(now%pri[i]==0)
			now/=pri[i];
		}
	}
	return;
}
ll ans=0;
void dfs(ll pos,ll nowcnt,ll nowsum,ll maxx)
{
	if(nowsum>maxx) return;
	if(pos==cnt+1)
	{
		if(!nowcnt) return;
		ll flag=(nowcnt&1) ? 1 : -1;
		ans+=flag*(maxx/nowsum);
		return;
	}
	dfs(pos+1,nowcnt,nowsum,maxx);
	dfs(pos+1,nowcnt+1,nowsum*mypri[pos],maxx);
}
ll solve(ll now)
{
	ans=0;
	dfs(1,0,1,now);
	return now-ans;
}
int check(ll now,ll temp)
{
	return (solve(now)-temp>=k);
}
int main()
{
	shai();
	int _;
	scanf("%d",&_);
	while(_--)
	{
		scanf("%I64d%I64d%I64d",&a,&p,&k);
		cnt=0;
		getpri(p);
		ll temp=solve(a);
		ll l=a+1,r=12345678910,res;
		while(l<=r)
		{
			ll mid=((l+r)>>1);
			if(check(mid,temp))
				res=mid,r=mid-1;
			else
				l=mid+1;
		}
		printf("%I64d\n",res);
	}
}