JZOJ 3504. 运算符(calc)

题目

Description

鸡腿想到了一个很高（sha）明（bi）的运算符，那就是’!’，没错就是感叹号。他给了如下的定义：

1、n!k = n!(k-1) * (n-1)!k (n> 0 and k > 0)

2、n!k = 1 (n = 0)

3、n!k = n (k = 0)

现在鸡腿告诉你n和k你能告诉他n!k的不同约数个数有多少个吗？只要对1,000,000,009取模就可以了哦！

 

Input

一行，输入两个正整数n，k。

Output

一行，输出一个整数表示答案。

 

Sample Input

输入1：

3 1

输入2：

100 2

Sample Output

输出1：

4

输出2：

321266186

 

Data Constraint

对于30%的数据0 <n ≤ 10, 0 <k ≤ 10；

对于100%的数据0 <n ≤ 1000, 0 <k ≤ 100。

 

分析

 

• 这道题显然我们需要先找规律
• 这里一个矩阵，代表
• 1 1 1 1 1
• 2 2 2 2 2
• 6 12 24 48
• 24 288....
• 所以当前数等于左边乘上面
• 左边一列是阶乘
• 然后我们要求的是约数个数
• 就要分解质因数
• 我们举个例子，24=4*6； 4=2^2,6=2*3 24=2^3*3
• 所以我们知道两个数相乘，质因子是相加的
• 所以设f[i][j][k]为到第i行j列第k个质因子转移显然

 

代码

#include<iostream>
using namespace std;
const long long mod=1000000009;
int f[1001][101][170],prime[1001],cnt;
bool vis[1001];
int n,kk;
int main ()
{
    for (int i=2;i<=1000;i++)
      if (!vis[i])
      {
           prime[++cnt]=i;
           for (int j=i+i;j<=1000;j+=i)
           vis[j]=1;
      }
    cin>>n>>kk;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=cnt;j++)
        {
            long long now=prime[j],t=i/now,t1=0;
            if (t==0) break;
            while (t) t1=(t1+t)%mod,t/=now;
            f[i][1][j]=t1%mod;
        }
    for (int i=2;i<=n;i++)
      for (int j=2;j<=kk;j++)
        for (int k=1;k<=cnt;k++)
          f[i][j][k]=(f[i-1][j][k]+f[i][j-1][k])%mod;
    long long ans=1;
    for (int k=1;k<=cnt;k++)
       ans=ans*(f[n][kk][k]+1)%mod;
    cout<<ans;
}