洛谷 P3904 三只小猪

题目背景

你听说过三只小猪的故事吗？这是一个经典的故事。很久很久以前，有三只小猪。第一只小猪用稻草建的房子，第二个小猪用木棍建的房子，第三个小猪则使用砖做为材料。一只大灰狼想吃掉它们并吹倒了稻草和木棍建的房子。但是砖盖的房子很结实，狼最终也没有破坏掉，最后小猪们战胜了大灰狼并把它尾巴烧掉了。

题目描述

为了自己的安全，小猪们又见了一个新砖房。但是现在问题出现了，怎样把三个小猪分配到两个房子里呢？第三只小猪是三只小猪中最聪明的一只，为了不浪费任何一个房子，它总共考虑了三种方案，如下图

“但是将来怎么办呢？”第三只小猪知道将来随着成员的增多，它们将会盖更多的房子。它想知道给定了房子和猪的数目后，房子的分配方案有多少，但这个问题对于它来说，很明显有点难了，你能帮小猪解决这个问题吗？

输入输出格式

输入格式：

输入文件piggy.in，仅有一行，包含两个整数n和m，分别表示小猪的数目和房间数（1≤n≤50，0≤m≤50）。

输出格式：

输出文件piggy.out，仅一个整数，表示将n只小猪安置在m个房间且没有房间空闲的方案数。

输入输出样例

输入样例#1：

4 2

输出样例#1：

7

输入样例#2：

6 7

输出样例#2：

0

题解：高精+ Stirling数
粘过上篇题解代码就A了....重题...
递推公式s[i][j]=s[i-1][j]*j+s[i-1][j-1]
代码：

#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 ;
}