【2020牛客多校一：数学积分】J ：Easy Integration

传送门

【难度】
4/10
数学题/或找规律题

【题意】

求出$I_n$
其中$I_n={\int }_0^1(x-x^2{)}^{n}dx\text{In=∫01(x−x2)ndx}{I}_n={\int }_0^1(x-x^2{)}^{n}dx$

【数据范围】

T组样例
$T\le 1e5$
$n\le 1e6$

【输入样例】

1
2
3

【输出样例】

166374059
432572553
591816295

【解释】

$I_1={\int }_0^1(x-x^2)dx=\frac{1}{6}\text{I1=∫01(x−x2)dx=61}{I}_1={\int }_0^1(x-x^2)dx=\frac{1}{6}$

【思路】

据说是wallis integral ，但是看了公式页许久也不知道为什么

我们利用分部积分做。
即：
$${\int }_a^{b}u{v}^{\prime }dx=uv{|}_a^b{\int }_a^b{u}^{\prime }vdx$$
接下来先化简公式再进行运算：
$$\begin{gathered} \begin{array}{rl}{I}_n& ={\int }_0^1(x-x^2{)}^{n}dx\\ & ={\int }_0^1{x}^n(1-x{)}^{n}dx\\ & =\frac{1}{n+1}{x}^{n+1}(1-x{)}^n{|}_0^1+{\int }_0^{1}n(1-x{)}^{n-1}\frac{1}{n+1}{x}^{n+1}dx\\ & ={\int }_0^{1}n(1-x{)}^{n-1}\frac{1}{n+1}{x}^{n+1}dx（左边项为0）\end{array} \\ ⋮(n次分部积分后) \\ \begin{array}{rl}{I}_n& =\frac{n\ast (n-1)\ast (n-2)\ast \cdots \ast 2\ast 1}{(n+1)\ast (n+2)\ast \cdots \ast (2n)}{\int }_0^1(1-x{)}^0{x}^{2n}dx\\ & =\frac{n\ast (n-1)\ast (n-2)\ast \cdots \ast 2\ast 1}{(n+1)\ast (n+2)\ast \cdots \ast (2n)}\frac{1}{2n+1}{x}^{2n+1}{|}_0^1\\ & =\frac{n\ast (n-1)\ast (n-2)\ast \cdots \ast 2\ast 1}{(n+1)\ast (n+2)\ast \cdots \ast (2n+1)}\\ & =\frac{(n!{)}^2}{(2n+1)!}\end{array} \end{gathered}$$

接下来我们只要先预处理出阶乘数组与阶乘逆元数组

【核心代码】

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
const int MAX = 2e6+50;
const ll  MOD = 998244353;

ll qpow(ll a,ll n){/* */ll res = 1LL;while(n){if(n&1)res=res*a%MOD;a=a*a%MOD;n>>=1;}return res;}
ll qpow(ll a,ll n,ll p){ll res = 1LL;while(n){if(n&1)res=res*a%p;a=a*a%p;n>>=1;/**/}return res;}
ll inv(ll a){/* */return qpow(a,MOD-2);}
ll inv(ll a,ll p){return qpow(a,p-2,p);}

ll fac[MAX];
ll inv_fac[MAX];
void init(int n){
    fac[0]=1LL;
    for(int i=1;i<=n;++i){
        fac[i]=(fac[i-1]*i)%MOD;
    }
    inv_fac[n]=inv(fac[n]);
    for(int i=n-1;i>=1;--i){
        inv_fac[i]=(inv_fac[i+1]*(i+1))%MOD;
    }
}
int main()
{
    init(MAX-50);
    ll n;
    while(cin >> n){
        cout << (fac[n]*fac[n])%MOD*inv_fac[2*n+1]%MOD << endl;
    }
    return 0;
}