「CF1091D」New Year and the Permutation Concatenation【排列组合】

D. New Year and the Permutation Concatenation

time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard output

Let $n$ be an integer. Consider all permutations on integers 1 to $n$ in lexicographic order, and concatenate them into one big sequence $p$. For example, if $n=3$, then $p=[1,2,3,1,3,2,2,1,3,2,3,1,3,1,2,3,2,1].$ The length of this sequence will be $n\ast n!$.

Let $1\le i\le j\le n\ast n!$ be a pair of indices. We call the sequence $(p_i,p_{i+1},\dots ,p_{j-1},p_j)$ a subarray of $p$. Its length is defined as the number of its elements, i.e., $j-i+1$. Its sum is the sum of all its elements, i.e.,${\sum }_{k=i}^j{p}_k$.

You are given $n$. Find the number of subarrays of $p$ of length $n$ having sum $\frac{n\ast (n+1)}{2}$. Since this number may be large, output it modulo $998244353$ (a prime number).

Input

The only line contains one integer $n(1\le n\le 10^6)$, as described in the problem statement.

Output

Output a single integer — the number of subarrays of length $n$ having sum $\frac{n\ast (n+1)}{2}$, modulo 998244353.

Examples

input

3

output

9

input

4

output

56

input

10

output

30052700

Note

In the first sample, there are 16 subarrays of length 3. In order of appearance, they are:

$[1,2,3],[2,3,1],[3,1,3],[1,3,2],[3,2,2],[2,2,1],[2,1,3],[1,3,2],[3,2,3],[2,3,1],[3,1,3],[1,3,1],[3,1,2],[1,2,3],[2,3,2],[3,2,1].$

Their sums are $6,6,7,6,7,5,6,6,8,6,7,5,6,6,7,6.$ As $\frac{n\ast (n+1)}{2}$=6, the answer is 9.

• 题意：

  将一个含有1～n的所有排列按字典序排序并依次首尾拼接，问子序列（连续）的个数使得其元素和为$\frac{n\ast (n+1)}{2}$

• 解法：

  • 至于为什么这些子序列的长度都是n还不会证明。。
  • 首先你得知道$c++$库函数$next\_permutation$的实现原理：

    •   /** Tips: next permuation based on the ascending order sort
        * sketch :
         * current: 3   7  6  2  5  4  3  1  .
         *                    |  |     |     |
         *          find i----+  j     k     +----end
         * swap i and k :
         *          3   7  6  3  5  4  2  1  .
         *                    |  |     |     |
         *               i----+  j     k     +----end
         * reverse j to end :
         *          3   7  6  3  1  2  4  5  .
         *                    |  |     |     |
         *          find i----+  j     k     +----end
         * */
      

    • 具体方法为：

      • 从后向前查找第一个相邻元素对$(i,j)$，并且满足$A[i]<A[j]（i<j）$

        易知，此时从j到end必然是降序。可以用反证法证明，请自行证明

      • 在$[j,end)$中寻找一个最小的$k$使其满足$A[i]<A[k]$

        由于$[j,end)$是降序的，所以必然存在一个k满足上面条件；并且可以从后向前查找第一个满足$A[i]<A[k]$关系的$k$，此时的$k$必是待找的$k$

      • 将$A[i]$与$A[k]$交换

        此时，$i$处变成比$i$大的最小元素，因为下一个全排列必须是与当前排列按照升序排序相邻的排列，故选择最小的元素替代$i$

        易知，交换后的$[j,end)$仍然满足降序排序。因为在$(k,end)$中必然小于$i$，在$[j,k)$中必然大于$k$，并且大于$i$

      • 逆置$[j,end)$

        由于此时$[j,end)$是降序的，故将其逆置。最终获得下一全排序

      • 注意：如果在步骤a)找不到符合的相邻元素对，即此时$i=begin$，则说明当前$[begin,end)$为一个降序顺序，即无下一个全排列，$STL$的方法是将其逆置成升序

  • $next\_permutation$函数找到当前序列的一个比当前序列字典序大的最小字典序的排列，也就是对应了题目要求，那么本题我们可以这么想：对于可以构成 $\frac{n\ast (n+1)}{2}$的连续子序列，他的位置有如下两种情况：
    • 恰好为第$i$个排列，可以发现所有的$n!$个排列都满足要求
    • 或者是第$i$个排列的后$k$个与第$i+1$个排列的前$n-k$个排列拼接而成
    • 现在考虑后者：显然要使这$n$个数恰好是一个排列的话，第i个排列的前$n-k$个一定与第$i+1$个排列的前$n-k$个相同，那么联系$nex{t}_{p}ermutation$的过程，不难发现，最长递减后缀长度$len<k$才能满足，接下来统计会有多少区间
    • 正难则反！直接统计不行的会方便很多：后$k$位从$n$个数中选出$k$个递减依次放置在最后棉，方案数为$C_n^k$，前面$n-k$个直接把剩下的排列就行了，方案数为$(n-k)!$，所以答案就是$$n\ast n!-\sum _{k=1}^{n-1}{C}_n^k\ast (n-k)!=n\ast n!-\sum _{k=1}^{n-1}\frac{n!}{k!}$$
    • 有大佬看出来这个题答案是下面这个式子：$$ans[n]=(ans[n-1]+(n-1)!-1)\ast n$$

• 附代码

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
const ll mod=998244353;
const int maxn=1e6+10;

int n;
ll fac[maxn];

ll quick_pow(ll a,ll b)
{
	ll res=1LL;
	while(b){
		if(b&1) res=res*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return res;
}

int main()
{
	scanf("%d",&n);fac[0]=1LL;
	for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod;

	ll ans=1LL*n*fac[n]%mod;
	for(int i=1;i<n;i++) ans=(ans-(fac[n]*quick_pow(fac[i],mod-2)%mod))%mod;
	printf("%lld\n",(ans+mod)%mod);
}