线性筛、gcd+题目

线性筛的时间复杂度为O(n)

题目链接：质数问题

/**/
#include<bits/stdc++.h>
using namespace std;
int t,n,k,cnt,res;
bool st[1005];//1表示是合数，0表示质数
int prime[1005];//prime[]数组存的是质因子 
void getprime(int x){
	cnt=0;
	//memset(x,1,sizeof(x));//假定全部都是合数 
	st[0]=st[1]=0;//0,1既不是合数也不是质数
	for(int i=2;i<=n;i++){
		if(st[i]==0){
			cnt++;
			prime[cnt]=i;
		} 
		for(int j=1;j<=cnt&&i*prime[j]<=n;j++){
			st[i*prime[j]]=1;
			if(i%prime[j]==0) break;
			//保证一个数只被筛选一次，例如12是在i=6时被筛掉的，而不是在i=4，4*3筛掉的 
		}
	} 
}

int main(){
	cin>>t;
	while(t--){
		res=0; 
		cin>>n>>k;
		getprime(n);
		for(int i=2;i<=n;i++){
			if(st[i]==0){
				for(int j=1;j<=cnt;j++){
					if(prime[j]+prime[j+1]+1==i){
						res++;
						break;
					}
				}
			}
		} 
		if(res>=k) cout<<"YES"<<endl;
		else cout<<"NO"<<endl;
	}
}
 

题目链接：消灭老鼠

/**/
#include<bits/stdc++.h>
using namespace std;
int n,x,y;
set<pair<int,int>>s;
int gcd(int a,int b){
	if(b==0) return a;
	return gcd(b,a%b);
}
int main(){
	cin>>n>>x>>y;
	int ans=0;
	while(n--){
		int x1,y1;
		cin>>x1>>y1;
		if(x1==x) ans=1;//斜率不存在的情况 
		else{
			int yy=y1-y,xx=x1-x;
			double k1=(y1-y)/gcd(yy,xx);//把斜率换为最简形式 
			double k2=(x1-x)/gcd(yy,xx);
			s.insert({k1,k2});//用点对存斜率更方便 
		}
	}
	cout<<s.size()+ans;
}

快速幂的核心思想：每次把指数分为两半，底数做平方运算

eg.3^10

(3*3)^5=9^5

题目链接：幂次方

/*快速幂的实质上是倍增*/
#include<bits/stdc++.h>
using namespace std;
long long n,x;//要用 long long ，不然会WA 
long long solve(long long p,long long b){
	long long res=1;
	while(p){
		if(p%2) res=res*b%233333;
		p/=2;
		b=(b*b)%233333;
	}
	return res;
}
int main(){
	cin>>x>>n;
	cout<<solve(n,x);
}