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
Sample Output
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;
}