标签：tmp PAT Progression int anstmp 质数 vis 2021 maxd

#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
const int maxn = 1e5+10;
int n,m,primes[maxn],num,vis[maxn];
void primes_table() { //质数表
	vis[0]=1,vis[1]=1;
	for(int i=2; i<=m; i++) {
		if(vis[i]!=0)continue;
		primes[num++]=i;
		for(int j=i+i; j<=m; j+=i) {
			vis[j]=1;
		}
	}
}
int main() {
	scanf("%d%d",&n,&m);
	//质数表
	primes_table();
	// 找到质数表中小于m的最大质数的位置
	if(num<n)printf("%d\n",primes[num-1]); //没有找到数列，打印最后一个
	else if(num==n) {
		for(int i=0; i<num; i++) {
			if(i!=0)printf(" ");
			printf("%d",primes[i]);
		}
	} else { //num>n 寻找解决方案
		vector<int> tmp,anstmp;
		vector<vector<int>> ans;
		int mmaxd = -1; // 所有符合答案的序列中最大差 
		for(int i=0; i<num-1; i++) {
			int maxd=-1; // 以primes[i]为起点满足条件的序列的最大差 
			for(int j=i+1; j<num; j++) {
				int d=primes[j]-primes[i];
				if(d<maxd||d<mmaxd)continue;//剪枝，如果差小于已经找到的最大差，跳过 
				int c=primes[j];//数列第一个数
				tmp.clear(); //每次查找清楚tmp
				tmp.push_back(primes[i]); 
				while(c<m) { //在范围内
					if(vis[c]==1)break;
					tmp.push_back(c);
					c+=d;
					if(tmp.size()>n) { //如果有多个，选择第一个元素最大的
						tmp.erase(tmp.begin());
						anstmp=tmp;
						maxd=d; //更新已找到的等差质数序列的最大差 
					} else if(tmp.size()==n) {
						maxd=d; //更新已找到的等差质数序列的最大差 
						anstmp=tmp;
					}
				}
			}
			if(maxd>mmaxd)mmaxd=maxd;
		}
		if(anstmp.size()==0) {
			printf("%d",primes[num-1]);
		} else {
			for(int i=0; i<n; i++) {
				if(i!=0)printf(" ");
				printf("%d",anstmp[i]);
			}
		}
	}
	return 0;
}

标签：tmp,PAT,Progression,int,anstmp,质数,vis,2021,maxd
来源： https://www.cnblogs.com/houzm/p/14529667.html

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