2023牛客五一集训派对Day1 J-Easy Integration

没没没没没超疯还是疯了

已于 2023-05-03 21:03:35 修改

阅读量83

点赞数 1

分类专栏：算法学习文章标签：算法 c++ 开发语言

于 2023-05-03 20:51:18 首次发布

本文链接：https://blog.csdn.net/m0_63033104/article/details/130476982

版权

算法学习专栏收录该内容

14 篇文章 0 订阅

订阅专栏

2023牛客五一集训派对Day1 J-Easy Integration

链接：https://ac.nowcoder.com/acm/contest/55992/J

头一回做一个需要求积分的题目。爽的嘞~~蒟蒻小姐积分不会求，求出积分也不会写代码 😢

题目描述

Given $n$ , find the value of $\int_{0}^1 (x - x^2)^n \mathrm{d} x$ .

It can be proved that the value is a rational number $\frac{p}{q}$ .

Print the result as $\cdot q^{-1}) \bmod 998244353$ .

输入描述

The input consists of several test cases and is terminated by end-of-file.

Each test case contains an integer $n$ .

* $\leq n \leq 10^6$
* The number of test cases does not exceed $10^5$ .

输出描述

For each test case, print an integer which denotes the result.

示例

输入：
1
2
3
输出：
166374059
432572553
591816295

备注

For $n = 1$ , $\int_{0}^1 (x - x^2) \mathrm{d} x = \frac{x^2}{2} - \frac{x^3}{3} |_0^1 = \frac{1}{6}$

题解

先计算公式：

在这里插入图片描述

得出结论： $\int_{0}^1 (x - x^2)^n \mathrm{d} x=\frac{n!·n!}{(2n-1)!}$

此时题意转化为：给定整数 $n$ ，求 $\frac{n!·n!}{(2n-1)!} \bmod 998244353 $ 的值。

注意⚠：

求阶乘用数组 $f$ 记录，类似于前缀和的做法。
使用快速幂求逆元将原本的除法转换为乘法（除法计算会有精度缺失等问题）。

代码

#include <iostream>
#include <cmath>
using namespace std;
const int N=1e6+7, MOD=998244353;
long long f[2*N]; //f[i]表示i的阶乘
long long fun(long long base, long long power) //快速幂模板
{
    long long result=1;
    while(power>0)
    {
        if(power&1) result=result*base%MOD;
        power>>=1;
        base=(base*base)%MOD;
    }
    return result;
}
int main()
{
    int n;
    long long ans=1;
    for(long long i=1; i<=2*N+1; ++i)
    {
        ans*=i;
        ans%=MOD;
        f[i]=ans;
    }
    while(cin>>n)
    {
        long long a, b;
        a=f[n]%MOD, b=f[2*n+1]%MOD;
        cout<<(a*a)%MOD*fun(b, MOD-2)%MOD<<endl; //快速幂求逆元将原本的除法转换为乘法
    }
    return 0;
}