bzoj2755【SCOI2012】喵星人的入侵

输入格式

第一行为三个整数n，m，K，分别表示地图的长和宽，以及最多能放置的炮塔数量。 
接下来的n行，每行包含m个字符，‘#’表示地图上原有的障碍，‘.’表示该处为空地，
数据保证在原地图上存在S到T的路径。

---

输出格式

输出在合理布阵下，喵星人采取最优策略后，会受到的最大伤害。 
注意必须保证在布阵结束后喵星人仍然可以沿一条或以上的路径从起点S到达终点T，
否则他们组织更大规模的侵略。

---

提示

【数据范围】 
对于30%的数据，保证：
1<=N,M<=6
对于100%的数据，保证：
1<=N<=6，1<=M<=20，1<=K<=15，且从S到T的路径必定存在。

---

• 题解：

  • 由于可以无限放置障碍，所以最后一定只会有一条路径；
  • 插头$dp$求出路径即可；
  • 注意和普通路径不同，路径不能组成四连通块；

 #include<bits/stdc++.h>
 #define il inline
 #define rg register
 using namespace std;
 const int N=,M=,sz=1e5;
 int n,m,K,c1[N],c2[N],cur,ans;
 char s[M][M];
 struct HASH{
     int o,hd[sz],v[sz],w[sz],nt[sz];
     il void init(){for(int i=;i<=o;++i)hd[v[i]%sz]=,v[i]=w[i]=nt[i]=;o=;}
     il void upd(int x,int y){
         for(rg int i=hd[x%sz];i;i=nt[i])if(v[i]==x){
             if(w[i]<y)w[i]=y;
             return;
         }
         nt[++o]=hd[x%sz],hd[x%sz]=o,v[o]=x,w[o]=y;
     }
 }f[][M];
 il void upd(int&x,int y){if(x<y)x=y;}
 il void decode(int x){
     for(rg int i=;i<=m;++i)c1[i]=x&,x>>=;
     for(rg int i=;i<=m;++i)c2[i]=x&,x>>=;
 }
 il int encode(){
     int x=;
     for(rg int i=m;~i;--i)x=(x<<)^c2[i];
     for(rg int i=m;~i;--i)x=(x<<)^c1[i];
     return x;
 }
 il int le(int x){
     for(rg int i=x,y=;~i;--i)if(c1[i]&&c1[i]!=){
         y+=(c1[i]&)?:-;
         if(!y)return i;
     }return ;
 }
 il int ri(int x){
     for(rg int i=x,y=;i<=m;++i)if(c1[i]&&c1[i]!=){
         y+=(c1[i]&)?:-;
         if(!y)return i;
     }return ;
 }
 il int cal1(int x){
     int l=max(x-,),r=min(x+,m),re=;
     for(rg int i=l;i<=r;++i)if(c2[i]==)re++;
     return re;
 }
 il int cal2(int x){
     int l=max(x-,),r=min(x+,m),re=;
     for(rg int i=l;i<=r;++i)if(c2[i]==)re++;
     return re;
 }
 il bool check(int j){
     for(rg int i=;i<=m;++i)if(i!=j&&i!=j+&&c1[i])return false;
     return true;
 }
 int main(){
     #ifndef ONLINE_JUDGE
     freopen("C.in","r",stdin);
     freopen("C.out","w",stdout);
     #endif
     scanf("%d%d%d",&n,&m,&K);
     for(rg int i=;i<n;++i)scanf("%s",s[i]);
     if(n<m){
         for(rg int i=;i<n;++i)
         for(rg int j=i+;j<m;++j){
             swap(s[i][j],s[j][i]);
         }
         swap(n,m);
     }
     f[cur=][].upd(,);
     for(rg int i=;i<n;++i)
     for(rg int j=;j<m;++j){
         for(rg int k=;k<=K;++k){
             for(rg int l=;l<=f[cur][k].o;++l){
                 int w=f[cur][k].w[l];
                 decode(f[cur][k].v[l]);
                 if(!j){
                     if(c1[m])continue;
                     for(rg int i=m;i;--i)c1[i]=c1[i-],c2[i]=c2[i-];
                     c1[]=c2[]=;
                 }
                 /*
                 {
                     printf("%d %d %d:\n",i,j,k);
                     for(int i=0;i<=m;++i)printf("%d ",c1[i]);
                     puts("");
                     for(int i=0;i<=m;++i)printf("%d ",c2[i]);
                     puts("");
                     printf("%d\n",w);
                     printf("#\n");
                 }*/

                 int &p=c1[j],&q=c1[j+],&r=c2[j],w1=w+cal1(j),w2=w+cal2(j);
                 int t1=j&&(c2[j-]&),t2=j!=m&&(c2[j+]&);
                 if(s[i][j]=='#'){
                     if(p||q)continue;
                     r=,f[cur^][k].upd(encode(),w);
                 }else if(!p&&!q){
                     if(s[i][j]=='.'){
                         p=,q=;
                         r=,f[cur^][k].upd(encode(),w);
                         if(k<K)r=,f[cur^][k+].upd(encode(),w1);
                     }
                     if(t1||t2)continue;
                     if(s[i][j]=='S'||s[i][j]=='T'){
                         p=,q=,r=,f[cur^][k].upd(encode(),w2);
                         p=,q=,r=,f[cur^][k].upd(encode(),w2);
                     }else p=,q=,r=,f[cur^][k].upd(encode(),w2);
                 }else if(!p||!q){
                     if(t1&&t2)continue;
                     if(s[i][j]=='S'||s[i][j]=='T'){
                         if(p+q!=){
                             if(p+q==)c1[ri(j)]=;else c1[le(j+)]=;
                             p=q=,r=,f[cur^][k].upd(encode(),w2);
                         }
                         else if(check(j))p=,q=,r=,f[cur^][k].upd(encode(),w2);
                     }else{
                         r=,swap(p,q),f[cur^][k].upd(encode(),w2);
                         r=,swap(p,q),f[cur^][k].upd(encode(),w2);
                     }
                 }else{
                     if(p==&&q==){
                         p=q=,r=,f[cur^][k].upd(encode(),w2);
                     }else if(p==||q==){
                         if(p+q-==)c1[ri(j)]=;
                         else if(p+q-==)c1[le(j+)]=;
                         p=q=,r=,f[cur^][k].upd(encode(),w2);
                     }else {
                         if(p==&&q==)continue;
                         if(p==q){
                             if(p==)c1[ri(j+)]=;
                             else c1[le(j)]=;
                         }
                         p=q=,r=,f[cur^][k].upd(encode(),w2);
                     }
                 }
             }
             f[cur][k].init();
         }
         cur^=;
     }
     for(int i=;i<=K;++i)
     for(int j=;j<=f[cur][i].o;++j){
         decode(f[cur][i].v[j]);
         if(check(m+))upd(ans,f[cur][i].w[j]);
     }
     printf("%d\n",ans);
     return ;
 }