记录编号 436890 评测结果 AAAAAAAAAAAAAAAAAAAA
题目名称 [NOIP 2016]组合数问题 最终得分 100
用户昵称 GravatarFisher. 是否通过 通过
代码语言 C++ 运行时间 0.218 s
提交时间 2017-08-12 22:21:01 内存使用 31.25 MiB
显示代码纯文本
  1. #include <iostream>
  2. #include <cstdio>
  3. #include <algorithm>
  4. using namespace std;
  5. const int maxn=2010;
  6. const int maxt=1e4+10;
  7. int d[maxn][maxn];
  8. int sum[maxn][maxn];
  9. int wen[maxt][3],tot;
  10. int t,n,m,k;
  11. int maxx,maxy;
  12. inline int lowbit(int x){return x&(-x);}
  13. inline void add(int x,int y){
  14. for(int i=x;i<=maxx;i+=lowbit(i)){
  15. for(int j=y;j<=maxy;j+=lowbit(j)){
  16. sum[i][j]++;
  17. }
  18. }
  19. }
  20. inline int find(int x,int y){
  21. int ans=0;
  22. for(int i=x;i>0;i-=lowbit(i)){
  23. for(int j=y;j>0;j-=lowbit(j)){
  24. ans+=sum[i][j];
  25. }
  26. }
  27. return ans;
  28. }
  29. int main(){
  30. freopen("problem.in","r",stdin);
  31. freopen("problem.out","w",stdout);
  32. scanf("%d%d",&t,&k);
  33. d[0][0]=d[1][0]=d[1][1]=1;
  34. maxx=maxy=-0x7fffffff;
  35. while(t--){
  36. scanf("%d%d",&n,&m);
  37. maxx=max(maxx,n);maxy=max(maxy,m);
  38. wen[++tot][1]=n,wen[tot][2]=m;
  39. }
  40. for(int i=2;i<=maxx;i++){
  41. for(int j=0;j<=maxy;j++){
  42. if(j-1>=0&&j<=i){
  43. d[i][j]=(d[i-1][j]+d[i-1][j-1])%k;
  44. if(d[i][j]==0)add(i,j);
  45. }
  46. if(j==0){
  47. d[i][j]=(d[i-1][j])%k;
  48. if(d[i][j]==0)add(i,j);
  49. }
  50. }
  51. }
  52. for(int i=1;i<=tot;i++){
  53. printf("%d\n",find(wen[i][1],wen[i][2]));
  54. }
  55. return 0;
  56. }