CF1204E Natasha, Sasha and the Prefix Sums

C202044zxy

于 2020-09-09 09:47:44 发布

阅读量216

点赞数 1

分类专栏：组合计数

本文链接：https://blog.csdn.net/C202044zxy/article/details/108461395

版权

组合计数专栏收录该内容

15 篇文章 0 订阅

订阅专栏

一、题目

点此看题

二、解法

看到这个前缀和比较容易转化成二维平面上移动的模型，每次向右走一格， $+ 1$ 多向上走一格， $- 1$ 多向下走一格。我们可以枚举这个最大值，那么答案就是（ $f_i$ 表示最大值为 $i$ 的方案数）：
$\sum_{i=1}^n if_i=\sum_{i=1}^n g_i$ 其中 $g_i$ 表示最大值至少为 $i$ 的方案数，不难发现他们等价，至少的方案数会好求一些。

起终点我们是知道的，要求最大值至少为 $i$ 意味着路径上要经过 $y = i$ 这条直线，分为两种情况讨论：起终点在直线异侧，起终点在直线同侧。第一种情况是很好算的，因为一定会经过直线，就是 $C (n + m, n)$ ，第二种情况我们还是结合图形分析：
在这里插入图片描述
我们把 $(0, 0)$ 翻折上去，那么从 $(0, 2 i)$ 到终点一定对应了 $(0, 0)$ 到终点且经过 $y = i$ 的所有方案，通过解方程可知它上移了 $n - i$ 步，下移了 $m + i$ 步，那么这样的方案就是 $C (n + m, n - i)$

#include <cstdio>
const int M = 4005;
const int MOD = 998244853;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,ans,fac[M],inv[M];
void init(int n)
{
	inv[0]=inv[1]=fac[0]=fac[1]=1;
	for(int i=2;i<=n;i++) inv[i]=1ll*(MOD-MOD/i)*inv[MOD%i]%MOD;
	for(int i=2;i<=n;i++) inv[i]=1ll*inv[i]*inv[i-1]%MOD;
	for(int i=2;i<=n;i++) fac[i]=1ll*fac[i-1]*i%MOD;
}
int C(int n,int m)
{
	return 1ll*fac[n]*inv[m]%MOD*inv[n-m]%MOD;
}
signed main()
{
	init(4000);
	n=read();m=read();
	for(int i=1;i<=n;i++)
	{
		if(i<=n-m) ans=(ans+C(n+m,n))%MOD;
		else ans=(ans+C(n+m,n-i))%MOD;
	}
	printf("%d\n",ans);
}