CF 375C Circling Round Treasures [DP(spfa) 状压 射线法]

C - Circling Round Treasures

题意：

在一个$n*m$的地图上，有一些障碍，还有a个宝箱和b个炸弹。你从(sx，sy)出发，走四连通的格子。你需要走一条闭合的路径，可以自交，且围出来的复杂多边形内不能包含任何炸弹。
你围出来的复杂多边形中包含的宝箱的价值和减去步数就是你的收益。
求最大收益。$n,m \le 20，a + b \le 8$

---

太坑了课件上的题意有问题一开始没说步数...

后来改上然后$inq[][][]$最后没从$2$改成$S$然后我就看着玄学的$inq$好长时间

和上题一样，只不过从一个点变成多个点，状压一下就好了、

更新答案的时候枚举状态计算是否有炸弹和价值和

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int N=,S=<<,M=N*N*S+,INF=1e9;
inline int read(){
    char c=getchar();int x=,f=;
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x*f;
}
int n,m,g[N][N],sx,sy,val[N];
struct Object{
    int x,y,id;
    Object(int a=,int b=,int c=):x(a),y(b),id(c){}
    bool operator <(const Object &r)const{return id<r.id;}
}a[N];
int p,t;
int dx[]={,-,, },
    dy[]={, ,,-};
char s[N];
int d[N][N][S];
struct Grid{
    int x,y,s;
    Grid(int a=,int b=,int c=):x(a),y(b),s(c){}
}q[M];
int head,tail,inq[N][N][S];
inline void lop(int &x){if(x==M) x=;}
inline bool isInter(int a,int b,int x1,int y1,int x2,int y2){//1-->2 isInter (a,b)--right-->
    if(x1<a&&x2==a&&y2>b) return ;
    if(x2<a&&x1==a&&y1>b) return ;
    return ;
}
int HandleInter(Grid u,Grid t){
    int s=u.s;
    for(int i=;i<;i++)
        s^= isInter(a[i].x , a[i].y , u.x , u.y , t.x , t.y)<<i;
    return s;
}
void spfa(){
    d[sx][sy][]=;
    head=tail=;
    q[tail++]=Grid(sx,sy,);inq[sx][sy][]=;
    while(head!=tail){
        Grid u=q[head++];lop(head);
        int x=u.x,y=u.y,s=u.s;//printf("\nnow %d %d %d\n",x,y,p);
        inq[x][y][s]=;
        for(int i=;i<;i++){
            int nx=x+dx[i],ny=y+dy[i];
            if(nx<||nx>m||ny<||ny>n||g[nx][ny]) continue;

            Grid t(nx,ny);
            int ns=t.s=HandleInter(u,t);//printf("lok %d %d %d\n",nx,ny,np);
            if(d[nx][ny][ns]>d[x][y][s]+){//srintf("new %d %d %d\n",nx,ny,ns);
                d[nx][ny][ns]=d[x][y][s]+;
                if(!inq[nx][ny][ns])
                    q[tail++]=t,lop(tail),inq[nx][ny][ns]=;
            }
        }
    }
    int ans=;
    for(int i=;i<S;i++) if(d[sx][sy][i]<INF){
        int _=,bomb=;
        for(int j=;j<;j++) if(i&(<<j)){
            if(a[j].id==) {bomb=;continue;}
            else _+=val[j];
        }
        if(bomb) continue;
        ans=max(ans,_-d[sx][sy][i]);
    }
    printf("%d",ans);
}
int main(){
    //freopen("in","r",stdin);
    m=read();n=read();
    memset(d,,sizeof(d));
    for(int i=;i<=m;i++){
        scanf("%s",s+);
        for(int j=;j<=n;j++){
            g[i][j]=(s[j]!='.');
            if(s[j]=='S') g[i][j]=,sx=i,sy=j;
            else if(s[j]=='B') a[p++]=Object(i,j,);
            else if(s[j]>=''&&s[j]<='') a[p++]=Object(i,j,s[j]-''-),t++;
        }
    }
    sort(a,a+p);
    for(int i=;i<t;i++) val[i]=read();
    spfa();
}