牛客网暑期ACM多校训练营（第三场） H-diff-prime pairs

最新推荐文章于 2018-09-07 22:23:00 发布

原创最新推荐文章于 2018-09-07 22:23:00 发布 · 335 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

ACM

当前文章被收录于：

ACM

56 篇文章 0 人学习

订阅专栏查看详情

当前文章被以下社区和专栏收录：

本文介绍了一种利用素数筛法结合扫描线算法解决特定数学问题的方法。通过预处理得到所有小于等于Max的素数，并统计每个素数倍数的相关数据，最后计算出特定条件下整数对的数量。

//扫描线做法
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define  Max  int(1e7)
bool com[Max+10];
ll primes, prime[Max/10+100];
ll  ans[Max+100];
void solve() 
{ 
    primes = 0; 
    memset(com,false,sizeof(com)); 
    com[0] = com[1] = true; 
    for (ll i = 2; i <= Max; ++i) 
    { 
        if (!com[i]) 
            prime[++primes] = i; 
        for (ll j = 1; j <= primes && i*prime[j] <= Max; ++j) 
        { 
            com[i*prime[j]] = true; 
            if (!(i % prime[j])) 
                break; 
        } 
    } 
}
int main() {
    ll n;
    solve();
    scanf("%lld",&n);
    ans[1]=0;
 
        for(ll i=1;i<=primes&&prime[i]<=n;i++)
           {
           	ans[n/prime[i]+1]-=1;//用来后面求每一个倍数里面有多少个数
           	ans[1]++;
		   }   
		  ll num=0;  
		for(ll i=1;i<=n;i++)
		   {
		   	ans[i]=ans[i]+ans[i-1];
		   	num+=1ll*ans[i]*(ans[i]-1)/2;
		   }   
    	printf("%lld\n",2*num);
	
	return 0;
}