zoukankan      html  css  js  c++  java
  • 卡特兰数 + 大数

    The "Harry Potter and the Goblet of Fire" will be on show in the next few days. As a crazy fan of Harry Potter, you will go to the cinema and have the first sight, won’t you?

    Suppose the cinema only has one ticket-office and the price for per-ticket is 50 dollars. The queue for buying the tickets is consisted of m + n persons (m persons each only has the 50-dollar bill and n persons each only has the 100-dollar bill).

    Now the problem for you is to calculate the number of different ways of the queue that the buying process won't be stopped from the first person till the last person.
    Note: initially the ticket-office has no money.

    The buying process will be stopped on the occasion that the ticket-office has no 50-dollar bill but the first person of the queue only has the 100-dollar bill.

    InputThe input file contains several test cases. Each test case is made up of two integer numbers: m and n. It is terminated by m = n = 0. Otherwise, m, n <=100.
    OutputFor each test case, first print the test number (counting from 1) in one line, then output the number of different ways in another line.
    Sample Input
    3 0
    3 1
    3 3
    0 0
    Sample Output
    Test #1:
    6
    Test #2:
    18
    Test #3:
    180

    题意 : 经典的卡特兰数问题,最终的答案是 C(m+n, m) - C(m+n, m+1) , 再进一步化简可得 (m+n)!(m+1-n)/(m+1) , 一个大数乘法就OK
    代码示例 :
    #define ll long long
    const int maxn = 1e6+5;
    const int mod = 1e9+7;
    const double eps = 1e-9;
    const double pi = acos(-1.0);
    const int inf = 0x3f3f3f3f;
    
    int m, n;
    vector<int>a, b, c;
    stack<int>s;
    
    void bigsum(){
        c.clear();
        c.assign(a.size()+b.size()-1, 0);
        for(int i = 0; i < b.size(); i++){
            int k = i;
            for(int j = 0; j < a.size(); j++){
                c[k] += b[i]*a[j];
                k++;
            }
        }
        
        for(int i = c.size()-1; i >= 0; i--){
            if (c[i] > 9){
                if (i != 0){c[i-1] += c[i]/10; c[i] %= 10;}
                else {
                    int tem = c[i]/10;
                    c[i] %= 10;
                    c.insert(c.begin(), tem);
                }
            }
        }
    }
    
    int main() {
        //freopen("in.txt", "r", stdin);
        //freopen("out.txt", "w", stdout);
        int kase = 1;
        
        while(~scanf("%d%d", &m, &n) && n+m){
            printf("Test #%d:
    ", kase++);
            if (m < n) printf("0
    ");
            else{
                b.clear();
                b.push_back(1);
                for(int i = 2; i <= m+n; i++){
                    if (i == m+1) continue;
                    a.clear();
                    int x = i;
                    while(x) s.push(x%10), x /= 10;
                    while(!s.empty()){
                        a.push_back(s.top());
                        s.pop();
                    } 
                    bigsum();
                    b = c;       
                }
                
                a.clear();
                int x = m+1-n;
                if (x != m+1) {
                    while(x) s.push(x%10), x /= 10;
                    while(!s.empty()){
                        a.push_back(s.top());
                        s.pop();
                    } 
                    bigsum();
                }
                for(int i = 0; i < c.size(); i++){
                    printf("%d", c[i]);
                }
                printf("
    ");
            }
        }
        return 0;
    }
    


    东北日出西边雨 道是无情却有情
  • 相关阅读:
    Oracle EBS-SQL (MRP-7):检查MRP计划运行报错原因之超大数据查询2.sql
    Oracle EBS-SQL (MRP-6):检查MRP计划运行报错原因之超大数据查询1.sql
    Oracle EBS-SQL (MRP-5):重起MRP Manager.sql
    Oracle EBS-SQL (INV-12):检查待定事物处理1.sql
    Oracle EBS-SQL (OM-6):打开订单.sql
    Oracle EBS-SQL (INV-11):检查子库存会计信息.sql
    Oracle EBS-SQL (OM-5):检查订单事务处理类型.sql
    Oracle EBS-SQL (OM-4):检查发运网络.sql
    Oracle EBS-SQL (PO-18):检查工作台下达的PR在系统找不到.sql
    Oracle EBS-SQL (PO-17):检查供货比例不为100%.sql
  • 原文地址:https://www.cnblogs.com/ccut-ry/p/8944889.html
Copyright © 2011-2022 走看看