[codeforces 1312D] Count the Arrays 排列组合+快速幂+乘法逆元+二项式定理

Educational Codeforces Round 83 (Rated for Div. 2)   比赛人数8705

[codeforces 1312D]  Count the Arrays   排列组合+快速幂+乘法逆元+二项式定理

总目录详见https://blog.csdn.net/mrcrack/article/details/103564004

也在线测评地址https://codeforces.ml/contest/1312/problem/D

Problem	Lang	Verdict	Time	Memory
D - Count the Arrays	GNU C++11	Accepted	31 ms	0 KB

比赛结束前30秒，程序编好，样例还差最后一组没有测试，顾不了那么多了，马上提交，数据一路AC，突然遇到Time limit exceeded on test 6，此时比赛也已结束，心想，可能是比赛结束，没有测评完毕。该网站，比赛结束，会很快做数据回滚工作，也即，将所有数据重新测试一遍，那就等睡醒后再来看。

再看，奇迹没有发生，还是Time limit exceeded on test 6。重新读题，造了几个数据(数据如下)，发现问题，很快修改，提交AC.总共不到2分钟，要是比赛能有5分钟富余，该题还是能拿下的。

Input:
2 2
Output:
0

Input:
2 3
Output:
0

不过比赛中能AC与比赛一结束就能AC,能力还是有差距的。很高兴，现在到了比赛一结束就能AC，比起以往的无头苍蝇，水平还是高了不少。在比赛结束前,多经历AC与不AC，能极大的刺激参赛者的比赛反应，经过这样的试炼，参赛者水平会高起来的。

该题的正确做法如下，很高兴能在比赛中想出，回头看看，在接触该题时，没有想到思维能这么深入，想出前，与想出后，思维有着天壤之别。比赛时，在看到该题时，觉得自己挺渣的，之后，比赛中，想出公式，样例全通过，真不敢相信，这是自己独立做到的，这种感觉太美妙了。

公式C(m,n-1)*(n-2)*2^(n-3)比较感慨，这是比赛中独立想出的公式

n=5 m=7

m序列
1 2 3 4 5 6 7

讨论，n序列中最大值为4的情况
[1 4 2 3 1]
[1 2 4 3 1]
[1 2 3 4 1]

[2 4 3 2 1]
[1 2 4 3 2]
[1 2 3 4 2]

[3 4 3 2 1]
[1 3 4 3 2]
[1 2 3 4 2]

要组成序列n，需在m中挑n-1个数，可能情况为C(m,n-1)

在挑出的n-1个数中，除了最大的数外，其它的数都是可以充当 一对相等的数,可能情况为(n-1-1)=(n-2),即n-1个数里扣除最大的那个数。

考虑剩下的n-3个数，即n-1个数中，扣除最大的数，扣除一对相等的数，也即n-1-1-1=n-3:
这n-3个数，在最大数左边可能有0个,1个,2个,...,n-3个
0个对应C(n-3,0)
1个对应C(n-3,1)
2个对应C(n-3,2)
...
n-3个对应C(n-3,n-3)

这n-3个数的情况加起来
C(n-3,0)+C(n-3,1)+C(n-3,2)+...+C(n-3,n-3)=2^(n-3)对应二项式定理

综合上述情况，对应公式
C(m,n-1)*(n-2)*2^(n-3)

AC代码如下，别忘了n==2的特判

#include <stdio.h>
#define mod 998244353
#define LL long long
LL quick_pow(LL a,LL b){//快速幂
	LL ans=1;
	while(b){
		if(b&1)ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}
LL C(LL a,LL b){//计算组合数
	LL i,ans=1,d=1;
	for(i=a;i>=a-b+1;i--)ans=ans*i%mod;//分子
	for(i=b;i>=1;i--)d=d*i%mod;//分母
	ans=ans*quick_pow(d,mod-2)%mod;
	return ans;
}
int main(){
	LL n,m,ans;
	scanf("%lld%lld",&n,&m);
	if(n==2){//特判,提交时太匆忙，没考虑到，造成了Time limit exceeded on test 6
		printf("0\n");
		return 0;
	}
	ans=C(m,n-1)*(n-2)%mod*quick_pow(2,n-3)%mod;
	printf("%lld\n",ans);
	return 0;
}