洛谷3861八月月赛A题解

链接

用f[i][j]表示乘积为i的，包含的最大数小于等于j时的方案总数

我们先考虑所用的数为1到n的情况

最后的答案就是f[n][n]-1

转移时考虑f[i][j]可以转移到的状态

显然f[i][j]可以转移到f[i*k][k](k>j)，可以转移到f[i][k]（k>j）

我们在考虑，可以发现只有i,j都为n的因数时才会对答案产生贡献

而1e12中因数最多的数有大约7000个因数，所以我们就可以愉快的O（7000^2）水过去了

# include<iostream>
# include<cstdio>
# include<cstring>
# include<cmath>
# include<algorithm>
using namespace std;
typedef long long LL;
const int mod = ;
const int mn = ;
LL f[mn][mn],a[mn],n;
int cnt,t;
void work(LL x)
{
    LL m=sqrt(x*1.0);
    for(int i=;i<=m;i++)
    {
        if(x%i==)
        {
            if(1ll*i*i==x)
                a[++cnt]=i;
            else {
                a[++cnt]=i;
                a[++cnt]=x/i;
            }
        }
    }
    sort(a+,a++cnt);
    a[cnt+]=x+;//增加虚拟节点
    f[][]=;
    //printf("%d",cnt);
    for(int i=;i<=cnt;i++)
    {
        int l=i+;
        for(int j=;j<cnt;j++)
        {
            if(f[i][j]==) continue;
            if(a[i]*a[j+]<=x)
            {
                while(a[l+]<=a[i]*a[j+]) l++;
                if(a[l]==a[i]*a[j+])
                    f[l][j+]=(f[l][j+]+f[i][j])%mod;
            }
            f[i][j+]=(f[i][j+]+f[i][j])%mod;
        }
    }
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        cnt=;
        //memset(f,0,sizeof(f));
        scanf("%lld",&n);
        work(n);
        printf("%lld\n",f[cnt][cnt]-);
        if(t)
        {
            for(int i=;i<=cnt;i++)
                for(int j=;j<=cnt;j++)
                    f[i][j]=;
        }
    }
    return ;
}