程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> C++ >> C++入門知識 >> 1485 [HNOI2009] 有趣的數列 題解

1485 [HNOI2009] 有趣的數列 題解

編輯:C++入門知識

1485: [HNOI2009]有趣的數列 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 484  Solved: 272 [Submit][Status] Description  我們稱一個長度為2n的數列是有趣的,當且僅當該數列滿足以下三個條件:     (1)它是從1到2n共2n個整數的一個排列{ai};     (2)所有的奇數項滿足a1<a3<…<a2n-1,所有的偶數項滿足a2<a4<…<a2n;     (3)任意相鄰的兩項a2i-1與a2i(1≤i≤n)滿足奇數項小於偶數項,即:a2i-1<a2i。     現在的任務是:對於給定的n,請求出有多少個不同的長度為2n的有趣的數列。因為最後的答案可能很大,所以只要求輸出答案 mod P的值。 Input 輸入文件只包含用空格隔開的兩個整數n和P。輸入數據保證,50%的數據滿足n≤1000,100%的數據滿足n≤1000000且P≤1000000000。 Output 僅含一個整數,表示不同的長度為2n的有趣的數列個數mod P的值。 Sample Input 3 10 Sample Output 5           對應的5個有趣的數列分別為(1,2,3,4,5,6),(1,2,3,5,4,6),(1,3,2,4,5,6),(1,3,2,5,4,6),(1,4,2,5,3,6)   【轉化】就是求卡特蘭數。 【初始代碼】 [cpp]   #include<cstdio>   using namespace std;   typedef long long ll;   ll prime[200005],a[200005];   bool f[2000005];   ll temp,n,p,i,j,cnt,mod;   ll pow(ll a,ll b)   {     ll ans;     for (ans=1;b;b>>=1,a=a*a%mod)        if (b&1) ans=ans*a%mod;     return ans;   }   int main()   {     scanf("%lld%lld",&n,&mod);     for (i=2;i<=n*2;i++)     {       if (!f[i]) prime[++cnt]=i;       for (j=1;j<=cnt&&prime[j]*i<=n*2;j++)         f[prime[j]*i]=true;     }     for (i=n+2;i<=n*2;i++)       for (j=1,p=i;j<=cnt&&p;j++)         while (p%prime[j]==0) a[j]++,p/=prime[j];     for (i=2;i<=n;i++)       for (j=1,p=i;j<=cnt&&p;j++)         while (p%prime[j]==0) a[j]--,p/=prime[j];     temp=1;     for (j=1;j<=cnt;j++)       if (a[j]) temp=(temp*pow(prime[j],a[j]))%mod;     printf("%lld",temp);for (;;);     return 0;   }     用歐拉篩法,O(n)的效率求出每個質數。然後枚舉階乘,像質數表一樣把一個數給分解。但是效率很低。 【優化1】如果一個數是合數,我們可以把它的某個因子記下來。然後我們同樣從開始枚舉階乘,而且是倒著枚舉。對於每個數,如果它是合數,我就把它分解。比如,設f[n]為結果中含有n因子的個數。u是n的一個約數。那麼我們可以f[u]+=f[n],f[n/u]+=f[n]。這樣就不用多次用快速冪了。直到n是質數為止。 【優化2】開始可以把1--n的f[i]設為-1,把n+2--2*n(注意,最後要除n+1,所以從n+2開始)的f[i]設為1.這樣只需1次循環。 【AC代碼】 [cpp]   #include<cstdio>   using namespace std;   typedef long long ll;   ll prime[200005],a[2000005],come[2000005];   ll temp,n,p,i,j,cnt,mod;   ll pow(ll a,ll b)   {     ll ans;     for (ans=1;b;b=b/2,a=a*a%mod)        if (b&1) ans=ans*a%mod;     return ans;   }   int main()   {     scanf("%lld%lld",&n,&mod);     for (i=2;i<=n*2;i++)     {       if (!come[i]) prime[++cnt]=i;          for (j=1;j<=cnt&&prime[j]*i<=n*2;j++)         come[prime[j]*i]=i;     }     temp=1;     for (i=2;i<=n;i++) a[i]=-1;     for (i=n+2;i<=2*n;i++) a[i]=1;     for (i=n*2;i>1;i--)       if (come[i])       {         a[come[i]]+=a[i];         a[i/come[i]]+=a[i];       }       else         temp=temp*pow(i,a[i])%mod;              printf("%lld",temp);     return 0;   }    

  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved