首页 > 代码库 > 51nod 1434 区间LCM (质因数分解)
51nod 1434 区间LCM (质因数分解)
分析:考虑从1到n所有数的质因数分解,记录每个质数的最高次数,同理从n+1循环到2n,如果循环到m时每个质因子的次数都不低于所记录的,则跳出循环,结果即为m。先预处理质数,复杂度为O(nlongn)。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm> 5 using namespace std; 6 const int maxn=1000005; 7 int prinum[maxn],len=0,n,m,power[maxn],order[maxn]; 8 int num[maxn]; 9 void CalPri(){ 10 memset(order,-1,sizeof(order)); 11 for(int i=0;i<maxn;i++)num[i]=i; 12 for(int i=2;i<maxn;i++){ 13 if(num[i]==0)continue; 14 prinum[len++]=i; 15 order[i]=len-1; 16 for(int j=2*i;j<maxn;j+=i){ 17 num[j]=0; 18 } 19 } 20 } 21 int solve(){ 22 int countp=0; 23 if(order[n]!=-1)return 2*n; 24 memset(power,0,sizeof(power)); 25 for(int k=0;k<len&&prinum[k]<=n;k++){ 26 long long p=prinum[k],p0=prinum[k]; 27 countp++; 28 while(p<=n){ 29 p*=p0; 30 power[k]++; 31 } 32 // int k0=k; 33 // for(int i=0;i<len&&prinum[i]*prinum[i]<=k0;i++){ 34 // if(k0%prinum[i])continue; 35 // int p=prinum[i],cnt=0; 36 // if(power[i]==0)countp++; 37 // while(k0%p==0){ 38 // k0/=p;cnt++; 39 // } 40 // power[i]=max(cnt,power[i]); 41 // } 42 // if(k0!=1){ 43 // if(power[order[k0]]==0)countp++; 44 // power[order[k0]]=max(1,power[order[k0]]); 45 // } 46 } 47 for(int k=n+1;k<=2*n;k++){ 48 int k0=k; 49 for(int i=0;i<len&&prinum[i]*prinum[i]<=k0;i++){ 50 if(k0%prinum[i])continue; 51 int p=prinum[i],cnt=0; 52 while(k0%p==0){ 53 k0/=p;cnt++; 54 } 55 if(power[i]>0&&power[i]<=cnt){ 56 countp--; 57 power[i]=0; 58 } 59 } 60 if(k0!=1&&power[order[k0]]==1){ 61 countp--;power[order[k0]]=0; 62 } 63 if(!countp){ 64 return k; 65 } 66 } 67 } 68 int main(){ 69 //freopen("e:\\in.txt","w",stdout); 70 CalPri(); 71 int t; 72 cin>>t; 73 while(t--){ 74 cin>>n; 75 cout<<solve()<<endl; 76 } 77 return 0; 78 }
51nod 1434 区间LCM (质因数分解)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。