bzoj2793: [Poi2012]Vouchers

最新推荐文章于 2018-09-05 12:02:00 发布

转载最新推荐文章于 2018-09-05 12:02:00 发布 · 116 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/Sakits/p/7744027.html

本文介绍了一个筛选特定范围内幸运数的算法实现。通过记录每个幸运数倍数的位置，并使用调和级数来降低复杂度至O(NlogN)，该算法有效地找出所有符合条件的幸运数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

　　傻逼题...

　　直接记录x倍数已经扫到哪里了，每次接着上次的扫就好，因为幸运数范围只有100w，扫到100w就停，根据调和级数复杂度为O(NlogN)。

　　记得开long long！！！T T

#include<iostream> 
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath> 
#include<algorithm> 
#define ll long long
using namespace std;
const int maxn=1000010, inf=1e9;
int n, m, x, cnt;
int now[maxn];
ll num, ans[maxn];
bool lucky[maxn], used[maxn];
inline void read(int &k)
{
    int f=1; k=0; char c=getchar();
    while(c<'0' || c>'9') c=='-'&&(f=-1), c=getchar();
    while(c<='9' && c>='0') k=k*10+c-'0', c=getchar();
    k*=f;
}
int main()
{
    read(n);
    for(int i=1;i<=n;i++) read(x), lucky[x]=1; 
    for(int i=1;i<=1000000;i++) now[i]=i;
    read(m); num=0;
    for(int i=1;i<=m;i++)
    {
        read(x);
        for(int j=1;j<=x;j++)
        {
            for(int k=now[x];k<=1000000;k+=x, now[x]=k)
                if(!used[k]) break;
            if(now[x]>1000000)break;
            used[now[x]]=1;
            if(lucky[now[x]]) ans[++cnt]=num+j;
        }
        num+=x;
    }    
    printf("%d\n", cnt);
    for(int i=1;i<=cnt;i++) printf("%lld\n", ans[i]); 
}