题意:给出多个n(1<=n<=1e6),求满足phi(x)>=n的最小的x之和。
分析:先预处理出1~1e6的欧拉函数,然后建立一颗线段树维护最大值,对于每个n询问大于等于n的最左边下标。
#pragma comment(linker,"/STACK:1024000000,1024000000") #include <cstdio> #include <cstring> #include <string> #include <cmath> #include <limits.h> #include <iostream> #include <algorithm> #include <queue> #include <cstdlib> #include <stack> #include <vector> #include <set> #include <map> #define LL long long #define mod 100000000 #define inf 0x3f3f3f3f #define eps 1e-6 #define N 1001000 #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 #define PII pair<int,int> using namespace std; inline LL read() { char ch=getchar();LL x=0,f=1; while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch<=‘9‘&&ch>=‘0‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f; } int phi[N+5],prime[N/10],mx[N<<2]; void init() { for(int i=2;i<=N;i++)phi[i]=i; int tot=0; for(int i=2;i<=N;i++) { if(phi[i]==i) { for(int j=i;j<=N;j+=i) phi[j]=phi[j]/i*(i-1); } } } void Pushup(int rt) { int ls=rt<<1,rs=ls|1; mx[rt]=max(mx[ls],mx[rs]); } void build(int l,int r,int rt) { if(l==r) { mx[rt]=phi[l]; return; } int m=(l+r)>>1; build(lson); build(rson); Pushup(rt); } int query(int x,int l,int r,int rt) { if(l==r)return l; int m=(l+r)>>1; if(mx[rt<<1]>=x)return query(x,lson); else return query(x,rson); } int main() { int T,n,cas=1; init(); build(1,N,1); T=read(); while(T--) { n=read(); LL ans=0; for(int i=1;i<=n;i++) { int x=read(); ans+=query(x,1,N,1); } printf("Case %d: %lld Xukha\n",cas++,ans); } }
原文:http://www.cnblogs.com/lienus/p/4298925.html