标签：LGP8474 log int 题解 sqrt prod dp mod

很萌萌的数数题。

考虑设 \(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

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。