感受一下dp优化的魅力

求1!+2!+3!+...+n!，由于答案很大，答案对998244353求余即可。

输入格式

输入存在多组测试数据，对于每组测试数据输入一个正整数n(1<=n<=10^6)

输出格式

对于每组测试数据输出一行表示答案。

输入样例

1
5

输出样例

1
153

方法1 输入一个数从1开始计算阶乘

 

# include <iostream>
#include<algorithm>

using namespace std;

const int MOD=998244353;
typedef unsigned long long ull;
int main(void)
{
    ull n ;
    ull i = 0;
    ull m = 1;
    ull sum=0;
    //printf("请输入n的值:");
    while(cin>>n)
    {
        for (i=1; i<=n; ++i)
        {
            m = ((m%MOD) * (i%MOD))%MOD;
            sum = (sum +m)%MOD;
        }
        cout<<sum<<endl;
    }

    return 0;
}

方法2：预处理阶乘

 

#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<stack>
#include<algorithm>
#include<string>
using namespace std;

typedef unsigned long long ull;
//const ull MOD=0x3f3f3f3f3f;
const ull MOD=998244353;
const double PI=3.1415926;

ull dp[1000010];
//string a,b;
int main()
{
    int n;
    dp[1]=1;
    while(cin>>n)
    {
        ull sum=0;
        for(int i=1;i<=n;i++)
        {
            if(dp[i]==0)
            {
                dp[i]=((dp[i-1]%MOD)*(i%MOD))%MOD;
            }
            sum=(sum+dp[i])%MOD;
        }
        cout<<sum<<endl;
    }

    return 0;
}

方法3：预处理阶乘  和   阶乘之和

 

#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<stack>
#include<algorithm>
#include<string>
using namespace std;

typedef unsigned long long ull;
//const ull MOD=0x3f3f3f3f3f;
const ull MOD=998244353;
const double PI=3.1415926;

ull dp[1000010];
ull jie[1000010];
//string a,b;
int main()
{
    int n;
    dp[1]=1;//前n阶数之和
    jie[1]=1;
    while(cin>>n)
    {
        int t=n;
        int i;
        if(dp[n]==0)
        {
            while(dp[n]==0)
            {
                n--;
            }
            //cout<<n<<endl;
            for( i=n+1;i<=t;i++)
            {
                jie[i]=((jie[i-1]%MOD)*(i%MOD))%MOD;
                dp[i]=(dp[i-1]+jie[i])%MOD;
                //cout<<jie[i]<<" "<<dp[i]<<endl;
            }
        }
        /*
        ull sum=0;
        for(int i=1;i<=n;i++)
        {
            if(dp[i]==0)
            {
                dp[i]=((dp[i-1]%MOD)*(i%MOD))%MOD;
            }
            sum=(sum+dp[i])%MOD;
        }
        */
        cout<<dp[t]<<endl;
    }

    return 0;
}