比赛 NOI2015Day1 评测结果 RRRRRRRRRR
题目名称 程序自动分析 最终得分 0
用户昵称 石家庄二中教练 运行时间 0.008 s
代码语言 C++ 内存使用 2.84 MiB
提交时间 2015-08-01 12:40:51
显示代码纯文本
  1. #include<iostream>
  2. #include<cstdio>
  3. #include<algorithm>
  4. #include<cstring>
  5. #include<bitset>
  6. using namespace std;
  7. typedef long long LL;
  8. const int SIZEN=100010;
  9. const int SIZEH=70;
  10. typedef bitset<SIZEN> BTS;
  11. int bitcnt(LL n){
  12. if(n==0) return 1;
  13. int ans=0;
  14. while(n){
  15. ans++;
  16. n/=2;
  17. }
  18. return ans;
  19. }
  20. int N,M=0;
  21. LL S;
  22. LL HA[SIZEN];
  23. bool Gauss_Jordan(BTS A[SIZEH],int h,int n,bool b[SIZEN]){
  24. static int pivot[SIZEH];
  25. int i,r=0;
  26. for(i=0;i<h;i++){
  27. while(r<n){
  28. bool flag=false;
  29. for(int k=i;k<h;k++){
  30. if(A[k][r]==1){
  31. flag=true;
  32. pivot[i]=r;
  33. swap(A[k],A[i]);
  34. break;
  35. }
  36. }
  37. if(flag) break;
  38. r++;
  39. }
  40. if(r==n) break;
  41. for(int t=0;t<h;t++){
  42. if(t==i) continue;
  43. if(A[t][r]==1) A[t]^=A[i];
  44. }
  45. }
  46. memset(b,0,sizeof(b));
  47. for(int t=i;t<h;t++) if(A[t][n]==1) return false;
  48. for(int t=0;t<i;t++) b[pivot[t]]=A[t][n];
  49. return true;
  50. }
  51. BTS A[SIZEH],A1[SIZEH];
  52. bool b[SIZEN];
  53. bool test(int h,int k,int p){
  54. memcpy(A1,A,sizeof(A));
  55. A1[h][N]=p;
  56. for(int i=0;i<N;i++) A1[h][i]=(HA[i]>>k)&1;
  57. if(Gauss_Jordan(A1,h+1,N,b)){
  58. memcpy(A,A1,sizeof(A1));
  59. return true;
  60. }
  61. return false;
  62. }
  63. void work(void){
  64. int h=0;
  65. for(int i=M-1;i>=0;i--){
  66. if(((S>>i)&1)==0){
  67. if(!test(h,i,1)) test(h,i,0);
  68. h++;
  69. }
  70. }
  71. for(int i=M-1;i>=0;i--){
  72. if(((S>>i)&1)==1){
  73. if(!test(h,i,0)) test(h,i,1);
  74. h++;
  75. }
  76. }
  77. for(int i=0;i<N;i++) printf("%d ",2-(int)b[i]);printf("\n");
  78. /*LL X1=0,X2=0;
  79. for(int i=0;i<N;i++){
  80. if(b[i]) X1^=HA[i];
  81. else X2^=HA[i];
  82. }
  83. cout<<X1+X2<<endl;*/
  84. }
  85. void read(void){
  86. scanf("%d",&N);
  87. S=0;
  88. for(int i=0;i<N;i++){
  89. scanf("%I64d",&HA[i]);
  90. M=max(M,bitcnt(HA[i]));
  91. S^=HA[i];
  92. }
  93. }
  94. int main(){
  95. //freopen("t1.in","r",stdin);
  96. read();
  97. work();
  98. return 0;
  99. }