2020牛客多校第一场J题[Easy Integration(逆元+数学推导打表)

牛客多校第一场J题 Easy Integration

先解决数学问题解积分表达式，听说这是高数题目，（水平还需要提高，这种式子算起来还是不算太难的,但是当时死活不知道咋求，放假放傻了…高中爷orz
$$\begin{array}{rl}& {\int }_0^1{x}^n(1-x{)}^{n}dx\\ =& \frac{1}{n+1}{\int }_0^1(1-x{)}^{n}d{x}^{n+1}\\ =& \frac{1}{n+1}{x}^{n+1}(1-x{)}^n-\frac{n}{n+1}{\int }_0^1{x}^{n+1}(1-x{)}^{n-1}(-1)dx\\ \because & 前一项积分上下限01均等于0所以可省略，经过n次约简得到\\ =& \frac{n!}{(n+1)(n+2)\times (2n-1)(2n)}{\int }_0^1{x}^{2n}dx\\ =& \frac{n{!}^2}{(2n+1)!}\end{array}$$

最后这个阶乘式子表达为最终的有理式$\frac{p}{q}$形式

询问1e5,n最大1e6

所以打阶乘表解决问题，用long long直接乘法+%p

p是素数，由费马小定理和逆元知识得到：$a{b}^{-1}\equiv a{b}^{p-2}(modp)$.

逆元，快速幂，数学，打表

最后AC，补题成功

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const ll mod = 998244353;
ll F[2000020],F_1[2000020];
ll n,up,down;
ll pow(ll a,ll b){
	ll ans = 1;
	for(;b;b >>= 1){
		if(b & 1){
			ans = (ans*a)%mod;
		}
		a = a*a%mod;
	}
	return ans;
}


int main(){
	//预处理阶乘数组
	F[0]=1;
	for(int i=1;i<=2000005;i++){
		F[i] = F[i-1] * i % mod; 
		// F_1[i] = pow(F[i],mod-2);
	}

	while(scanf("%lld",&n)!=EOF){
		up = F[n]*F[n] % mod;
		down = pow(F[2*n+1],mod-2)%mod;
		printf("%lld\n",up*down%mod);
	}
	return 0;
}