HDU 2049 错位排序 + 高精度

考新郎

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 14155    Accepted Submission(s): 5281


Problem Description
国庆期间,省城HZ刚刚举行了一场盛大的集体婚礼,为了使婚礼进行的丰富一些,司仪临时想出了有一个有意思的节目,叫做"考新郎",具体的操作是这样的:


首先,给每位新娘打扮得几乎一模一样,并盖上大大的红盖头随机坐成一排;
然后,让各位新郎寻找自己的新娘.每人只准找一个,并且不允许多人找一个.
最后,揭开盖头,如果找错了对象就要当众跪搓衣板...

看来做新郎也不是容易的事情...

假设一共有N对新婚夫妇,其中有M个新郎找错了新娘,求发生这种情况一共有多少种可能.
 

Input
输入数据的第一行是一个整数C,表示测试实例的个数,然后是C行数据,每行包含两个整数N和M(1<M<=N<=20)。
 

Output
对于每个测试实例,请输出一共有多少种发生这种情况的可能,每个实例的输出占一行。
 

Sample Input
  
  
2 2 2 3 2
 

Sample Output
  
  
1 3
 
全错位排列的解法在百度百科中有: http://baike.baidu.com/view/1926671.htm

#include<iostream>
#include<vector>
using namespace std;

struct bignum{
    vector<int> bit;
    void multi(int k);
    void div(int k);
    void inc();
    void dec();
    void display()
    {
        int idx = bit.size();
        while(idx--)  cout<<bit[idx];
        cout<<endl;
    }
}B[21];

void bignum::inc()
{
    vector<int>::iterator pos=bit.begin();
    for(pos=bit.begin();pos!=bit.end();pos++){
        if(*pos != 9){
            (*pos)++;
            return;
        }else{
            *pos = 0;
        }
    }
    bit.push_back(1);
}
void bignum::dec()
{
    vector<int>::iterator pos;
    for(pos=bit.begin();pos!=bit.end();pos++){
        if(*pos != 0){
            (*pos)--;
            return ;
        }else{
            *pos = 9;
        }
    }
}
void bignum::multi(int k)
{
    int more = 0;
    vector<int>::iterator pos;
    for(pos=bit.begin();pos!=bit.end();pos++)
    {
        *pos = (*pos) * k + more;
        more = *pos / 10;
        *pos %= 10;
    }
    while(more){
        bit.push_back(more%10);
        more /= 10;
    }
}
void bignum::div(int k)
{
    int more = 0;
    vector<int>::reverse_iterator p;
    for(p=bit.rbegin();p!=bit.rend();p++)
    {
        *p = (*p+more*10);
        more = *p % k;
        *p /= k;
    }
    for(p = bit.rbegin();*p==0;p++) bit.pop_back();
}

void init()
{
    B[1].bit.push_back(0);    
    for(int i=2;i<=20;i++){
        B[i] = B[i-1];
        B[i].multi(i);
        if(i&1)  B[i].dec();
        else     B[i].inc();
    }
}

void Solve(int n,int m)
{
    bignum ans = B[m];
    for(int i=1;i<=m;i++,n--){
        ans.multi(n);
        ans.div(i);
    }
    ans.display();
}

int main()
{
    init();
    int t;  cin>>t;
    while(t--){
        int m,n;   cin>>n>>m;
        Solve(n,m);
    }    
//    system("pause");
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值