CSPS模拟88-91
感觉自己好菜啊，考得一次不如一次了。。。压力好大，++滚粗感。
模拟88。
T1，sbt，发现离散化后数据范围变为6000，直接跑暴力即可。%%%机房众神斜率优化。
T2，大模拟，考场上只会乱搞骗分。本人菜鸡，只会大力分类讨论。。。
1 #include<bits/stdc++.h> 2 using namespace std; 3 const double eps=5e-5; 4 double f[100][100],f1[100][100],f2[100][100]; 5 double dp[110][10][10][10][10]; 6 double ans[6][10]; 7 int a[6][10],tt[10]; 8 int n; 9 char s1[20],s2[20],s3[20]; 10 double gg; 11 inline void prans(int id) 12 { 13 double gg=1.0; 14 memset(ans,0,sizeof(ans)); 15 for(int j=1;j<=8;++j) 16 for(int k=1;k<=8;++k) 17 for(int o=1;o<=8;++o) 18 for(int l=1;l<=8;++l) 19 { 20 gg-=dp[id][j][k][o][l]; 21 ans[1][a[1][j]]+=dp[id][j][k][o][l]; 22 ans[2][a[2][k]]+=dp[id][j][k][o][l]; 23 ans[3][a[3][o]]+=dp[id][j][k][o][l]; 24 ans[4][a[4][l]]+=dp[id][j][k][o][l]; 25 } 26 printf("%.2lf\n",gg*100); 27 for(int i=1;i<=4;++i) 28 { 29 for(int j=1;j<=8;++j) 30 { 31 printf("%.2lf ",ans[i][j]*100); 32 } 33 puts(""); 34 } 35 } 36 inline void init() 37 { 38 f[1][0]=f[1][1]=f[1][2]=1.0/3; 39 for(int i=1;i<=80;++i){ 40 for(int j=0;j<=80;++j){ 41 if(j)f1[i][j]=f1[i][j-1]; 42 f1[i][j]+=f[i][j]; 43 f[i+1][j]+=f[i][j]*f[1][0]; 44 f[i+1][j+1]+=f[i][j]*f[1][1]; 45 f[i+1][j+2]+=f[i][j]*f[1][2]; 46 } 47 for(int j=80;~j;--j)f2[i][j]=f2[i][j+1]+f[i][j]; 48 } 49 } 50 void work1(int m,int i,int num,int isw,double g=1.0,int fj=1,int fk=1,int fo=1,int fl=1,int jj=8,int kk=8,int oo=8,int ll=8) 51 { 52 if(m==1) 53 { 54 if(isw) 55 { 56 for(int j=fj;j<=jj;++j) 57 for(int k=fk;k<=kk;++k) 58 for(int o=fo;o<=oo;++o) 59 for(int l=fl;l<=ll;++l) 60 for(int t=0;t<=2*num;++t) 61 dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[num][t]*g; 62 } 63 else 64 { 65 for(int j=fj;j<=jj;++j) 66 for(int k=fk;k<=kk;++k) 67 for(int o=fo;o<=oo;++o) 68 for(int l=fl;l<=ll;++l) 69 dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*g; 70 } 71 } 72 else if(m==2) 73 { 74 if(isw) 75 { 76 for(int j=fj;j<=jj;++j) 77 for(int k=fk;k<=kk;++k) 78 for(int o=fo;o<=oo;++o) 79 for(int l=fl;l<=ll;++l) 80 for(int t=0;t<=2*num;++t) 81 dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[num][t]*g; 82 } 83 else 84 { 85 for(int j=fj;j<=jj;++j) 86 for(int k=fk;k<=kk;++k) 87 for(int o=fo;o<=oo;++o) 88 for(int l=fl;l<=ll;++l) 89 dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*g; 90 } 91 } 92 else if(m==3) 93 { 94 if(isw) 95 { 96 for(int j=fj;j<=jj;++j) 97 for(int k=fk;k<=kk;++k) 98 for(int o=fo;o<=oo;++o) 99 for(int l=fl;l<=ll;++l) 100 for(int t=0;t<=2*num;++t) 101 dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[num][t]*g; 102 } 103 else 104 { 105 for(int j=fj;j<=jj;++j) 106 for(int k=fk;k<=kk;++k) 107 for(int o=fo;o<=oo;++o) 108 for(int l=fl;l<=ll;++l) 109 dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*g; 110 } 111 } 112 else 113 { 114 if(isw) 115 { 116 for(int j=fj;j<=jj;++j) 117 for(int k=fk;k<=kk;++k) 118 for(int o=fo;o<=oo;++o) 119 for(int l=fl;l<=ll;++l) 120 for(int t=0;t<=2*num;++t) 121 dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[num][t]*g; 122 } 123 else 124 { 125 for(int j=fj;j<=jj;++j) 126 for(int k=fk;k<=kk;++k) 127 for(int o=fo;o<=oo;++o) 128 for(int l=fl;l<=ll;++l) 129 dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*g; 130 } 131 } 132 } 133 //////////////////////////////////////////////////分界线work1&2/////////////////////////////////////////////////////////////// 134 void work2(int m,int i,int num,int isw,double g=1.0,int fj=1,int fk=1,int fo=1,int fl=1,int jj=8,int kk=8,int oo=8,int ll=8) 135 { 136 if(m==1) 137 { 138 if(isw) 139 { 140 for(int j=fj;j<=jj;++j) 141 for(int k=fk;k<=kk;++k) 142 for(int o=fo;o<=oo;++o) 143 for(int l=fl;l<=ll;++l) 144 for(int t=0;t<=2*num;++t) 145 dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[num][t]*g; 146 } 147 else 148 { 149 for(int j=fj;j<=jj;++j) 150 for(int k=fk;k<=kk;++k) 151 for(int o=fo;o<=oo;++o) 152 for(int l=fl;l<=ll;++l) 153 dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*g; 154 } 155 } 156 else if(m==2) 157 { 158 if(isw) 159 { 160 for(int j=fj;j<=jj;++j) 161 for(int k=fk;k<=kk;++k) 162 for(int o=fo;o<=oo;++o) 163 for(int l=fl;l<=ll;++l) 164 for(int t=0;t<=2*num;++t) 165 dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[num][t]*g; 166 } 167 else 168 { 169 for(int j=fj;j<=jj;++j) 170 for(int k=fk;k<=kk;++k) 171 for(int o=fo;o<=oo;++o) 172 for(int l=fl;l<=ll;++l) 173 dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*g; 174 } 175 } 176 else if(m==3) 177 { 178 if(isw) 179 { 180 for(int j=fj;j<=jj;++j) 181 for(int k=fk;k<=kk;++k) 182 for(int o=fo;o<=oo;++o) 183 for(int l=fl;l<=ll;++l) 184 for(int t=0;t<=2*num;++t) 185 dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[num][t]*g; 186 } 187 else 188 { 189 for(int j=fj;j<=jj;++j) 190 for(int k=fk;k<=kk;++k) 191 for(int o=fo;o<=oo;++o) 192 for(int l=fl;l<=ll;++l) 193 dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*g; 194 } 195 } 196 else 197 { 198 if(isw) 199 { 200 for(int j=fj;j<=jj;++j) 201 for(int k=fk;k<=kk;++k) 202 for(int o=fo;o<=oo;++o) 203 for(int l=fl;l<=ll;++l) 204 for(int t=0;t<=2*num;++t) 205 dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[num][t]*g; 206 } 207 else 208 { 209 for(int j=fj;j<=jj;++j) 210 for(int k=fk;k<=kk;++k) 211 for(int o=fo;o<=oo;++o) 212 for(int l=fl;l<=ll;++l) 213 dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*g; 214 } 215 } 216 } 217 218 int main() 219 { 220 // freopen("da.in","r",stdin); 221 init(); 222 for(int i=1,x;i<=4;++i){ 223 scanf("%d%d",&a[i][8],&tt[i]); 224 for(int j=8;j;--j){a[i][j-1]=a[i][j]/10;a[i][j]%=10;} 225 } 226 dp[1][tt[1]][tt[2]][tt[3]][tt[4]]=1.0; 227 scanf("%d",&n); 228 /* 229 t1:当前处理点 230 t2:受影响点 231 i:当前处理的操作 232 l2:数字串长度，对？处理 233 tag:标记是否有问号 234 tg2:标记是否有= 235 */ 236 for(int i=1,t1,t2,l2,num,tag,tg2;i<=n;++i)//work2(int m,int i,int num,int isw) 237 { 238 // if(i==5)prans(i); 239 t1=l2=tag=num=tg2=0; 240 scanf("%s%s",s1,s2+1); 241 l2=strlen(s2+1); 242 if(s1[1]=='i')t1=1; 243 else if(s1[1]=='p')t1=2; 244 else if(s1[1]=='a')t1=3; 245 else t1=4; 246 switch(s2[1]) 247 { 248 case '<': 249 { 250 scanf("%d",&tg2); 251 if(l2==2)++tg2; 252 scanf("%s%s",s1,s2+1); 253 if(s1[1]=='i')t2=1; 254 else if(s1[1]=='p')t2=2; 255 else if(s1[1]=='a')t2=3; 256 else t2=4; 257 258 l2=strlen(s2+1); 259 for(int j=2;j<l2;++j)num=num*10+s2[j]-'0'; 260 if(s2[l2]=='?')tag=1; 261 else num=num*10+s2[l2]-'0'; 262 for(int j=1;j<=8;++j) 263 for(int k=1;k<=8;++k) 264 for(int o=1;o<=8;++o) 265 for(int l=1;l<=8;++l){ 266 if(t1==1) 267 { 268 for(int tmd=0;tmd<=2*a[1][j];++tmd) 269 { 270 if(tmd>=tg2) 271 { 272 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 273 continue; 274 } 275 if(t2==1) 276 { 277 if(tag){ 278 for(int t=0;t<=num*2;++t){ 279 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 280 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 281 } 282 } 283 else 284 { 285 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 286 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 287 } 288 } 289 else if(t2==2) 290 { 291 if(tag){ 292 for(int t=0;t<=num*2;++t){ 293 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 294 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 295 } 296 } 297 else 298 { 299 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 300 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 301 } 302 } 303 else if(t2==3) 304 { 305 if(tag){ 306 for(int t=0;t<=num*2;++t){ 307 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 308 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 309 } 310 } 311 else 312 { 313 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 314 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 315 } 316 } 317 else 318 { 319 if(tag){ 320 for(int t=0;t<=num*2;++t){ 321 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 322 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 323 } 324 } 325 else 326 { 327 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 328 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 329 } 330 } 331 } 332 } 333 else if(t1==2) 334 { 335 for(int tmd=0;tmd<=2*a[2][k];++tmd) 336 { 337 if(tmd>=tg2) 338 { 339 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 340 continue; 341 } 342 if(t2==1) 343 { 344 if(tag){ 345 for(int t=0;t<=num*2;++t){ 346 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 347 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 348 } 349 } 350 else 351 { 352 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 353 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 354 } 355 } 356 else if(t2==2) 357 { 358 if(tag){ 359 for(int t=0;t<=num*2;++t){ 360 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 361 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 362 } 363 } 364 else 365 { 366 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 367 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 368 } 369 } 370 else if(t2==3) 371 { 372 if(tag){ 373 for(int t=0;t<=num*2;++t){ 374 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 375 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 376 } 377 } 378 else 379 { 380 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 381 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 382 } 383 } 384 else 385 { 386 if(tag){ 387 for(int t=0;t<=num*2;++t){ 388 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 389 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 390 } 391 } 392 else 393 { 394 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 395 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 396 } 397 } 398 } 399 } 400 else if(t1==3) 401 { 402 for(int tmd=0;tmd<=2*a[3][o];++tmd) 403 { 404 if(tmd>=tg2) 405 { 406 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 407 continue; 408 } 409 if(t2==1) 410 { 411 if(tag){ 412 for(int t=0;t<=num*2;++t){ 413 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 414 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 415 } 416 } 417 else 418 { 419 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 420 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 421 } 422 } 423 else if(t2==2) 424 { 425 if(tag){ 426 for(int t=0;t<=num*2;++t){ 427 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 428 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 429 } 430 } 431 else 432 { 433 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 434 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 435 } 436 } 437 else if(t2==3) 438 { 439 if(tag){ 440 for(int t=0;t<=num*2;++t){ 441 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 442 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 443 } 444 } 445 else 446 { 447 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 448 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 449 } 450 } 451 else/////////////////////////////////////////////////////////////////here 452 { 453 if(tag){ 454 for(int t=0;t<=num*2;++t){ 455 // cout<<t<<endl; 456 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 457 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 458 } 459 } 460 else 461 { 462 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 463 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 464 } 465 } 466 } 467 } 468 else 469 { 470 for(int tmd=0;tmd<=2*a[4][l];++tmd) 471 { 472 if(tmd>=tg2) 473 { 474 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 475 continue; 476 } 477 if(t2==1) 478 { 479 if(tag){ 480 for(int t=0;t<=num*2;++t){ 481 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 482 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 483 } 484 } 485 else 486 { 487 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 488 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 489 } 490 } 491 else if(t2==2) 492 { 493 if(tag){ 494 for(int t=0;t<=num*2;++t){ 495 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 496 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 497 } 498 } 499 else 500 { 501 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 502 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 503 } 504 } 505 else if(t2==3) 506 { 507 if(tag){ 508 for(int t=0;t<=num*2;++t){ 509 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 510 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 511 } 512 } 513 else 514 { 515 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 516 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 517 } 518 } 519 else 520 { 521 if(tag){ 522 for(int t=0;t<=num*2;++t){ 523 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 524 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 525 } 526 } 527 else 528 { 529 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 530 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 531 } 532 } 533 } 534 } 535 } 536 break; 537 } 538 case '>': 539 { 540 scanf("%d",&tg2); 541 if(l2==2)--tg2; 542 scanf("%s%s",s1,s2+1); 543 if(s1[1]=='i')t2=1; 544 else if(s1[1]=='p')t2=2; 545 else if(s1[1]=='a')t2=3; 546 else t2=4; 547 548 l2=strlen(s2+1); 549 for(int j=2;j<l2;++j)num=num*10+s2[j]-'0'; 550 if(s2[l2]=='?')tag=1; 551 else num=num*10+s2[l2]-'0'; 552 553 for(int j=1;j<=8;++j) 554 for(int k=1;k<=8;++k) 555 for(int o=1;o<=8;++o) 556 for(int l=1;l<=8;++l) 557 if(t1==1) 558 { 559 for(int tmd=0;tmd<=2*a[1][j];++tmd) 560 { 561 if(tmd<=tg2) 562 { 563 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 564 continue; 565 } 566 if(t2==1) 567 { 568 if(tag){ 569 for(int t=0;t<=num*2;++t){ 570 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 571 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 572 } 573 } 574 else 575 { 576 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 577 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 578 } 579 } 580 else if(t2==2) 581 { 582 if(tag){ 583 for(int t=0;t<=num*2;++t){ 584 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 585 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 586 } 587 } 588 else 589 { 590 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 591 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 592 } 593 } 594 else if(t2==3) 595 { 596 if(tag){ 597 for(int t=0;t<=num*2;++t){ 598 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 599 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 600 } 601 } 602 else 603 { 604 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 605 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 606 } 607 } 608 else 609 { 610 if(tag){ 611 for(int t=0;t<=num*2;++t){ 612 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 613 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]*f[num][t]; 614 } 615 } 616 else 617 { 618 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 619 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[1][j]][tmd]; 620 } 621 } 622 } 623 } 624 else if(t1==2) 625 { 626 for(int tmd=0;tmd<=2*a[2][k];++tmd) 627 { 628 if(tmd<=tg2) 629 { 630 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 631 continue; 632 } 633 if(t2==1) 634 { 635 if(tag){ 636 for(int t=0;t<=num*2;++t){ 637 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 638 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 639 } 640 } 641 else 642 { 643 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 644 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 645 } 646 } 647 else if(t2==2) 648 { 649 if(tag){ 650 for(int t=0;t<=num*2;++t){ 651 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 652 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 653 } 654 } 655 else 656 { 657 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 658 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 659 } 660 } 661 else if(t2==3) 662 { 663 if(tag){ 664 for(int t=0;t<=num*2;++t){ 665 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 666 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 667 } 668 } 669 else 670 { 671 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 672 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 673 } 674 } 675 else 676 { 677 if(tag){ 678 for(int t=0;t<=num*2;++t){ 679 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 680 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]*f[num][t]; 681 } 682 } 683 else 684 { 685 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 686 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[2][k]][tmd]; 687 } 688 } 689 } 690 } 691 else if(t1==3) 692 { 693 for(int tmd=0;tmd<=2*a[3][o];++tmd) 694 { 695 if(tmd<=tg2) 696 { 697 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 698 continue; 699 } 700 if(t2==1) 701 { 702 if(tag){ 703 for(int t=0;t<=num*2;++t){ 704 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 705 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 706 } 707 } 708 else 709 { 710 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 711 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 712 } 713 } 714 else if(t2==2) 715 { 716 if(tag){ 717 for(int t=0;t<=num*2;++t){ 718 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 719 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 720 } 721 } 722 else 723 { 724 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 725 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 726 } 727 } 728 else if(t2==3) 729 { 730 if(tag){ 731 for(int t=0;t<=num*2;++t){ 732 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 733 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 734 } 735 } 736 else 737 { 738 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 739 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 740 } 741 } 742 else 743 { 744 745 if(tag){ 746 for(int t=0;t<=num*2;++t){ 747 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 748 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]*f[num][t]; 749 } 750 } 751 else 752 { 753 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 754 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[3][o]][tmd]; 755 } 756 } 757 } 758 } 759 else 760 { 761 for(int tmd=0;tmd<=2*a[4][l];++tmd) 762 { 763 if(tmd<=tg2) 764 { 765 dp[i+1][j][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 766 continue; 767 } 768 if(t2==1) 769 { 770 if(tag){ 771 for(int t=0;t<=num*2;++t){ 772 if(s2[1]=='+')dp[i+1][min(j+t,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 773 else dp[i+1][max(j-t,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 774 } 775 } 776 else 777 { 778 if(s2[1]=='+')dp[i+1][min(j+num,8)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 779 else dp[i+1][max(j-num,0)][k][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 780 } 781 } 782 else if(t2==2) 783 { 784 if(tag){ 785 for(int t=0;t<=num*2;++t){ 786 if(s2[1]=='+')dp[i+1][j][min(k+t,8)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 787 else dp[i+1][j][max(k-t,0)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 788 } 789 } 790 else 791 { 792 if(s2[1]=='+')dp[i+1][j][min(k+num,8)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 793 else dp[i+1][j][max(k-num,0)][o][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 794 } 795 } 796 else if(t2==3) 797 { 798 if(tag){ 799 for(int t=0;t<=num*2;++t){ 800 if(s2[1]=='+')dp[i+1][j][k][min(o+t,8)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 801 else dp[i+1][j][k][max(o-t,0)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 802 } 803 } 804 else 805 { 806 if(s2[1]=='+')dp[i+1][j][k][min(o+num,8)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 807 else dp[i+1][j][k][max(o-num,0)][l]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 808 } 809 } 810 else 811 { 812 if(tag){ 813 for(int t=0;t<=num*2;++t){ 814 if(s2[1]=='+')dp[i+1][j][k][o][min(l+t,8)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 815 else dp[i+1][j][k][o][max(l-t,0)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]*f[num][t]; 816 } 817 } 818 else 819 { 820 if(s2[1]=='+')dp[i+1][j][k][o][min(l+num,8)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 821 else dp[i+1][j][k][o][max(l-num,0)]+=dp[i][j][k][o][l]*f[a[4][l]][tmd]; 822 } 823 } 824 } 825 } 826 break; 827 } 828 case '+': 829 { 830 for(int j=2;j<l2;++j)num=num*10+s2[j]-'0'; 831 if(s2[l2]=='?')tag=1; 832 else num=num*10+s2[l2]-'0'; 833 work1(t1,i,num,tag); 834 break; 835 } 836 case '-': 837 { 838 for(int j=2;j<l2;++j)num=num*10+s2[j]-'0'; 839 if(s2[l2]=='?')tag=1; 840 else num=num*10+s2[l2]-'0'; 841 work2(t1,i,num,tag); 842 break; 843 } 844 } 845 } 846 prans(n+1); 847 }
T3，dp or 记搜，考场上像一个sb一样打了状压，其实值域只有0，1，2，3。记录每一种有多少个以及上一次的操作即可。注意打高精
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<vector> 6 #define N 13 7 using namespace std; 8 int n,a[N],sta; 9 int bk[N]; 10 struct Big_int{ 11 int a[60]; 12 inline void init(int x){while(x){a[++a[0]]=x%10;x/=10;}} 13 inline bool empty(){return !a[0];} 14 15 inline void print() 16 { 17 for(int i=a[0];i;--i)printf("%d",a[i]); 18 puts(""); 19 } 20 }dp[40][N][N][N][4]; 21 inline void add( Big_int &c,const Big_int &x,const int y) 22 { 23 if(!y)return; 24 c.a[0]=max(c.a[0],x.a[0]);int t=0; 25 for(int i=1;i<=c.a[0];++i){ 26 c.a[i]+=x.a[i]*y+t; 27 t=c.a[i]/10;c.a[i]%=10; 28 } 29 while(t)c.a[++c.a[0]]=t%10,t/=10; 30 return; 31 } 32 inline void work2() 33 { 34 int al=0,cnt=0; 35 for(int i=1;i<=n;++i)al+=a[i]; 36 dp[1][bk[1]][bk[2]+1][bk[3]-1][2].init(bk[3]); 37 dp[1][bk[1]+1][bk[2]-1][bk[3]][1].init(bk[2]); 38 dp[1][bk[1]-1][bk[2]][bk[3]][0].init(bk[1]); 39 for(int o=1;o<al;++o) 40 for(int i=0;i<=n;++i) 41 for(int j=0;j<=n;++j) 42 for(int k=0;k<=n;++k) 43 for(int p=0;p<=2;++p){ 44 if(dp[o][i][j][k][p].empty())continue; 45 // printf("dp[%d][%d][%d][%d][%d]=",o,i,j,k,p);dp[o][i][j][k][p].print(); 46 if(k)add(dp[o+1][i][j+1][k-1][2],dp[o][i][j][k][p],k); 47 if(j)add(dp[o+1][i+1][j-1][k][1],dp[o][i][j][k][p],j-(p==2)); 48 if(i)add(dp[o+1][i-1][j][k][0],dp[o][i][j][k][p],i-(p==1)); 49 } 50 dp[al][0][0][0][0].print(); 51 } 52 int main() 53 { 54 cin>>n; 55 for(int i=1;i<=n;++i)cin>>a[i],++bk[a[i]]; 56 work2();return 0; 57 }
模拟89。
T1，正解为图论，然而被各种乱搞*爆，%%%cbx大神打表发现用2，3，5，7筛即可。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<vector> 5 #include<algorithm> 6 #define N 1000050 7 using namespace std; 8 int ans[N]; 9 int n; 10 bool zhi[N]; 11 int pri[N],tot; 12 inline int Min(int x,int y){if(x>y)return y;return x;} 13 int main() 14 { 15 // freopen("da.in","r",stdin); 16 // freopen("my.out","w",stdout); 17 ans[0]=1;ans[1]=0; 18 for(int i=2;i<=1000049;++i)ans[i]=i; 19 for(int i=2;i<=1000000;++i){ 20 for(int j=i+1;j<=i+4;++j)ans[i]=Min(ans[i],ans[j]+j-i); 21 if(!zhi[i])pri[++tot]=i; 22 for(int j=1;j<=tot&&1ll*i*pri[j]<=1000000;++j) 23 { 24 ans[i*pri[j]]=Min(ans[i*pri[j]],ans[i]+pri[j]); 25 zhi[i*pri[j]]=1; 26 } 27 } 28 scanf("%d",&n); 29 printf("%d\n",ans[n]); 30 return 0; 31 }
T2，杜教筛and整除分块，考场上打线筛莫比乌斯能拿到40。
1 #include<bits/stdc++.h> 2 #define N 30000500 3 using namespace std; 4 long long n,ans; 5 bool zhi[N]; 6 int pri[N],tot; 7 char mu[N]; 8 int sum[N]; 9 unordered_map<int,int>H; 10 inline void init(int n) 11 { 12 mu[1]=sum[1]=1; 13 for(int i=2;i<=n;++i){ 14 if(!zhi[i])pri[++tot]=i,mu[i]=-1; 15 sum[i]=sum[i-1]+mu[i]; 16 for(int j=1;j<=tot&&1ll*i*pri[j]<=n;++j) 17 { 18 zhi[i*pri[j]]=1; 19 if(i%pri[j])mu[i*pri[j]]=-mu[i]; 20 else {mu[i*pri[j]]=0;break;} 21 } 22 } 23 } 24 inline int getsm(int x) 25 { 26 if(x<=30000000)return sum[x]; 27 if(H.find(x)!=H.end())return H[x]; 28 int ret=1; 29 for(int l=2,r;l<=x;l=r+1) 30 { 31 r=x/(x/l); 32 ret-=getsm(x/l)*(r-l+1); 33 } 34 return H[x]=ret; 35 } 36 inline void fenkuai() 37 { 38 int p=sqrt(n); 39 for(int l=2,r;l<=p;l=r+1){ 40 r=sqrt(n/(n/l/l)); 41 ans+=(getsm(r)-getsm(l-1))*(n/l/l); 42 } 43 } 44 int main() 45 { 46 cin>>n;ans=n; 47 init(30000000); 48 fenkuai(); 49 cout<<ans<<endl; 50 return 0; 51 }
T3，正解为线段树维护单调栈，按题意模拟打treap可以拿到60（然而考场上删除函数打炸了，只有15分）
1 #include<bits/stdc++.h> 2 #define N 200050 3 using namespace std; 4 int n,m,A,lim,w; 5 int lsh[N]; 6 int a[N]; 7 struct node{int op,x,y;}q[N]; 8 inline int py(int x){return A-x+1;} 9 struct Segment_tree{ 10 int ma[N<<2],len[N<<2],mi[N<<2],po[N<<2]; 11 void build(int g,int l,int r) 12 { 13 len[g]=r-l+1; 14 if(l==r)return; 15 const int m=l+r>>1; 16 build(g<<1,l,m);build(g<<1|1,m+1,r); 17 } 18 inline void upd(int g) 19 { 20 if(ma[g<<1]>ma[g<<1|1]){ma[g]=ma[g<<1];po[g]=po[g<<1];} 21 else{ma[g]=ma[g<<1|1];po[g]=po[g<<1|1];}mi[g]=mi[g<<1]; 22 } 23 inline int ask(int g,int l,int r) 24 { 25 if(ma[g]<=w||r<=lim)return 0; 26 if(l>lim) 27 { 28 if(mi[g]>w){w=ma[g];return len[g];} 29 const int m=l+r>>1; 30 if(ma[g<<1]>w){ 31 int ret=ask(g<<1,l,m);if(ma[g<<1|1]>ma[g<<1])w=ma[g<<1|1]; 32 return ret+len[g]-len[g<<1]; 33 } 34 else return ask(g<<1|1,m+1,r); 35 } 36 const int m=l+r>>1; 37 if(m<=lim){return ask(g<<1|1,m+1,r);} 38 else{ 39 if(ma[g<<1]<=w)return ask(g<<1|1,m+1,r); 40 int ret=ask(g<<1,l,m); 41 if(w==ma[g<<1]){ 42 if(ma[g<<1|1]>ma[g<<1])w=ma[g<<1|1]; 43 return ret+len[g]-len[g<<1]; 44 } 45 else{return ret+ask(g<<1|1,m+1,r);} 46 } 47 } 48 inline int getma(int g,int l,int r,int x,int y) 49 { 50 if(l>y||r<x)return 0; 51 if(l>=x&&r<=y)return po[g]; 52 const int m=l+r>>1; 53 const int a1=getma(g<<1,l,m,x,y),a2=getma(g<<1|1,m+1,r,x,y); 54 if(a[a1]>a[a2])return a1;return a2; 55 56 } 57 void change(int g,int l,int r,int pos,int ww) 58 { 59 if(l==r){ 60 mi[g]=ma[g]=ww;po[g]=l; 61 if(w){len[g]=1,po[g]=l;} 62 else len[g]=po[g]=0; 63 return; 64 } 65 const int m=l+r>>1; 66 if(pos<=m) 67 { 68 change(g<<1,l,m,pos,ww);upd(g); 69 lim=m;w=ma[g<<1]; 70 len[g]=len[g<<1]+ask(g<<1|1,m+1,r); 71 } 72 else 73 { 74 change(g<<1|1,m+1,r,pos,ww);upd(g); 75 lim=m;w=ma[g<<1]; 76 len[g]=len[g<<1]+ask(g<<1|1,m+1,r); 77 } 78 } 79 }tr1,tr2; 80 inline void init(){ 81 sort(lsh+1,lsh+A+1); 82 A=unique(lsh+1,lsh+A+1)-lsh-1; 83 } 84 inline int getd(int x) 85 { 86 lim=x;w=a[x]; 87 int ret=tr1.ask(1,1,A); 88 lim=py(x);w=a[x]; 89 return ret+tr2.ask(1,1,A); 90 91 } 92 int main() 93 { 94 95 scanf("%d",&n); 96 for(int i=1;i<=n;++i) 97 { 98 scanf("%d",&q[i].op); 99 if(q[i].op==0){ 100 scanf("%d%d",&q[i].x,&q[i].y); 101 lsh[++A]=q[i].x; 102 } 103 else if(q[i].op==1)scanf("%d",&q[i].x); 104 else scanf("%d%d",&q[i].x,&q[i].y); 105 } 106 init(); 107 for(int i=1,op,x,y,t;i<=n;++i) 108 { 109 // printf("\n\n\ni:%d\n",i); 110 op=q[i].op; 111 x=lower_bound(lsh+1,lsh+A+1,q[i].x)-lsh; 112 if(op==0){ 113 a[x]=q[i].y; 114 tr1.change(1,1,A,x,a[x]); 115 tr2.change(1,1,A,py(x),a[x]); 116 } 117 else if(q[i].op==1) 118 { 119 a[x]=0; 120 tr1.change(1,1,A,x,a[x]); 121 tr2.change(1,1,A,py(x),a[x]); 122 } 123 else 124 { 125 y=lower_bound(lsh+1,lsh+A+1,q[i].y)-lsh; 126 if(x>y)swap(x,y);t=tr1.getma(1,1,A,x,y); 127 printf("%d\n",getd(x)+getd(y)-2*getd(t)); 128 } 129 } 130 }
模拟90。
T1：luode dijkstra
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<queue> 6 #define N 505 7 using namespace std; 8 const int dx[]={1,-1,0,0}; 9 const int dy[]={0,0,1,-1}; 10 int a[N][N],b[N][N],n,m; 11 int sx,sy,tx,ty; 12 priority_queue<pair<int,int> >q; 13 inline void bfs(int x,int y) 14 { 15 q.push(make_pair(b[sx][sy],x*N+y)); 16 int d,xx,yy; 17 while(q.size()) 18 { 19 x=q.top().second;d=q.top().first;q.pop(); 20 y=x%N;x/=N;if(d!=b[x][y])continue; 21 for(int i=0;i<=3;++i) 22 { 23 xx=x+dx[i];yy=y+dy[i]; 24 if(xx<=0||xx>n||yy<=0||yy>m)continue; 25 if(d<=a[xx][yy]+b[xx][yy])continue; 26 b[xx][yy]=d-a[xx][yy]; 27 q.push(make_pair(b[xx][yy],xx*N+yy)); 28 } 29 } 30 } 31 inline void pr() 32 { 33 for(int i=1;i<=n;++i) 34 { 35 for(int j=1;j<=m;++j) 36 { 37 printf("%d ",b[i][j]); 38 } 39 puts(""); 40 } 41 } 42 int main() 43 { 44 scanf("%d%d",&n,&m); 45 for(int i=1;i<=n;++i) 46 for(int j=1;j<=m;++j) 47 scanf("%d",&a[i][j]); 48 scanf("%d%d",&sx,&sy); 49 scanf("%d%d%d",&b[sx][sy],&tx,&ty); 50 bfs(sx,sy); 51 printf("%d\n",b[tx][ty]); 52 // pr(); 53 }
T2：状压dp，怎么压都行，然而考场上打sb了，WA30，kuku。。
题解状压的思路还是不错的，积累一下。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<queue> 6 #define N 105 7 using namespace std; 8 const int inf=1000000007; 9 int n,m,bin[40],dp[N][1<<16],ct[1<<16],sta[N],ans=inf; 10 inline void init(int n) 11 { 12 for(int i=1;i<=bin[n];++i) 13 ct[i]=ct[i-(i&-i)]+1; 14 } 15 int main() 16 { 17 for(int i=0;i<=20;++i)bin[i]=1<<i; 18 memset(dp,0x3f,sizeof(dp)); 19 scanf("%d%d",&n,&m); 20 for(int i=1;i<=n;++i) 21 for(int j=0,x;j<m;++j){ 22 scanf("%d",&x); 23 if(x)sta[i]|=bin[j]; 24 } 25 init(m);dp[0][0]=0; 26 for(int j=0;j<n;++j) 27 { 28 for(int i=0,cnt,ok;i<bin[m];++i) 29 { 30 if(dp[j][i]>=inf)continue; 31 if(!sta[j+1]){dp[j+1][0]=min(dp[j+1][0],dp[j][i]);continue;} 32 for(int st=sta[j+1],ok;st;st=(st-1)&sta[j+1]){ 33 ok=1;cnt=ct[st]; 34 if((sta[j+1]&bin[0])&&!(st&bin[0]))continue; 35 for(int k=0;k<m-1;++k){ 36 if(sta[j+1]&bin[k])continue; 37 if((sta[j+1]&bin[k+1])&&!(st&bin[k+1])){ok=0;break;} 38 } 39 if(!ok)continue; 40 for(int k=0,s1,s2;k<m;++k){ 41 if((i&bin[k])&&(st&bin[k])) 42 { 43 for(s1=k+1;s1<m;++s1)if(!(sta[j]&bin[s1])||(i&bin[s1]))break; 44 for(s2=k+1;s2<m;++s2)if(!(sta[j+1]&bin[s2])||(st&bin[s2]))break; 45 if(s1==s2)--cnt; 46 } 47 } 48 dp[j+1][st]=min(dp[j+1][st],dp[j][i]+cnt); 49 } 50 } 51 } 52 for(int i=0,cnt,ok;i<bin[m];++i)ans=min(ans,dp[n][i]); 53 printf("%d\n",ans); 54 }
T3：%%%kx
对所有斜率进行离散化，对矩形的左边界和下边界分别维护一棵线段树，下标为斜率离散化后的下标。然后标记永久化一下就可以了。
注意各种各样的特判（y==0和x==0）。
1 #include<bits/stdc++.h> 2 #define N 100050 3 using namespace std; 4 int n,m,ls; 5 const double eps=1e-9; 6 const int inf=1000000007; 7 double lsh[N<<2]; 8 struct PA{ 9 double first; 10 int second; 11 friend bool operator <(const PA &a,const PA &b) 12 { 13 if(a.first==b.first)return a.second<b.second; 14 return a.first<b.first; 15 } 16 }; 17 inline PA mmp(int x,int y){PA a;a.first=x;a.second=y;return a;} 18 PA kkk; 19 struct node{double x,xx,y,yy;int op;}q[N]; 20 struct Segment_tree{ 21 PA mi[N<<2]; 22 inline void build(int g,int l,int r) 23 { 24 mi[g]=mmp(inf,0);if(l==r)return; 25 const int m=l+r>>1; 26 build(g<<1,l,m);build(g<<1|1,m+1,r); 27 } 28 inline PA getmi(int g,int l,int r,int pos) 29 { 30 if(l==r)return mi[g]; 31 const int m=l+r>>1; 32 PA ret; 33 if(pos<=m)ret=getmi(g<<1,l,m,pos); 34 else ret=getmi(g<<1|1,m+1,r,pos); 35 return min(ret,mi[g]); 36 } 37 inline void add(int g,int l,int r,int x,int y) 38 { 39 if(l>y||r<x)return; 40 if(l>=x&&r<=y){mi[g]=min(mi[g],kkk);return;} 41 const int m=l+r>>1; 42 add(g<<1,l,m,x,y);add(g<<1|1,m+1,r,x,y); 43 } 44 }T1,T2; 45 int main(){ 46 scanf("%d",&n); 47 for(int i=1;i<=n;++i) 48 { 49 scanf("%d",&q[i].op); 50 if(q[i].op==1){ 51 scanf("%lf%lf%lf%lf",&q[i].x,&q[i].y,&q[i].xx,&q[i].yy); 52 if(q[i].x)lsh[++ls]=q[i].y/q[i].x,lsh[++ls]=q[i].yy/q[i].x; 53 if(q[i].xx)lsh[++ls]=q[i].y/q[i].xx; 54 } 55 else{ 56 scanf("%lf%lf",&q[i].x,&q[i].y); 57 if(q[i].x)lsh[++ls]=q[i].y/q[i].x; 58 } 59 } 60 sort(lsh+1,lsh+ls+1);ls=unique(lsh+1,lsh+ls+1)-lsh-1; 61 double x,y,xx,yy;int t1,t2; 62 PA k1,k2,spj=mmp(inf,0),spj2=mmp(inf,0); 63 T1.build(1,1,ls);T2.build(1,1,ls); 64 q[0].x=inf;q[0].y=inf; 65 for(int i=1,op;i<=n;++i) 66 { 67 op=q[i].op;x=q[i].x;y=q[i].y;xx=q[i].xx;yy=q[i].yy; 68 if(op==1) 69 { 70 kkk=mmp(x,-i); 71 if(!y)spj2=min(spj2,mmp(x,-i)); 72 if(x) 73 { 74 t1=lower_bound(lsh+1,lsh+ls+1,y/x-eps)-lsh; 75 t2=lower_bound(lsh+1,lsh+ls+1,yy/x-eps)-lsh; 76 } 77 else 78 { 79 spj=min(spj,mmp(y,-i)); 80 t1=t2=ls; 81 } 82 T1.add(1,1,ls,t1,t2); 83 kkk=mmp(y,-i); 84 if(xx)t2=lower_bound(lsh+1,lsh+ls+1,y/xx-eps)-lsh; 85 else t2=ls; 86 T2.add(1,1,ls,t2,t1); 87 } 88 else 89 { 90 if(!q[i].y){printf("%d\n",-spj2.second);continue;} 91 if(!q[i].x){printf("%d\n",-spj.second);continue;} 92 t1=lower_bound(lsh+1,lsh+ls+1,y/x-eps)-lsh; 93 k1=T1.getmi(1,1,ls,t1);k2=T2.getmi(1,1,ls,t1); 94 yy=max(q[-k1.second].x*y/x,q[-k1.second].y);xx=max(q[-k2.second].x*y/x,q[-k2.second].y); 95 if(fabs(yy-xx)<1e-8){printf("%d\n",-min(k1.second,k2.second));} 96 else if(yy>xx)printf("%d\n",-k2.second); 97 else printf("%d\n",-k1.second); 98 } 99 } 100 }
模拟91。
T1，下面说一下心态问题，我当时秒切
发现不同的数在根号级别，然后它就成了sbt。开链表维护存在的权值，每次暴扫即可。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #define N 300050 6 using namespace std; 7 int n,m,tot; 8 int fa[N],sz[N]; 9 int pd[N]; 10 11 struct node{int pre,ne,cnt,sum;}a[N]; 12 int getfa(int x){return x==fa[x]?x:fa[x]=getfa(fa[x]);} 13 14 inline void del(int id) 15 { 16 pd[id]=0; 17 a[a[id].pre].ne=a[id].ne; 18 a[a[id].ne].pre=a[id].pre; 19 } 20 inline void ins2(int x,int id) 21 { 22 a[x].pre=a[id].pre; 23 a[a[id].pre].ne=x; 24 a[id].pre=x;a[x].ne=id; 25 pd[x]=1; 26 } 27 inline void ins(int x) 28 { 29 for(int i=a[0].ne;i<=n+1;i=a[i].ne) 30 if(i>x){ins2(x,i);return;} 31 } 32 33 int main() 34 { 35 // freopen("cards105.in","r",stdin);//diff -b -B cards105.out my.out 36 // freopen("my.out","w",stdout); 37 scanf("%d%d",&n,&m);tot=1; 38 for(int i=1;i<=n;++i)fa[i]=i,sz[i]=1; 39 40 a[0].ne=a[n+1].pre=1; 41 a[1].pre=0;a[1].ne=n+1; 42 a[1].cnt=n;pd[1]=1; 43 44 int op,x,y,sum; 45 long long ans; 46 while(m--) 47 { 48 scanf("%d",&op); 49 if(op==1) 50 { 51 scanf("%d%d",&x,&y); 52 if(getfa(x)==getfa(y))continue; 53 x=getfa(x);y=getfa(y); 54 --a[sz[x]].cnt; 55 if(!a[sz[x]].cnt)del(sz[x]); 56 --a[sz[y]].cnt; 57 if(!a[sz[y]].cnt)del(sz[y]); 58 59 fa[y]=x;sz[x]+=sz[y]; 60 ++a[sz[x]].cnt; 61 if(!pd[sz[x]])ins(sz[x]); 62 } 63 else 64 { 65 scanf("%d",&x); 66 for(int i=a[0].ne;i<=n;i=a[i].ne) 67 { 68 a[i].sum=a[i].cnt+a[a[i].pre].sum; 69 } 70 ans=0; 71 if(!x){ 72 ans=1ll*a[a[n+1].pre].sum*(a[a[n+1].pre].sum-1);ans>>=1; 73 } 74 else 75 { 76 for(int i=a[n+1].pre,j=a[n+1].pre;i>x;i=a[i].pre) 77 { 78 while(i-j<x)j=a[j].pre; 79 ans+=1ll*a[j].sum*a[i].cnt; 80 } 81 } 82 printf("%lld\n",ans); 83 } 84 } 85 }
T2，考场上大部分时间都扔给了T2，然而还是爆0了。
其实还是挺难的，考虑按概率dp，最后用概率乘位置得到期望。
考虑如何归并相同的数。
设dp[dep][i][j]，表示在归并的第dep层，原数组中位置为i的值在当前归并的区间内排名为j的概率，在l==r时显然有dp[dep][l][1]=1。
对于非叶子
首先，已有两个子区间各自的排名，考虑预处理一些东西辅助转移。
设f[i][j]表示某点在原区间排名为j，在新区间排名为i的概率，
f数组可以递推转移（考虑前一个位置来自相同区间or不同区间，具体细节看代码）
但我们发现只有一个f数组还不够。因为不能确定另一个子区间是否已为空
设h[i][j]表示某点在原区间排名为j，在新区间排名为i，且另一个子区间的数全部在当前数之前的概率。
转移与f数组类似（同样考虑前一个位置来自相同区间or不同区间，具体细节看代码）
然后就可以kx地转移了。
1 #include<bits/stdc++.h> 2 #define N 550 3 #define pb push_back 4 #define int long long 5 const int mod=998244353; 6 const int inv2=499122177; 7 using namespace std; 8 int n,a[N],dp[12][N][N],f[N][N],h[N][N],pd[5500],ans[N]; 9 vector<int>v[N<<2]; 10 inline void init(int n) 11 { 12 f[1][1]=inv2;h[1][1]=1; 13 for(int i=2;i<=n;++i) 14 for(int j=1;j<=i;++j){ 15 f[i][j]=(f[i-1][j-1]+f[i-1][i-j])%mod; 16 f[i][j]=1ll*inv2*f[i][j]%mod; 17 h[i][j]=(h[i-1][j-1]+f[i-1][i-j])%mod; 18 } 19 } 20 inline void merge(int g,int l,int r,int dep,int val) 21 { 22 v[g].clear(); 23 if(l==r){if(a[l]==val)dp[dep][l][1]=1,v[g].pb(l);return;} 24 const int m=l+r>>1,lc=g<<1,rc=g<<1|1;; 25 merge(lc,l,m,dep+1,val);merge(rc,m+1,r,dep+1,val); 26 for(int i=0;i<v[lc].size();++i)v[g].pb(v[lc][i]); 27 for(int j=0;j<v[rc].size();++j)v[g].pb(v[rc][j]); 28 for(int o=0,t;o<v[lc].size();++o){t=v[lc][o]; 29 for(int i=1;i<=v[lc].size();++i) 30 for(int j=0;j<=v[rc].size();++j) 31 if(j==v[rc].size())(dp[dep][t][i+j]+=1ll*dp[dep+1][t][i]*h[i+j][i]%mod)%=mod; 32 else (dp[dep][t][i+j]+=1ll*dp[dep+1][t][i]*f[i+j][i]%mod)%=mod; 33 } 34 for(int o=0,t;o<v[rc].size();++o){t=v[rc][o]; 35 for(int i=1;i<=v[rc].size();++i) 36 for(int j=0;j<=v[lc].size();++j) 37 if(j==v[lc].size())(dp[dep][t][i+j]+=1ll*dp[dep+1][t][i]*h[i+j][i]%mod)%=mod; 38 else (dp[dep][t][i+j]+=dp[dep+1][t][i]*f[i+j][i]%mod)%=mod; 39 } 40 v[lc].clear();v[rc].clear(); 41 } 42 signed main() 43 { 44 scanf("%lld",&n);init(n); 45 for(int i=1;i<=n;++i)scanf("%lld",&a[i]); 46 for(int i=1;i<=n;++i) 47 if(!pd[a[i]])pd[a[i]]=1,merge(1,1,n,1,a[i]); 48 for(int k=0;k<=1000;++k)v[k].clear(); 49 for(int i=1;i<=n;++i)v[a[i]].push_back(i); 50 for(int k=0,t=1;k<=1000;++k){ 51 for(int i=0,p;i<v[k].size();++i){p=v[k][i]; 52 for(int j=1;j<=v[k].size();++j)(ans[p]+=1ll*dp[1][p][j]*(j+t-1)%mod)%=mod; 53 } 54 t+=v[k].size(); 55 } 56 for(int i=1;i<=n;++i)printf("%lld ",ans[i]); 57 puts(""); 58 }
1 #include<bits/stdc++.h> 2 #define N 550 3 #define pb push_back 4 const int mod=998244353; 5 const double inv2=0.5; 6 using namespace std; 7 int n; 8 int a[N],b[N],c[N]; 9 double dp[12][N][N]; 10 double f[N][N]; 11 double h[N][N]; 12 double ans[N]; 13 vector<int>v[N<<2]; 14 inline void init(int n) 15 { 16 f[1][1]=inv2;h[1][1]=1; 17 for(int i=2;i<=n;++i){ 18 for(int j=1;j<=i;++j){ 19 f[i][j]=f[i-1][j-1]+f[i-1][i-j]; 20 f[i][j]= inv2*f[i][j] ; 21 h[i][j]=( h[i-1][j-1]+ f[i-1][i-j]); 22 } 23 } 24 // for(int i=1;i<=n;++i,puts("")) 25 // for(int j=1;j<=i;++j)printf("%.2lf ",f[i][j]); 26 } 27 int pd[1050]; 28 inline void merge(int g,int l,int r,int dep,int val) 29 { 30 v[g].clear(); 31 if(l==r){ 32 if(a[l]==val){ 33 dp[dep][l][1]=1; 34 v[g].pb(l); 35 }return; 36 } 37 const int m=l+r>>1,lc=g<<1,rc=g<<1|1;; 38 merge(lc,l,m,dep+1,val);merge(rc,m+1,r,dep+1,val); 39 for(int i=0;i<v[lc].size();++i)v[g].pb(v[lc][i]); 40 for(int j=0;j<v[rc].size();++j)v[g].pb(v[rc][j]); 41 for(int o=0,t;o<v[lc].size();++o){t=v[lc][o]; 42 for(int i=1;i<=v[lc].size();++i) 43 for(int j=0;j<=v[rc].size();++j) 44 if(j==v[rc].size())dp[dep][t][i+j]+=dp[dep+1][t][i]*h[i+j][i]; 45 else dp[dep][t][i+j]+=dp[dep+1][t][i]*f[i+j][i]; 46 } 47 for(int o=0,t;o<v[rc].size();++o){t=v[rc][o]; 48 for(int i=1;i<=v[rc].size();++i) 49 for(int j=0;j<=v[lc].size();++j) 50 if(j==v[lc].size())dp[dep][t][i+j]+=dp[dep+1][t][i]*h[i+j][i]; 51 else dp[dep][t][i+j]+=dp[dep+1][t][i]*f[i+j][i]; 52 }v[lc].clear();v[rc].clear(); 53 } 54 int main() 55 { 56 scanf("%d",&n);init(n); 57 for(int i=1;i<=n;++i)scanf("%d",&a[i]),b[i]=a[i]; 58 for(int i=1;i<=n;++i){ 59 if(!pd[a[i]]){pd[a[i]]=1; 60 merge(1,1,n,1,a[i]); 61 } 62 } 63 for(int k=0;k<=1000;++k)v[k].clear(); 64 for(int i=1;i<=n;++i)v[a[i]].push_back(i); 65 int t=1; 66 for(int k=0;k<=1000;++k){ 67 for(int i=0,p;i<v[k].size();++i){p=v[k][i]; 68 for(int j=1;j<=v[k].size();++j)ans[p]+=dp[1][p][j]*(j+t-1); 69 } 70 t+=v[k].size(); 71 } 72 for(int i=1;i<=n;++i)printf("%lf ",ans[i]); 73 puts(""); 74 }
T3，直接粘题接&代码，此题给我的启示：学习暴力，学习卡常，++前置，加fread，重载max或直接用if，学会乱搞，减小枚举上界，然后就可以卡常AC此题
1 #include<bits/stdc++.h> 2 #define N 1000050 3 char xch,xB[1<<15],*xS=xB,*xTT=xB; 4 #define getc() (xS==xTT&&(xTT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xTT)?0:*xS++) 5 inline int read() 6 { 7 int x=0,f=1;char ch=getc(); 8 while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getc();} 9 while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getc();} 10 return x*f; 11 } 12 using namespace std; 13 int n,K; 14 int a[N]; 15 int ans[N]; 16 int an2[N]; 17 inline void work1() 18 { 19 for(int i=1;i<=n;++i)ans[i]=-1; 20 int ma,mi,an,o; 21 for(int i=n;i>=1;--i){ 22 ma=mi=an=o=a[i]; 23 for(int k=i+1;k<=n;++k){ 24 if(a[k]>ma)ma=a[k]; 25 if(a[k]<mi)mi=a[k]; 26 an&=a[k];o|=a[k]; 27 if(o+mi-ma-an>=K)ans[i]=k-i+1; 28 } 29 for(int k=min(n,i+ans[i]-1);k>i;--k)if(ans[k]<ans[i])ans[k]=ans[i]; 30 } 31 for(int i=1;i<=n;++i)printf("%d ",ans[i]); 32 puts(""); 33 } 34 inline void work2() 35 { 36 for(int i=1;i<=n;++i)ans[i]=-1; 37 int ma,mi,an,o; 38 for(int i=n;i>=1;--i){ 39 ma=mi=an=o=a[i]; 40 const int lim=min(n,i+700); 41 for(int k=i+1;k<=lim;++k){ 42 if(a[k]>ma)ma=a[k]; 43 if(a[k]<mi)mi=a[k]; 44 an&=a[k];o|=a[k]; 45 if(o+mi-ma-an>=K)ans[i]=k-i+1; 46 } 47 for(int k=min(n,i+ans[i]-1);k>i;--k)if(ans[k]<ans[i])ans[k]=ans[i]; 48 } 49 for(int i=1;i<=n;++i)printf("%d ",ans[i]); 50 puts(""); 51 } 52 int main() 53 { 54 scanf("%d%d",&n,&K); 55 for(int i=1;i<=n;++i)a[i]=read(); 56 if(n>30000){work2();return 0;} 57 work1();return 0; 58 }
以下是正解：
正解代码：
1 #include<bits/stdc++.h> 2 #define N 1000050 3 char xch,xB[1<<15],*xS=xB,*xTT=xB; 4 #define getc() (xS==xTT&&(xTT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xTT)?0:*xS++) 5 inline int read() 6 { 7 int x=0,f=1;char ch=getc(); 8 while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getc();} 9 while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getc();} 10 return x*f; 11 } 12 using namespace std; 13 int n,K,a[N],ans[N],bin[35],lg[N]; 14 int st1[21][N],st2[21][N],st3[21][N],st4[21][N];//min max and or 15 inline void init() 16 { 17 for(int i=0;i<=30;++i)bin[i]=1<<i;lg[1]=0; 18 for(int i=2;i<=n;++i)lg[i]=lg[i>>1]+1; 19 for(int i=1;i<=n;++i)st1[0][i]=st2[0][i]=st3[0][i]=st4[0][i]=a[i]; 20 for(int j=1;j<=20;++j) 21 for(int i=1;i<=n;++i){ 22 st1[j][i]=min(st1[j-1][i],st1[j-1][i+bin[j-1]]); 23 st2[j][i]=max(st2[j-1][i],st2[j-1][i+bin[j-1]]); 24 st3[j][i]=st3[j-1][i]&st3[j-1][i+bin[j-1]]; 25 st4[j][i]=st4[j-1][i]|st4[j-1][i+bin[j-1]]; 26 } 27 } 28 inline int Min(int l,int r){return min(st1[lg[r-l+1]][l],st1[lg[r-l+1]][r-bin[lg[r-l+1]]+1]);} 29 inline int Max(int l,int r){return max(st2[lg[r-l+1]][l],st2[lg[r-l+1]][r-bin[lg[r-l+1]]+1]);} 30 inline int And(int l,int r){return st3[lg[r-l+1]][l]&st3[lg[r-l+1]][r-bin[lg[r-l+1]]+1];} 31 inline int Or(int l,int r){return st4[lg[r-l+1]][l]|st4[lg[r-l+1]][r-bin[lg[r-l+1]]+1];} 32 inline int getval(int l,int r){return Or(l,r)-And(l,r);} 33 struct list{int ne,pre,l,r;}li[N];int tot=1,kkk; 34 inline void del(int id){ 35 int a1=li[id].pre,a2=li[id].ne; 36 li[a1].ne=a2;li[a2].pre=a1;li[a2].l=li[id].l; 37 } 38 inline bool check(int l,int r){ 39 // printf("you checked:%d %d returned :%d\n",l,r,Or(l,r)+Min(l,r)-And(l,r)-Max(l,r)); 40 return Or(l,r)+Min(l,r)-And(l,r)-Max(l,r)>=K; 41 } 42 int ma[N<<2]; 43 inline void add(int g,int l,int r,int x,int y,int w){ 44 if(l>y||r<x)return; 45 if(l>=x&&r<=y){ma[g]=max(ma[g],w);return;} 46 const int m=l+r>>1; 47 add(g<<1,l,m,x,y,w);add(g<<1|1,m+1,r,x,y,w); 48 } 49 inline void bl(int g,int l,int r) 50 { 51 ma[g]=max(ma[g],ma[g>>1]); 52 if(l==r){printf("%d ",ma[g]);return;} 53 const int m=l+r>>1; 54 bl(g<<1,l,m);bl(g<<1|1,m+1,r); 55 } 56 inline void erfen(int l,int r,const int alr) 57 { 58 // printf("l:%d r:%d\n",l,r); 59 int mid; 60 while(l+1<r) 61 { 62 mid=l+r>>1; 63 if(check(mid,alr))r=mid; 64 else l=mid; 65 } 66 add(1,1,n,r,alr,alr-r+1); 67 } 68 int main() 69 { 70 scanf("%d%d",&n,&K);memset(ma,0xFF,sizeof(ma)); 71 for(int i=1;i<=n;++i)a[i]=read(); 72 init();li[0].l=li[0].r=0; 73 for(int i=1;i<=n;++i) 74 { 75 // printf("now:%d \n",i); 76 li[kkk].ne=++tot;li[tot].pre=kkk;kkk=tot; 77 li[kkk].l=li[kkk].r=i; 78 for(int j=li[kkk].pre;j;j=li[j].pre) 79 if(getval(li[j].l,i)==getval(li[li[j].ne].l,i))del(j); 80 for(int j=li[0].ne;j;j=li[j].ne) 81 if(check(li[j].r,i)){erfen(li[j].l-1,li[j].r,i);break;} 82 } 83 bl(1,1,n); 84 puts(""); 85 }
革命尚未成功，吾辈仍需努力
