标签：prime ... Pairs frac gp int 素数 Diff

牛客 Diff-prime Pairs

题意:给定 N ( 1 ≤ N ≤ 1 0 7 ) N(1\le N \le 10^7) N(1≤N≤107)，求有多少对数 ( i , j ) (i,j) (i,j)满足 i g c d ( i , j ) \frac{i}{gcd(i,j)} gcd(i,j)i​和 j g c d ( i , j ) \frac{j}{gcd(i,j)} gcd(i,j)j​均为素数， ( i 1 , j 1 ) , ( i 2 , j 2 ) (i_1,j_1),(i_2,j_2) (i1​,j1​),(i2​,j2​)不同当且 i 1 ! = i 2   o r   j 1 ! = j 2 i_1!=i_2\ or\ j_1!=j_2 i1​!=i2​ or j1​!=j2​

Sol:

i g = p r i m e \frac{i}{g}=prime gi​=prime j g = p r i m e \frac{j}{g}=prime gj​=prime

i = p r i m e ∗ g i=prime*g i=prime∗g j = p r i m e ∗ g j=prime*g j=prime∗g

i ∈ [ 1 , n ] i\in [1,n] i∈[1,n] j ∈ [ 1 , n ] j\in [1,n] j∈[1,n]

这样的 i , j i,j i,j要么是素数对，要么是素数对的倍数

枚举 g g g，则 n g \frac{n}{g} gn​表示 p r i m e prime prime的范围，如果在 n g \frac{n}{g} gn​范围内有 k k k个素数

则说明有 ∣ ( g p 1 , g p 2 , . . . , g p k ) × ( g p 1 , g p 2 , . . . , g p k ) ∣ − k |(gp_1,gp_2,...,gp_k)\times (gp_1,gp_2,...,gp_k)|-k ∣(gp1​,gp2​,...,gpk​)×(gp1​,gp2​,...,gpk​)∣−k对组合是符合要求的

且 i , j ∈ g p 1 , g p 2 , . . . , g p k i,j\in {gp_1,gp_2,...,gp_k} i,j∈gp1​,gp2​,...,gpk​

所以枚举 g g g可以表示出所有素数对和他们的倍数

#include <bits/stdc++.h>
using namespace std;
const int N=1e7+10;
int prime[N],st[N],sum[N];
int cnt;
void get_primes(int n)
{
     for(int i=2;i<=n;i++)
     {
     	if(!st[i]) prime[++cnt]=i,sum[i]=1;

            for(int j=1;i*prime[j]<=n;j++)
            {
                st[i*prime[j]]=1;
                if(i%prime[j]==0) break;
            }
        
     }
}
int main()
{     
	 ios::sync_with_stdio(0);
	 cin.tie(0);
	 int n;
	 cin>>n;
	 get_primes(n);
  
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+sum[i];
	 long long ans=0;
     ans=1LL*(cnt-1)*cnt;
	 cout<<ans<<endl;
	 return 0;
}

标签：prime,...,Pairs,frac,gp,int,素数,Diff
来源： https://blog.csdn.net/weixin_43900484/article/details/120595324

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