首页 > 代码库 > ZOJ 2562 HDU 4228 反素数
ZOJ 2562 HDU 4228 反素数
反素数:
对于任何正整数x,起约数的个数记做g(x).例如g(1)=1,g(6)=4.
如果某个正整数x满足:对于任意i(0<i<x),都有g(i)<g(x),则称x为反素数.
ZOJ 2562 反素数
因为写了POJ 2886的线段树,然后里面有反素数,以前没遇到过,所以先搞这两题普及一下知识再说。
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<map> #include<queue> #include<set> #include<cmath> #include<bitset> #define mem(a,b) memset(a,b,sizeof(a)) #define lson i<<1,l,mid #define rson i<<1|1,mid+1,r #define INF 510010 #define maxn 400010 using namespace std; typedef long long ll; typedef unsigned long long ull; ll prime[20]={2,3,5,7,11,13,17,19,23,29,31,37,39,41,43,47,53}; ll n; ll bestcurr;//bestcurr 相同最大因数个数中值最小的数 ll largecnt;//largecnt:n范围内最大的因数个数 void getarcprime(ll curr,int cnt,int limit,int k) { if(curr>n) return ; if(largecnt<cnt)//此时枚举到的因数个数比之前记录的最大的因数个数要大,就替换最大因数个数 { largecnt=cnt; bestcurr=curr; } if(largecnt==cnt && bestcurr>curr)//替换最优值 bestcurr=curr; ll temp=curr; for(int i=1;i<=limit;i++) { temp=temp*prime[k]; if(temp>n) return; getarcprime(temp,cnt*(i+1),i,k+1); } } int main() { while(scanf("%lld",&n)!=EOF) { bestcurr=0; largecnt=0; getarcprime(1,1,50,0); printf("%lld\n",bestcurr); } return 0; }HDU 4228
这题就是上题的延伸吧,就是求出每个然后打表。
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<map> #include<queue> #include<set> #include<cmath> #include<bitset> #define mem(a,b) memset(a,b,sizeof(a)) #define lson i<<1,l,mid #define rson i<<1|1,mid+1,r #define INF 510010 #define maxn 400010 using namespace std; typedef long long ll; typedef unsigned long long ull; ll p[1010]; ll prime[30]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53}; void getartprime(ll cur,int cnt,int limit,int k) { //cur:当前枚举到的数; //cnt:该数的因数个数; //limit:因数个数的上限;2^t1*3^t2*5^t3……t1>=t2>=t3…… //第k大的素数 if(cur>(1LL<<60) || cnt>150) return ; if(p[cnt]!=0 && p[cnt]>cur)//当前的因数个数已经记录过且当时记录的数比当前枚举到的数要大,则替换此因数个数下的枚举到的数 p[cnt]=cur; if(p[cnt]==0)//此因数个数的数还没有出现过,则记录 p[cnt]=cur; ll temp=cur; for(int i=1; i<=limit; i++) //枚举数 { temp=temp*prime[k]; if(temp>(1LL<<60)) return; getartprime(temp,cnt*(i+1),i,k+1); } } int main() { int n; getartprime(1,1,75,0); for(int i=1; i<=75; i++) { if(p[i*2-1]!=0 && p[i*2]!=0) p[i]=min(p[i*2-1],p[i*2]); else if(p[i*2]!=0) p[i]=p[i*2]; else p[i]=p[i*2-1]; } while(scanf("%d",&n),n) printf("%I64d\n",p[n]); return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。