杭电 2200 Eddy's AC难题（排列组合）

Eddy's AC难题

Time Limit : 3000/1000ms (Java/Other)   Memory Limit : 32768/32768K (Java/Other)

Total Submission(s) : 10   Accepted Submission(s) : 3

Font: Times New Roman | Verdana | Georgia

Font Size: ← →

Problem Description

Eddy是个ACMer,他不仅喜欢做ACM题,而且对于Ranklist中每个人的ac数量也有一定的研究,他在无聊时经常在纸上把Ranklist上每个人的ac题目的数量摘录下来，然后从中选择一部分人(或者全部)按照ac的数量分成两组进行比较，他想使第一组中的最小ac数大于第二组中的最大ac数，但是这样的情况会有很多，聪明的你知道这样的情况有多少种吗?

特别说明：为了问题的简化，我们这里假设摘录下的人数为n人，而且每个人ac的数量不会相等，最后结果在64位整数范围内.

Input

输入包含多组数据，每组包含一个整数n,表示从Ranklist上摘录的总人数。

Output

对于每个实例，输出符合要求的总的方案数，每个输出占一行。

Sample Input

2
4

Sample Output

1
17

Author

Eddy 

/*
插板法，C(n,n)*(n-1)+C(n,n-1)*(n-2)+C(n,n-2)*(n-3)……C(n,2)*1
即为总个数                                           C(n,n/2)有两项
如果n是奇数，刚好是C(n,n/2)*(n/2-1)+C(n,n-n/2)*((n-n/2)-1)刚好是C(n,n/2)
如果n是偶数,刚好C(n,n-2)有两项
所以化简一下
=(n-1)+C(n,1)*(n-2)+C(n,2)*(n-2)+C(n,3)*(n-2)……C(n,n/2)*(n-2);
=(n-1)+(n-2)*2^(n/2);
Time:2015-4-14 17:37
*/
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
int main(){
    LL n,m;
    while(scanf("%lld",&n)!=EOF){
        LL ans=1+(n-2)*(1<<(n-1));
        printf("%lld\n",ans);
    }
return 0;
}

组合数的定义：C(n,m)=n!/( (n-m)!*m! )

计算组合数主要头疼的是溢出，long long 类型的数字算C(82,41)已经不行了。。。

一、普通算法

    由于溢出问题严重，所以算出三个阶乘再做除法的话，中间结果会溢出。

    首先做个小优化，利用 C(n,m) = C(n,n-m) ，如果m超过n的一半就让 m = n-m。

    这样处理之后，m一定是小于等于n-m的

    对于：1，2，3，…，m，m+1，…，n-m，…，n

    公式中 n! , (n-m)!, m!肯定有重叠的部分，所以把n!和(n-m)!中重叠部分消去，直接算(n-m+1)*(n-m+2)*…*n / m!（红色部分的连乘）就好，这样时间复杂度是O( m )，已经和n没有关系了，而且就算要计算 C(10000,3)也可以，因为中间结果最大是1000*999*998。

    还有个小优化，为了防止溢出，一边计算 (n-m+1)*(n-m+2)*…*n ，一边除掉1,2,3,...,m中能除的数（下面的算法是按顺序除的）

#include<stdio.h>
#define maxn 500
__int64 sum[maxn];


__int64 Comb(__int64 n,__int64 m)
    {
    int sum0=1;int k=1;
    for(int i=n-m+1;i<=n;i++)
    {
 sum0=sum0*i;
 while(sum0%k==0&&k<=m)
 {
   sum0/=k;
   k++;
 }
}
return sum0;
}
void dabiao()
{

int x;
for(int i=2;i<maxn;i++)
{
 __int64 num=0;
       for(int j=2;j<=i;j++)
       {
         if(j>i/2 ) x=i-j;因为组合数的性质
else x=j;
 
     num+=Comb(i,x)*(j-1);/*注：这个地方Comb(i,x)是从i个数中选取j个数进行组合，组合之后又j-1个空，进行插板。这儿不能用x,x是根据组合数性质算的组合，实际上                            是有j个元素*/
  }
     sum[i]=num;
}
     
}
int main()
{
__int64 m,i;
  dabiao();
while(scanf("%I64d",&m)!=EOF)
 {
     printf("%I64d\n",sum[m]);
 }
return 0;
}