bzoj2111,P2606-[ZJOI2010]排列计数【Lucas,组合计数,dp】
生活随笔
收集整理的這篇文章主要介紹了
bzoj2111,P2606-[ZJOI2010]排列计数【Lucas,组合计数,dp】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:
https://www.lydsy.com/JudgeOnline/problem.php?id=2111
https://www.luogu.org/problem/P2606
題目大意
長為nnn的序列PPP,然后要求Pi>P?i2?P_i>P_{\lfloor \frac{i}{2}\rfloor}Pi?>P?2i???。求排列個數。
解題思路
若用iii連向?i2?\lfloor \frac{i}{2}\rfloor?2i??,那么這是一棵以111為根的數,然后考慮樹形dpdpdp。
我們發現每棵子樹滿足子節點小于父節點的性質且兩棵不包含的子樹互不關聯,所以我們只需要分配編號即可
也就是fi=f2?i?f2?i+1?Csizei?1size2?if_{i}=f_{2*i}*f_{2*i+1}*C_{size_i-1}^{size_{2*i}}fi?=f2?i??f2?i+1??Csizei??1size2?i??
然后答案就是f1f_1f1?。因為在其中fif_ifi?可能會大于ppp,所以我們用LucasLucasLucas定理
時間復雜度O(nlog?n+min{n,p})O(n\log n+min\{n,p\})O(nlogn+min{n,p})
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=1e6+10; ll T,n,p,a[N],f[2*N],s[2*N]; ll power(ll x,ll b,ll p) {ll ans=1;while(b){if(b&1) ans=ans*x%p;x=x*x%p;b>>=1;}return ans; } ll C(ll n,ll m,ll p) {if(m>n) return 0ll;return a[n]*power(a[m],p-2,p)%p*power(a[n-m],p-2,p)%p; } ll lucas(ll n,ll m,ll p) {if(!m) return 1ll;return lucas(n/p,m/p,p)*C(n%p,m%p,p)%p; } int main() {scanf("%lld%lld",&n,&p);a[0]=1;for(ll i=1;i<=min(n,p);i++)a[i]=a[i-1]*i%p;for(ll i=n;i>=1;i--){s[i]=1+s[i*2]+s[i*2+1];if(i*2+1<=n)f[i]=f[i*2]*f[i*2+1]%p*lucas(s[i]-1,s[i*2],p)%p;else if(i*2<=n) f[i]=f[i*2];else f[i]=1; }printf("%lld",f[1]);return 0; }總結
以上是生活随笔為你收集整理的bzoj2111,P2606-[ZJOI2010]排列计数【Lucas,组合计数,dp】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 遍身罗绮者不是养蚕人是什么意思 遍身罗绮
- 下一篇: 学习手记(2019/7/05~2019/