【Test 2016-10-2】antiprime {搜索+剪枝+高精}

感觉自己搜索简直爆炸 >-<

---

• 【题目描述】
  将正整数x 的约数个数表示为g(x)。例如，g(1) = 1; g(4) = 3; g(6) = 4。 如果对于任意正整数y，当0 < y < x 时，x 都满足g(x) > g(y), 则称x 为反质数。整数1，2，4，6等都是反质数

• 【Sample Input】（一个数n）
  sample 1
  13
  sample 2
  100

• 【Sample Output】（一个数表示小于等于N的最大的反质数）
  sample 1
  12
  sample 2
  840

• 【数据范围】
  对于40% 的数据，N<=10^9
  对于70% 的数据，N<=10^25
  对于100% 的数据，N<=10^100

---

【题解】搜索啊搜索（+高精

首先答案的因数一定是从二开始一连串的素数。
一开始很容易贪心的想我们尽量乘不同的素数越多数越好呀！后来很伤心的发现自己错了（随意举个例子可以推翻）。不过这样一来我们可以算出最大素数的上界，也就是从2开始到250+的素数相乘就会超过10^100。所以筛素数的范围只要到260。
然后爆搜就能过70%了 QAQ（是我蠢啊蠢

{剪枝一} 搜索时枚举因数，对于一个因数y（y>2），设此时枚举2的指数为p（2^p），y的指数为q（y^q）若：2^x<=y<2^(x+1)，
则：q<(p+1)/x

[证明]：设当前枚举数为 M=2^p * a1^b1 * a2^b2 * …* y^q * …
所以 g(M)=(p+1) * (b1+1) * (b2+1) * … * (q+1) *…
这时，我们用 2^x 去替换y，
构造数 m=2^(p+x) * a1^b1 * a2^b2 * …* y^(q-1) * …
那么 m 的因数个数 g(m)=(p+x+1) * (b1+1) * (b2+1) * … * q *…
显然，m< M。所以如果 M 优于m，必然有g(M)>g(m)。
即：(p+1)*(q+1)>(p+x+1)q
化简得：p+1>x*q → ：q<(p+1)/x
//感觉出题人机智到爆炸

{剪枝二}只有出题人才懂的优化
2的指数最大是14… 我不知道为什么我不知道为什么我不知道为什么
//大概是因为如果写15TLE，写13WA吧

---

#include <cstdio>
#include <iostream>
#include <cstring>
#define LL long long
#define bit 10000
struct num{ int l,a[30];};
num n,ans,q;int k,tot,p[300],a[300];
LL anss;
bool bo[300];
    num operator*(num x,int y)
    {
        int i;
        for (i=1;i<=x.l;++i) x.a[i]*=y;
        for (i=1;i<=x.l;++i)
            if (x.a[i]>=bit)
            {
                x.a[i+1]+=x.a[i]/bit;x.a[i]%=bit;
            }
        for (;x.a[x.l+1];++x.l);
        return x;
    }
    bool operator<(num x,num y)
    {
        if (x.l<y.l) return true;
        if (x.l>y.l) return false;
        for (int i=x.l;i;--i)
            if (x.a[i]<y.a[i]) return true;
            else if (y.a[i]<x.a[i]) return false;
        return true;
    }
    void dfs(int x,LL y,int z,num w)
    {
        if (y>anss || w<ans && y==anss) ans=w,anss=y;
        num t=w;
        if (x>1) z=std::min(z,k/a[x]);
        for (int i=1;i<=z;++i)
        {
            if (x==1) k=i;
            t=t*p[x];
            if (t<n) dfs(x+1,y*(i+1),i,t);
            else break;
        }
    }
    void print(num x)
    {
        printf("%d",x.a[x.l]);
        for (int i=--x.l;i;--i) printf("%.4d",x.a[i]);
        printf("\n");
    }
int main()
{
    char s[105];scanf("%s\n",s+1);
    int len=strlen(s+1);
    for (int i=len,j;i>0;i-=4)
    {
        ++n.l;
        if (i>=4) j=i-3;else j=1;
        for (;j<=i;++j) n.a[n.l]=n.a[n.l]*10+s[j]-'0';
    }
    memset(bo,true,sizeof(bo));
    for (int i=2;i<=260;++i)
        if (bo[i])
            for (int j=i+i;j<=260;j+=i) bo[j]=false;
    for (int i=2;i<=260;++i)
        if (bo[i]) p[++tot]=i;
    for (int i=1,j=0;i<=tot;a[i++]=j-1)
        for (;(1<<j)<=p[i];++j);
    ans.l=1;ans.a[1]=0;
    q.l=1;q.a[1]=1;
    dfs(1,1,14,q);print(ans);
    return 0;
}

---

【题外话】

我感觉自己真是弱到搜索都写不好，写到TLE无数发。

%%% aufeas 给我讲题；%%% hjj 帮我看程序。