傳送門啦

15分暴力，但看題解說暴力分有30分。

就是找到公式，然后套公式。。

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> using namespace std;long long read(){char ch;bool f = false;while((ch = getchar()) < '0' || ch > '9')if(ch == '-') f = true;int res = ch - 48;while((ch = getchar()) >= '0' && ch <='9')res = res * 10 - ch + 48;return f ? res + 1 : res; }long long jc(long long a){//求階乘 if(a == 0) return 1;long long ans = 1;for(int i=1;i<=a;i++)ans *= i;return ans; //b = !a } long long C(long long n,long long m){return jc(n) / (jc(m) * jc(n - m)); } //組合數公式：Cn^m = !n / (!m * !(n - m)) long long t,k,n,m; long long sum,x;int main(){t = read(); k = read();while(t--){x = 0;n = read(); m = read();//sum = jc(n) / (jc(m) * jc(n - m));for(long long i=1;i<=n;i++){//int j = min(i , m);for(long long j=1;j<=min(i,m);j++){//sum = jc(i) / (jc(j) * jc(i - j));if(C(i,j) % k == 0)x++;}}printf("%lld\n",x);}return 0; }

15分，我現在用了組合數的遞推公式，按理說應該更快了，但。。（想不通，數據范圍在那里啊）

c[i][j]即為從i件物品中選j件的方案數。如果第i件物品不選，方案數就變為c[i-1][j]，如果選第i件物品，方案數就變為c[i-1][j-1]，總方案數就為兩種情況的方案數之和

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int maxn = 2005;long long read(){char ch;bool f = false;while((ch = getchar()) < '0' || ch > '9')if(ch == '-') f = true;int res = ch - 48;while((ch = getchar()) >= '0' && ch <='9')res = res * 10 - ch + 48;return f ? res + 1 : res; }long long t,k,n,m; long long sum,x,C[maxn][maxn];int main(){t = read(); k = read();while(t--){x = 0;C[1][0] = C[1][1] = 1;n = read(); m = read();for(long long i=2;i<=n;i++){C[i][0] = 1;for(long long j=1;j<=min(i,m);j++){C[i][j] = C[i-1][j] + C[i-1][j-1];if(C[i][j] % k == 0)x++;}}printf("%lld\n",x);}return 0; }

為了提高效率，我們可以進行進一步的優化，就是預處理出組合數從而求出所有區間的滿足條件的組合數個數，這里就要用到二維前綴和

楊輝三角

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int maxn = 2005;inline int read() {int x=0,f=1;char ch=getchar();while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f; }long long t,k,n,m; long long sum[maxn][maxn],x,C[maxn][maxn];void work(){for(int i=1;i<=2000;i++){C[i][0] = 1;C[i][i] = 1;}C[1][1] = 1;for(long long i=2;i<=2000;i++)for(long long j=1;j<i;j++){C[i][j] = (C[i-1][j] + C[i-1][j-1]) % k;}for(long long i=1;i<=2000;i++){for(long long j=1;j<=i;j++){sum[i][j] = sum[i-1][j] + sum[i][j-1] - sum[i-1][j-1];if(C[i][j] == 0)sum[i][j]++ ;}sum[i][i+1] = sum[i][i]; } }int main(){memset(C,0,sizeof(C));memset(sum,0,sizeof(sum));t = read(); k = read();work();while(t--){n = read(); m = read();m = min(n , m);printf("%lld\n",sum[n][m]);}return 0; }

還有一個事不得不說，我改了一下午竟然發現是自己的快讀打錯了：

修改后：

暴力 40分

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> using namespace std;inline int read() {int x=0,f=1;char ch=getchar();while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f; }long long jc(long long a){//求階乘 if(a == 0) return 1;long long ans = 1;for(int i=1;i<=a;i++)ans *= i;return ans; //b = !a } long long C(long long n,long long m){return jc(n) / (jc(m) * jc(n - m)); } //組合數公式：Cn^m = !n / (!m * !(n - m)) long long t,k,n,m; long long sum,x;int main(){t = read(); k = read();while(t--){x = 0;n = read(); m = read();//sum = jc(n) / (jc(m) * jc(n - m));for(long long i=1;i<=n;i++){//int j = min(i , m);for(long long j=1;j<=min(i,m);j++){//sum = jc(i) / (jc(j) * jc(i - j));if(C(i,j) % k == 0)x++;}}printf("%lld\n",x);}return 0; }

遞推公式 70

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int maxn = 2005;inline int read() {int x=0,f=1;char ch=getchar();while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f; }long long t,k,n,m; long long sum,x,C[maxn][maxn];int main(){t = read(); k = read();while(t--){x = 0;C[1][0] = C[1][1] = 1;n = read(); m = read();for(long long i=2;i<=n;i++){C[i][0] = 1;for(long long j=1;j<=min(i,m);j++){C[i][j] = C[i-1][j] + C[i-1][j-1];if(C[i][j] % k == 0)x++;}}printf("%lld\n",x);}return 0; }

轉載于:https://www.cnblogs.com/Stephen-F/p/9884660.html

總結

以上是生活随笔為你收集整理的洛谷P2822组合数问题的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。