运算符(calc)

最新推荐文章于 2023-10-16 10:27:27 发布

Chandery

最新推荐文章于 2023-10-16 10:27:27 发布

阅读量490

点赞数 1

本文链接：https://blog.csdn.net/cdy1206473601/article/details/54985619

版权

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

题解：首先把1~1000里的所有质数求出来，然后设
f[i,j,k]表示坐标为i，j的第k个质因数的指数，所有，把每个k=0,的（i,0）这个数分解质因数，算出每个质因数的指数，所以每次
f[i,j,k]:=(f[i-1,j,k]+f[i,j-1,k]) mod 1000000009；
所以最后运用公式：x的约数个数=(q1+1)(q2+1)….(qn+1)；
把(f[n,n,i]+1)乘起来。
代码：

        s:array[0..1000]of boolean;
        ss:array[0..168]of longint;
        b:array[0..1000,0..100]of longint;
        f:array[0..1000,0..100,0..168]of longint;
        i,j,n,k,l,t,m:longint;
        ans:int64;
function max(a,b:longint):longint;
begin
        if a>b then exit(a);
        exit(b);
end;
begin
        readln(n,k);
        fillchar(s,sizeof(s),true);
        for i:=2 to 1000 do
                for j:=2 to 1000 div i do
                        s[i*j]:=false;
        s[1]:=false;
        for i:=2 to 1000 do
                if s[i] then
                begin
                        inc(t);
                        ss[t]:=i;
                end;
        for i:=2 to n do
        begin
                j:=1;
                m:=i;
                while m>1 do
                begin
                        while m mod ss[j]=0 do
                        begin
                                inc(f[i,0,j]);
                                m:=m div ss[j];
                        end;
                        inc(j);
                end;
                b[i,0]:=j;
        end;
        for i:=1 to n do
                for j:=1 to k do
                begin
                        for l:=1 to max(b[i-1,j],b[i,j-1]) do
                                f[i,j,l]:=(f[i-1,j,l]+f[i,j-1,l]) mod 1000000009;
                        b[i,j]:=max(b[i-1,j],b[i,j-1]);
                end;
        ans:=1;
        for i:=1 to b[n,k] do
                ans:=ans*(f[n,k,i]+1) mod 1000000009;
        writeln(ans);
end.