NOIP 2016 提高组 复赛 第二天 第一题 组合数问题 problem 整除意味余数为0 杨辉三角 AC代码(二维前缀和(矩阵和))+95分代码+70分代码+55代码
总目录详见:NOIP 提高组 复赛 试题 目录 信奥 历年
在线测评地址:https://www.luogu.com.cn/problem/P2822
1.AC代码(杨辉三角+二维前缀和(矩阵和))
二维前缀和(矩阵和) 公式推导如下:
AC代码(杨辉三角+二维前缀和(矩阵和)) 如下:
#include <cstdio>
#include <algorithm>
#define maxn 2010
using namespace std;
int c[maxn][maxn],cnt[maxn][maxn];
int main(){
int t,k,n,m,i,j,b;
scanf("%d%d",&t,&k);
for(i=0;i<=2000;i++)
for(j=0;j<=2000;j++)
c[i][j]=-1;
for(i=0;i<=2000;i++)c[i][0]=c[i][i]=1;
for(i=1;i<=2000;i++)
for(j=1;j<i;j++)
c[i][j]=(c[i-1][j-1]+c[i-1][j])%k;//杨辉三角
for(i=0;i<=2000;i++)
for(j=0;j<=2000;j++){
cnt[i][j]=cnt[i-1][j]+cnt[i][j-1]-cnt[i-1][j-1];//二维前缀和(矩阵和)
if(c[i][j]==0)cnt[i][j]+=1;
}
while(t--){
scanf("%d%d",&n,&m);
b=min(n,m);
printf("%d\n",cnt[n][b]);
}
return 0;
}
2.95分代码(杨辉三角)
95分代码(杨辉三角)如下:
#include <cstdio>
#include <algorithm>
#define maxn 2010
using namespace std;
int c[maxn][maxn];
int main(){
int t,k,n,m,i,j,cnt,b;
scanf("%d%d",&t,&k);
for(i=0;i<=2000;i++)c[i][0]=c[i][i]=1;
for(i=1;i<=2000;i++)
for(j=1;j<i;j++)
c[i][j]=(c[i-1][j-1]+c[i-1][j])%k;//杨辉三角
while(t--){
scanf("%d%d",&n,&m);
cnt=0;
for(i=0;i<=n;i++){
b=min(i,m);
for(j=0;j<=b;j++)
if(c[i][j]==0)
cnt++;
}
printf("%d\n",cnt);
}
return 0;
}
3.70分代码(long long)
70分代码(long long)代码如下:
#include <cstdio>
#include <algorithm>
#define maxn 2010
#define LL long long
using namespace std;
LL c[maxn][maxn];
int main(){
int t,k,n,m,i,j,cnt,b;
scanf("%d%d",&t,&k);
for(i=0;i<=2000;i++)c[i][0]=c[i][i]=1;
for(i=1;i<=2000;i++)
for(j=1;j<i;j++)
c[i][j]=(c[i-1][j-1]+c[i-1][j]);//杨辉三角
while(t--){
scanf("%d%d",&n,&m);
cnt=0;
for(i=0;i<=n;i++){
b=min(i,m);
for(j=0;j<=b;j++)
if(c[i][j]%k==0)
cnt++;
}
printf("%d\n",cnt);
}
return 0;
}
4.55分代码(int)
55分代码(int) 如下:
#include <cstdio>
#include <algorithm>
#define maxn 2010
using namespace std;
int c[maxn][maxn];
int main(){
int t,k,n,m,i,j,cnt,b;
scanf("%d%d",&t,&k);
for(i=0;i<=2000;i++)c[i][0]=c[i][i]=1;
for(i=1;i<=2000;i++)
for(j=1;j<i;j++)
c[i][j]=(c[i-1][j-1]+c[i-1][j]);//杨辉三角
while(t--){
scanf("%d%d",&n,&m);
cnt=0;
for(i=0;i<=n;i++){
b=min(i,m);
for(j=0;j<=b;j++)
if(c[i][j]%k==0)
cnt++;
}
printf("%d\n",cnt);
}
return 0;
}