首页 > 代码库 > bzoj 1053: [HAOI2007]反素数ant 搜索
bzoj 1053: [HAOI2007]反素数ant 搜索
1053: [HAOI2007]反素数ant
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1497 Solved: 821
[Submit][Status]
Description
对于任何正整数x，其约数的个数记作g(x)。例如g(1)=1、g(6)=4。
如果某个正整数x满足：g(x)>g(i) 0<i<x，则称x为反质数。例如，整数1，2，4，6等都是反质数。
现在给定一个数N，你能求出不超过N的最大的反质数么？
Input
一个数N（1<=N<=2,000,000,000）。
Output
不超过N的最大的反质数。
Sample Input
1000
Sample Output
840
HINT
Source
這種題應該是通過看數據範圍和解的密度可以發現這是打表，而實際證明連打表都不用，直接暴力構造搜索就行了，搜索中避免如下數字：質因數不連續，質因數過大。這樣解的空間就很小了。
#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<ctime>#include<cmath>#include<algorithm>#include<set>#include<map>#include<vector>#include<string>#include<queue>using namespace std;#ifdef WIN32#define LL "%I64d"#else#define LL "%lld"#endif#define MAXN 11000000#define MAXV MAXN*2#define MAXE MAXV*2#define INF 0x3f3f3f3f#define INFL 0x3f3f3f3f3f3f3f3fLL#define MAXL 2100000000typedef long long qword;inline int nextInt(){ char ch; int x=0; bool flag=false; do ch=getchar(),flag=(ch==‘-‘)?true:flag; while(ch<‘0‘||ch>‘9‘); do x=x*10+ch-‘0‘; while (ch=getchar(),ch<=‘9‘ && ch>=‘0‘); return x*(flag?-1:1);}int n,m;int seq[MAXN],tops=-1;int prime[30]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,51,53,59,61,67,71,73,79,83,89,93,97,101,103,107};struct aaa{ int v,x;}a[MAXN];bool cmp_x(aaa a1,aaa a2){ return a1.x<a2.x;}int topa=-1;void search_m(qword now,int l,int tot){ int i; int x; if (l==30 || now>=MAXL)return ;// cout<<now<<" "<<tot<<endl; a[++topa].x=now; a[topa].v=tot; l++; x=1; for (i=1;i<31 && now<MAXL;i++) { now*=prime[l]; x++; search_m(now,l,tot*x); }}int main(){ freopen("input.txt","r",stdin); //freopen("output.txt","w",stdout); int i,j,k; int x,y,z; int tot; int mx=0; scanf("%d",&n); search_m(1,-1,1); sort(a,a+topa+1,cmp_x); for (i=1;i<=topa;i++) { if (a[i].v>mx) { mx=a[i].v; if (a[i].x>n) { cout<<x<<endl; return 0; } x=a[i].x; } } return 0;/* x=1; mx=0; for (i=1;i<=n;i++) { tot=0; for (j=1;j<=i;j++) { if (i%j==0)tot++; } if (tot>mx) { mx=tot; cout<<i<<" "<<tot<<" "<<endl; x=i; } } cout<<endl;*/ return 0;}
bzoj 1053: [HAOI2007]反素数ant 搜索
声明：以上内容来自用户投稿及互联网公开渠道收集整理发布，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任，若内容有误或涉及侵权可进行投诉： 投诉/举报 工作人员会在5个工作日内联系你，一经查实，本站将立刻删除涉嫌侵权内容。