BZOJ1226学校食堂Dining(状压DP)
Posted myx12345
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ1226学校食堂Dining(状压DP)相关的知识,希望对你有一定的参考价值。
题意:见题面
思路:设dp[i,sta,k]为前i个人已经吃完,从第i人到第i+b[i]人的吃饭状况是sta,前一个吃完的人离i的距离是k(可能为负)的最小值
\[ dp[i+1,sta>>1,k-1]=min(dp[i+1,sta>>1,k-1],dp[i,sta,k]) (sta and 1=1,如果i已经吃完)\]
\[dp[i,sta+1<<l,l]=min(dp[i,sta+1<<l,l],dp[i,sta,k]+w(i+k,i+l) (第i+l人吃,前一个吃的人是i+k)\]
\[w(i,j)=b[i] or b[j]-b[i] and b[j]=b[i] xor b[j]\]
注意转移的时候边界判断,是否超过l人里面最小的容忍范围
1 const oo=2000000000; 2 var dp:array[1..1100,0..300,-8..8]of longint; 3 t,b:array[1..1100]of longint; 4 cas,v,n,i,j,k,l,ans,r:longint; 5 6 function clac(x,y:longint):longint; 7 begin 8 if x=0 then exit(0); 9 exit(t[x] xor t[y]); 10 end; 11 12 function min(x,y:longint):longint; 13 begin 14 if x<y then exit(x); 15 exit(y); 16 end; 17 18 begin 19 assign(input,‘bzoj1226.in‘); reset(input); 20 assign(output,‘bzoj1226.out‘); rewrite(output); 21 readln(cas); 22 for v:=1 to cas do 23 begin 24 read(n); 25 for i:=1 to n do read(t[i],b[i]); 26 fillchar(dp,sizeof(dp),$7f); 27 dp[1,0,-1]:=0; 28 for i:=1 to n do 29 for j:=0 to (1<<8)-1 do 30 for k:=-8 to 7 do 31 if dp[i,j,k]<oo then 32 begin 33 if j and 1>0 then 34 dp[i+1,j>>1,k-1]:=min(dp[i+1,j>>1,k-1],dp[i,j,k]) 35 else 36 begin 37 r:=oo; 38 for l:=0 to 7 do 39 if j and (1<<l)=0 then 40 begin 41 if i+l>r then break; 42 r:=min(r,i+b[i+l]+l); 43 dp[i,j+(1<<l),l]:=min(dp[i,j+(1<<l),l],dp[i,j,k]+clac(i+k,i+l)); 44 end; 45 end; 46 end; 47 ans:=oo; 48 for i:=-8 to -1 do ans:=min(ans,dp[n+1,0,i]); 49 writeln(ans); 50 end; 51 close(input); 52 close(output); 53 end.
以上是关于BZOJ1226学校食堂Dining(状压DP)的主要内容,如果未能解决你的问题,请参考以下文章
bzoj1226: [SDOI2009]学校食堂Dining 状压dp
bzoj1226: [SDOI2009]学校食堂Dining