BNUOJ 34963 相似数列

Time Limit: 3000ms

Case Time Limit: 3000ms

Memory Limit: 65536KB

64-bit integer IO format:  %lld      Java class name:  Main

众所周知Liserious的一大爱好是在草稿纸上写随机数列，有一天他突然发现5,2,3,3,3和3,5,2,2,2很相似的样子。但是到底是哪里相似呢，李希尔瑞斯还在思考。

2.333年过去了。

Liserious在学习字母表的时候突然明白了，对于一个数列，他把第一个数字标成A，对于后面的数字，如果前面出现过，就标注成前面标的字母，否则标成未使用过的第一个字母（字典序）。这样如果两个序列标注的结果是一样的，那它们就是相似的。

李希尔瑞斯被自己惊呆了！

现在他想知道，对于N个数字的数列，总共有多少个互不相似的数列呢？同时，Liserious不希望得到的数列太过简单，所以他不会计算不同数字少于K个的数列。

 

Input

每行两个正整数N（0<N<=1000），K（0<K<=N）表示数列的长度，以及最少不同数字的个数。读入到EOF为止。

 

Output

对于每组数据，输出一个整数表示总共有多少互不相似的数列，模去1000000007。

 

Sample Input

1 1
2 1

Sample Output

1
2

一个矩阵问题，需要一定的线代功底。

推出公式就非常简单了。

当然题目描述有一点点小问题：用来替换的字母数是无限的。也就是ABC...XYZ后面还有好多好多不知名的字母。

#include <stdio.h>
#include <fstream>
#include <string.h>
#include <iostream>
#include <math.h>
#include <algorithm>
#include <vector>
#include <map>
#define PI 3.1415926
#define M 2005  //10^5
#define eps 1e-8
#define moo 1000000007
using namespace std;
long long a[1005][1005];

int main ()
{
    for(int i=1;i<=1000;i++)
        for(int j=1;j<=1000;j++)
    {
        a[i][j]=0;
    }  //之前这里忘了重置，WA了4次才发现
    for(int i=1;i<=1000;i++)
        a[i][1]=1;
    for(int i=2;i<=1000;i++)
        for(int j=2;j<=i;j++)
    {
        a[i][j]=(a[i-1][j-1]+a[i-1][j]*j)%moo;
    }
//以上为公式部分，注意因为是按字典序取字母替换，所以不需要乘组合数。
    int n,k;
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        long long ans=0;
        for(int i=k;i<=n;i++)
        {
            //if(i<=26)
                ans=(ans+a[n][i])%moo;
        }
        printf("%lld\n",ans);
    }
}