欧拉筛三道小练习题

A：

题目：

 输入输出：

样例：

题意：

假设你站在（0，0）处看，求你能看到的可见点个数：例如你在（0，0）能看到（0，1），但是（0，2）被遮挡住所以看不见，为不可见点。

解题方案：

数据范围并不大暴力就行，因为需要特判一下包含1和0的所有点，加起来正好是2n+1，最后加上就行了；

#include<iostream>
#include<string.h>
#include<stdio.h>
#include<math.h>
using namespace std;
bool ma[1010][1010];
int gcd(int a,int b) {
    if(b==0) return a;
    return gcd(b,a%b);
}
void run(int n){
	for(int i=2;i<=n;i++){
		for(int j=2;j<i;j++){
			if(i!=j&&gcd(i,j)==1){
				ma[i][j]=1;
			}
		}
	}
}
int main(){
	int t,len=0;
	run(1000);
	cin>>t;
	while(t--){
		len++;
		int n,ans=0;
		cin>>n;
		ans+=2*n+1;
		for(int i=2;i<=n;i++){
			for(int j=2;j<n;j++){
				if(ma[i][j]==1){
					ans+=2;
				}
			}
		}
		cout<<len<<" "<<n<<" "<<ans<<endl;	
	}
	
}

B：

题目：

 

 输入输出：

 

样例：

题意：

一句话就是让你求n位数组中比ai大的最小非质数的累加，因为数据范围给的并不小，所以直接用欧拉筛筛出质数然后找出最小非质数就行；

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;

bool isPrime[1000010];
//isPrime[i] == 1表示：i是素数
int Prime[1000010], cnt = 0;
//Prime存质数

void euler(int n)//筛到n
{
	memset(isPrime, 1, sizeof(isPrime));
	//以“每个数都是素数”为初始状态，逐个删去
	isPrime[1] = 0;//1不是素数
	
	for(int i = 2; i <= n; i++)
	{
		if(isPrime[i])//没筛掉 
			Prime[++cnt] = i; //i成为下一个素数
			
		for(int j = 1; j <= cnt && i*Prime[j] <= n/*不超上限*/; j++) 
		{
        	//从Prime[1]，即最小质数2开始，逐个枚举已知的质数，并期望Prime[j]是(i*Prime[j])的最小质因数
            //当然，i肯定比Prime[j]大，因为Prime[j]是在i之前得出的
			isPrime[i*Prime[j]] = 0;
            
			if(i % Prime[j] == 0)//i中也含有Prime[j]这个因子
				break; //重要步骤。见原理
		}
	}
}
int main(){
	int n,len=0;
	scanf("%d",&n);
	euler(1000010);
	while(n--){
		len++;
		int a,b;
		long long ans=0;
		scanf("%d",&a);
		for(int i=0;i<a;i++){
			cin>>b;
			for(int j=b+1;;j++){
				if(isPrime[j]){
					ans+=j;
					break;
				}
			}
		}
		cout<<"Case "<<len<<": "<<ans<<" Xukha"<<endl;
	}
}

C:

题目：

 

 输入输出：

样例：

 

题意：

求第k小的质数，值得注意的是范围巨大，非得用欧拉筛不可；

解题方案：

先求1到n中所有的质数，然后用Prime数组存起，再输出Prime[i]就行，询问太多不建议使用cin；

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;

bool isPrime[100000010];
//isPrime[i] == 1表示：i是素数
int Prime[6000010], cnt = 0;
//Prime存质数

void euler(int n)//筛到n
{
	memset(isPrime, 1, sizeof(isPrime));
	//以“每个数都是素数”为初始状态，逐个删去
	isPrime[1] = 0;//1不是素数
	for(int i = 2; i <= n; i++)
	{
		if(isPrime[i])//没筛掉 
			Prime[++cnt] = i; //i成为下一个素数
		for(int j = 1; j <= cnt && i*Prime[j] <= n/*不超上限*/; j++) 
		{
        	//从Prime[1]，即最小质数2开始，逐个枚举已知的质数，并期望Prime[j]是(i*Prime[j])的最小质因数
            //当然，i肯定比Prime[j]大，因为Prime[j]是在i之前得出的
			isPrime[i*Prime[j]] = 0;
			if(i % Prime[j] == 0)//i中也含有Prime[j]这个因子
				break; //重要步骤。见原理
		}
	}
}
int main(){
	int n,q;
	scanf("%d%d",&n,&q);
	euler(n);
	while(q--){
		int a;
		scanf("%d",&a);
		printf("%d\n",Prime[a]);
	}
}