当前位置: 代码迷 >> 综合 >> 【筛法】Codeforces1017F The Neutral Zone
  详细解决方案

【筛法】Codeforces1017F The Neutral Zone

热度:59   发布时间:2023-09-27 06:50:28.0

分析:

就是个算贡献的水题。。。但出题人相当无聊。。。把内存卡得爆小。。。只有16M。。。。

不得不说CF的评测机真心快。。。n=3?108n=3?108的情况下nlognnlogn居然还能卡进2s(卧槽不科学啊啊)

其实就是个埃氏筛,但是由于n=3?108n=3?108且内存只有16M,bitset也只能存大概108108位,出题人以为这样就能卡人了?太天真了。。。直接对[0,108],[108+1,2?108],[2?108+1,3?108][0,108],[108+1,2?108],[2?108+1,3?108]分别做一次,每次只用108108的内存。。。就能过。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<bitset>
#include<cmath>
#include<vector>
#define SF scanf
#define PF printf
#define MAXN 100000010
#define MAXM 30010
using namespace std;
typedef unsigned int ui;
ui n,a,b,c,d,ans,m;
bitset<MAXN> used;
bitset<MAXM> prep;
ui px(ui x){return a*x*x*x+b*x*x+c*x+d;
}   
void get_ans(ui x){ui x1=x;if(x==1)return ;//PF("[%u %u]",x,ans);while(x<=n){ans+=px(x1)*(n/x);if(1ll*x*x1>1ll*n)return ;x=x*x1;}
}
int main(){SF("%u%u%u%u%u",&n,&a,&b,&c,&d);m=ceil(sqrt(n)+0.5);prep[0]=1;for(ui i=2;i<=m;i++)if(prep[i]==0){for(ui j=2;j*i<=m;j++)prep[i*j]=1;}for(ui up=0;n>up;up+=100000000){used.set();used.flip();for(ui i=2;i<=m;i++)if(prep[i]==0)for(ui j=max(2u,(up+i-1)/i);i*j<=up+100000000;j++)used[i*j-up]=1;for(ui i=1;i<=100000000&&i+up<=n;i++)if(used[i]==0)get_ans(i+up);}PF("%u",ans);
}
  相关解决方案