其他分享
首页 > 其他分享> > LGP8474题解

LGP8474题解

作者:互联网

很萌萌的数数题。

考虑设 \(dp[n]\) 表示 \(n\) 的答案。

考虑对于一个长度为 \(n\) 的排列,令排列的所有元素 \(+1\),然后塞一个 \(1\) 进去。

容易发现,逆序对增加的数量和 \(1\) 塞的位置有关。如果 \(1\) 塞到 \(p[i]\),那么会增加 \(i-1\) 个逆序对。

所以就有 \(dp[n]=dp[n-1]\times(\sum_{i=1}^{n}2^{i-1})=dp[n-1]\times(2^n-1)=\prod_{i=1}^{n}(2^i-1)\)。

随便写一下就可以 \(O(n)\)。

接下来考虑 \(O(\sqrt{n}polylogn)\) 的做法。

我们知道快速阶乘算法。于是我们可以考虑像他那样分块。

设 \(B=\lfloor\sqrt{n}\rfloor\),那么如果我们求出了多项式 \(f(x)=\prod_{i=1}^{B}(2^ix-1)\) 就可以使用 chirp-Z 变换做到 \(O(\sqrt{n}\log n)\)。

考虑倍增。设 \(F_n(x)=\prod_{i=1}^{n}(2^ix-1)\),那么对于 \(F_n(x)\to F_{n+1}(x)\) 可以直接卷一项。

对于 \(F_n(x)\to F_{2n}(x)\),考虑 \(\prod_{i=n+1}^{2n}(2^ix-1)=\prod_{i=1}^{n}(2^i(2^nx)-1)\),可以对 \(F_n(x)\) 复合一个 \(2^nx\) 来解决。

复杂度 \(T(n)=T(n/2)+O(n\log n)\)。

总复杂度 \(O(\sqrt{n}\log n)\)。如果有 polylog 算法撅我。

\(O(n)\) 的代码:

#include<cstdio>
const int mod=1e9+7;
int n,ans(1);
signed main(){
	scanf("%d",&n);for(int pw2(2),i=1;i<=n;++i)ans=(pw2-1ll)*ans%mod,pw2*=2,pw2>=mod&&(pw2-=mod);printf("%d",ans);
}

标签:LGP8474,log,int,题解,sqrt,prod,dp,mod
来源: https://www.cnblogs.com/lmpp/p/16587957.html