HDOJ HDU 1133 Buy the Ticket

HDOJ 1133 Buy the Ticket

题目

点此查看 HDOJ 1133 Buy the Ticket

分类

找规律

题意

电影院售票
50 元 一张
m 人 有 50元 n 人 有100元
这些人排队买票
收营员要拿收到的 50 找 100 没有 50 队列就停止
求符合条件序列的个数

题解

m < n 时 不存在合法方案
当 m >= n 时
合法排列有定有 $m! \times n! \times C_{m+n}^m$ 种 排法
而有 $m! \times n! \times C_{m+n}^{m+1}$ 非法排法 (从 有 1 个50前边至少多1个100)
所以 有 合法排法 $m! \times n! \times ( C_{m+n}^m - C_{m+n}^{m+1} )$ = $(m+n)! \times {(m - n + 1) \over (m+1)}$

代码

#include <iostream>
#include <iomanip>
#include <cstring>

#define maxnum 500
#define maxint 10000
using namespace std;

struct bigint
{
    int m;
    int perint[maxnum];

    bigint()
    {
        for(int i = 0;i < maxnum;i++)
            perint[i] = 0;
    }
    bigint & operator*(int b)
    {
        int c = 0,cur;
        for(int i = 0;i <= m;i++)
        {
            cur = perint[i] * b + c;
            perint[i] = cur % maxint;
            c = cur / maxint;
        }
        if(c)
        {
            m++;
            perint[m] = c;
        }
        return *this;
    }
    bigint & operator/(int b)
    {
        int c = 0,cur;
        for(int i = m;i >=0 ;i--)
        {
            cur = c * maxint + perint[i];
            perint[i] = cur / b;
            c = cur % b;
        }
        return *this;
    }
    friend ostream & operator<<(ostream &out,bigint & b)
    {
        int s;
        for(s = b.m;s >= 0;s--)
            if(b.perint[s])
                break;
        out << b.perint[s];
        for(int i = s - 1;i >= 0;i--)
                out << setw(4) << setfill('0') << b.perint[i];
        return out;
    }
};

int main()
{
    int f,h,z = 1;
    bigint t;
    while(cin >> f >> h && (f || h))
    {
        if(h > f)
        {
            cout << "Test #" << z++ << ':' << endl; 
            cout << 0 << endl;
            continue;
        }
        t.perint[0] = 1;
        t.m = 0;
        for(int i = 2;i <= (f + h);i++)
            t * i;
        t * (f + 1 - h);
        t / (f + 1);
        cout << "Test #" << z++ << ':' << endl; 
        cout << t << endl;
        memset(t.perint,0,sizeof(t.perint));
    }
    return 0;
}