考新郎
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个新郎找错了新娘,求发生这种情况一共有多少种可能.
首先,给每位新娘打扮得几乎一模一样,并盖上大大的红盖头随机坐成一排;
然后,让各位新郎寻找自己的新娘.每人只准找一个,并且不允许多人找一个.
最后,揭开盖头,如果找错了对象就要当众跪搓衣板...
看来做新郎也不是容易的事情...
假设一共有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;
}