洛谷 P1655 小朋友的球

题目描述

@发源于 小朋友最近特别喜欢球。有一天他脑子抽了，从口袋里拿出了N个不同的球，想把它们放到M个相同的盒子里，并且要求每个盒子中至少要有一个球，他好奇有几种放法，于是尝试编程实现，但由于他天天不好好学习，只会上B站看游泳教练，于是他向你求助。

输入输出格式

输入格式：

多组数据，每行两个数N，M。

输出格式：

每组数据一行，表示方案数。

输入输出样例

输入样例#1：

4 2

1 1

输出样例#1：

7

1

说明

【样例解释】

N=4,M=2

1,2 3 4

2,1 3 4

3,1 2 4

4,1 2 3

1 2,3 4

1 3,2 4

1 4,2 3

对于20%的数据，满足1≤N，M≤10；

对于100%的数据，满足1≤N，M≤100，数据组数≤10。

题解：高精+第二类Stirling数

递推公式s[i][j]=s[i-1][j]*j+s[i-1][j-1]

s[i][j]的一个组合学解释是：将i个物体划分成j个非空的不可辨别的（可以理解为盒子没有编号）集合的方法数。

考虑最后一个物品的去向，要么单独一个盒子，要么随便选一个盒子。

初始化：s[i][1]=s[i][i]=1,m<n||n==0 return 0;

ps:Candy?的高精模板太好看了。%

代码：

#include<iostream>

#include<cstdio>

#include<cstring>

#define LL long long

#define B 10

using namespace std;

LL m,n;

struct Big{

    int a[], n;

    int& operator [](int x) {return a[x];}

    Big():n() {memset(a, , sizeof(a));}

    void ini(int x) {a[]=x; n=;}

}t[][];

Big operator *(Big a, int b) {

    int g=;

    for(int i=; i<=a.n; i++)

        g += a[i]*b, a[i] = g%, g/=;

    if(g) a[++a.n] = g;

    return a;

}

Big operator *(Big a, Big b) {

    Big c;

    for(int i=; i<=a.n; i++) {

        int g=;

        for(int j=; j<=b.n; j++)

            g += c[i+j-]+a[i]*b[j], c[i+j-] = g%, g/=;

        c[i+b.n] = g;

    }

    c.n = a.n + b.n;

    while(c.n> && c[c.n]==) c.n--;

    return c;

}

Big operator +(Big a, Big b) {

    int g=, n=max(a.n, b.n);

    for(int i=; i<=n; i++) {

        g += i<=a.n ? a[i] : ;

        g += i<=b.n ? b[i] : ;

        a[i] = g%, g/=;

    }

    a.n = n;

    if(g) a[++a.n] = g;

    return a;

}

Big operator -(Big a, Big b) {

    for(int i=; i<=b.n; i++) {

        if(a[i]<b[i]) a[i]+=, a[i+]--;

        a[i] -= b[i];

    }

    int p=b.n+;

    while(a[p]<) a[p]+=, a[++p]--;

    while(a.n> && a[a.n]==) a.n--;

    return a;

}

void Print(Big &a) {

    printf("%d", a[a.n]);

    for(int i=a.n-; i>=; i--) printf("%d", a[i]);

}

int main(){

    for(int i=;i<=;i++){

        t[i][].a[]=;t[i][i].a[]=t[i][].a[]=;

    }

    for(int i=;i<=;i++)

     for(int j=;j<=i-;j++)

      t[i][j]=t[i-][j-]+t[i-][j]*j;

    while(scanf("%lld%lld",&m,&n)!=EOF){

        if(n==||n>m){printf("0\n");continue;}

        Print(t[m][n]);

        printf("\n");

    }

    return ;

}