cf #473 div2 D(素筛)

题意:
给n个数,求一个字典序最小,要求这个序不能小于原序列,并保证两两互素
题解:
思路很简单,首先埃筛出1-1e7的每个数的最小素因子作为预处理(但是题目给的是1e5,1e7不浪费时间空间吗?后面会讲到)。
现在就是要找第一个与原序列不同的数,在这个数之后,后面的所有数都搞成未标记素数,以保证整体字典序最小。
把每个数的素因子查看,是不是已存在,所有素因子查看完毕并未有标记的可以保持原序继续寻找第一个关键点,否则找最小未标记素数放进去
关键点找到后,直接往后面塞素数,但是考虑最坏情况1e5个素数最大的至少是上了1e7,所以我才上封1e7
code:(写得有点乱)

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=1e7+5;
int a[maxn],b[maxn],c[maxn],d[maxn],vis[maxn],p[maxn/10],tot=0;
void setviso(int x){
    for(int i=x;i<maxn;i+=x)vis[i]=1;
}
int find_index(int x){
    while(vis[p[x]]&&x<tot)x++;
    return x;
}
void setvis(int x){
    for(int i=d[x];x!=1;x/=i,i=d[x])vis[i]=1,setviso(i);
}
int find_ok(int x){
    while(vis[x])x++;
    return x;
}
int main(){
    bool flag=false;
    int lastnum=0;
    int n;
    set<int>pc;
	scanf("%d",&n);
	for(int i=2;i<maxn;i++)//埃筛最小因子
        if(!c[i]){
            d[i]=i;p[tot++]=i;//埃筛质数
            for(int j=i+i;j<maxn;j+=i){
                c[j]=1;
                if(!d[j])d[j]=i;
            }
        }
	for(int i=1;i<=n;i++)scanf("%d",a+i);
	for(int i=1;i<=n;i++){
        int elem=a[i];
        if(flag){
            b[i]=p[find_index(lastnum)];
            vis[b[i]]=1;
            lastnum=find_index(lastnum);
            continue;
        }
        //pc存所有因子,判断是否标记
        pc.clear();
        for(int i=d[elem];elem!=1;elem/=i,i=d[elem])
            pc.insert(i);
        for(auto x:pc){
            if(vis[x]){//被标记去找最小素数
                flag=true;
                b[i]=find_ok(a[i]);
                break;
            }
        }
        if(!b[i])b[i]=a[i];
        setvis(b[i]);
    }
	for(int i=1;i<=n;i++)
        printf("%d ",b[i]);
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Macarons_i

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值