个人存档用

YanRC · 2025-11-13 17:19:44 · 个人记录

#include<bits/stdc++.h>
using namespace std;
int cnt;
vector<bool>es(int x)
{
    vector<bool>isp(1+x,true);
    for(int i=2;i*i<=x;i++)if(isp[i])for(int j=i*2;j<=x;j++)isp[j]=false;
    for(int i=1;i<=x;i++)cnt+=isp[i];
    return isp;
}
int main()
{
    cout<<setprecision(15);
    vector<bool>v;
    int list[14]={0,10,15,20,50,100,150,200,500,1000,1500,2000,5000,10000};
    for(int i=1;i<=13;i++)
    {
        cnt=0;
        v=es(list[i]);
        int c=0;
        for(int j=3;j<=list[i];j++)c+=(v[j]&&(v[j-2]||v[j+2]));
        cout<<"N="<<list[i]<<",TPC="<<1.0*c/cnt<<endl;
    }
    return 0;
}