bzoj1856

Posted 123456

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj1856相关的知识,希望对你有一定的参考价值。

卡特兰数 组合数优化

毫无头绪。。。

其实我们发现,我们可以把字符串的选择转换一下,我们建立坐标系,起初我们在原点(0,0),每次可以走(1,1)或(1,-1),希望最终到达(n+m,n-m),并且不经过y=-1

那么我们就可以用类似卡特兰数的方法解决这个问题,先算出随便走的方案数,我们一共要走n+m步,然后其中n步是(1,1),其他m步是(1,-1),那么一共的方案数就是

C(n+m,n),然后计算不合法的方案数,也就是经过y=-1的方案数,那么我们转换一下,把起点放在(0,-2),这样其实是对称了一下,把起点关于y=-1作对称点,那么我们现在解一下x,y,x+y=n+m,x-y=n-m+2,所以x=n+1,y=m-1,那么方案数就是C(n+m,n+1)或C(n+m,m-1),两者相减就是答案

这个转化很奥妙啊

技术分享
#include<bits/stdc++.h>
using namespace std;
const long long mod = 20100403ll;
long long n, m;
long long power(long long x, long long t)
{
    long long ret = 1ll;
    for(; t; t >>= 1ll, x = x * x % mod) if(t & 1ll) ret = ret * x % mod;
    return ret;
}
long long fac(long long x) { long long ret = 1ll; for(; x; --x) ret = ret * x % mod; return ret; }
long long inv(long long x) { return power(x, mod - 2ll); }
long long C(long long n, long long m) { return fac(n) * inv(fac(m)) % mod * inv(fac(n - m)) % mod; }
int main()
{
    scanf("%lld%lld", &n, &m);
    printf("%lld\n", ((C(n + m, n) - C(n + m, m - 1)) % mod + mod) % mod);
    return 0;
}
View Code

 

以上是关于bzoj1856的主要内容,如果未能解决你的问题,请参考以下文章

bzoj 1856: [Scoi2010]字符串

bzoj 1856: [Scoi2010]字符串

bzoj1856

bzoj1856 [Scoi2010]字符串

bzoj 1856

字符串(bzoj 1856)