[BZOJ3733]Iloczyn(搜索)

玄学分数算法：贪心：每次尽量选取素因子个数少的数，用bfs实现，但是WA82（好像挺高的）。

刚开始以为是细节出错，所以...

对拍了100000组极限数据也没错

最后在skyh的提示下用小素数相乘来造数据对拍出了一组样例如下：

1

164025000 10

经过手模后发现我的算法会提前用掉以后非常关键的2导致后面的数无法区分。

100分做法：

搜索，但需要强大的剪枝：

设f[x][y]为从x位向后选y个约数的最小乘积，在搜索过程中判断当前剩下的数小于f则return false。

之后就可以得到WA82的好成绩...

我们通过一些小数据可以发现：f的增长是急速的，很可能爆int导致答案不准，

因为我们并不需要f的准确值，所以在处理f数组时把f与n取min即可。

+67...又一次创了我的历史新高...

这道题从开始思考到AC竟然用了4天？！

#include<bits/stdc++.h>
#define AA puts("ALITA")
#define DD puts("DYBALA")
using namespace std;
const int N=1e5+10;
int cnt,ans,n,k,T,p[N],f[N][25];
int read()
{
        register int sum,k=1;register char s;
        while(s=getchar(),s<'0'||s>'9') if(s=='-') k=-1;sum=s-'0';
        while(s=getchar(),s>='0'&&s<='9') sum=sum*10+s-'0';
        return k*sum;
}
bool dfs(int x,int y,int z)
{
        if(!y) return true;
        if(x<p[z]||y>cnt-z+1||f[z][y]>x) return false;
        if(x%p[z]==0&&dfs(x/p[z],y-1,z+1)) return true;
        if(dfs(x,y,z+1)) return true;
        return false;
}
void init(int x)
{
        cnt=0;
        for(int i=1;i<=sqrt(x);i++)
        {
                if(x%i) continue;
                p[++cnt]=i;
                if(i*i!=x) p[++cnt]=x/i;
        }
        sort(p+1,p+cnt+1);
        f[cnt+1][0]=1;
        for(int i=cnt;i>=1;i--)
        {
                f[i][0]=1;
                for(int j=1;j<=min(k,cnt-i+1);j++)
                {
                        f[i][j]=min(n,f[i+1][j-1]*p[i]);
                }
        }
}
bool work()
{
        n=read();k=read();
        init(n);
        return dfs(n,k,1);
}
signed main()
{
        //freopen("1.in","r",stdin);
        //freopen("1.out","w",stdout);
        T=read();
        while(T--)
        {
                if(work()) puts("TAK");
                else puts("NIE");
        }
        return 0;
}