bzoj 1053: [HAOI2007]反素数ant 搜索

1053: [HAOI2007]反素数ant

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 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 2100000000
typedef 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;
}