codeforces Helvetic Coding Contest 2016 D2 The Wall (medium)

本文解析了CodeForces平台上的D2题,通过枚举砖块使用数量的方法,结合组合数学中的卢卡斯定理,给出了求解不同建筑方案数目的算法实现。此外还提到了一个O(1)复杂度的解法。

####题目链接:
http://codeforces.com/contest/690/problem/D2
####题目大意:
一共有n块砖,筑成一道宽度为不超过c的墙,有多少的建筑方案?(n块不必要都用)
####结题思路:
枚举用的砖数,用的砖数为i时的建筑方案数为 ( i i + c − 1 ) \left(^{i+c-1}_{i}\right) (ii+c1)
这有一个定理把n个相同的物品放入m个不同的盒子中(盒子可为空)的情况数 ( m − 1 n + m − 1 ) \left(^{n+m-1}_{m-1}\right) (m1n+m1)。这个公式可这样理解有 n + m − 1 n+m-1 n+m1个物品,把其中 m − 1 m-1 m1个变成分割线,原来的n个物品被分成m部分。
所以原题就变成 ∑ 1 n ( i i + c − 1 ) \sum_1^n\left(^{i+c-1}_{i}\right) 1n(ii+c1),求组合数用卢卡斯定理,这题就解决了。
另外还有o(1)的解法, a n s = ( m i n ( n , c ) m a x ( n , c ) ) − 1 ans=\left(^{max(n,c)}_{min(n,c)}\right)-1 ans=(min(n,c)max(n,c))1,至于为什么还不知道。

#include<cstdio>
#include<cstdlib>
#define LL long long
const int mod=1000003;
LL f[700010];  //N为组合数的底数 的范围    
void init(int p){  
    f[0] = 1;  
    for(int i = 1; i <= 700000; ++i)  
        f[i]  = f[i-1] * i % p;  
}  
LL pow_mod(LL a, LL x, int p){  
    LL ret = 1;  
    a %= p;  
    while(x){  
        if(x & 1){  
            ret = ret * a % p;  
            --x;  
        }  
        else{  
            a = a * a % p;  
            x >>= 1;    
        }  
    }  
    return ret;  
}  
LL Lucas(LL n, LL k, int p){  
    LL ret = 1;  
    while(n && k){  
        LL nn = n % p, kk = k % p;  
        if(nn < kk) return 0;  
        ret = ret * f[nn] * pow_mod(f[kk] * f[nn - kk] % p, p - 2, p) % p;  
        n /= p;  
        k /= p;  
    }  
    return ret;  
}  
int main(){
	int n,c;
	LL ans=0;
	init(mod);
	scanf("%d%d",&n,&c);
	for(int i=1;i<=n;++i){
		ans=(ans+Lucas(i+c-1,i,mod))%mod;
	}
	printf("%I64d\n",ans);
	//system("pause");
	return 0;
}
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值