[联合集训6-12] Mirror 状压DP_[联合集训6-26]-CSDN博客

本文链接：https://blog.csdn.net/DOFYPXY/article/details/80699030

因为光路是可逆的，不闭合的光路都会被算进答案，所以可以把问题转化为最大化闭合光路长度，答案就是 $4nm$ 减掉其两倍。
首先考虑如果两个闭合光路共用一个镜子，那么把这个镜子拿走答案不变，所以存在一种最优方案使得每行每列的镜子都是偶数个，而且一定是相邻两个配对。于是状压每一行镜子的奇偶性，并记录当前镜子数量，一列一列填镜子即可。
因为镜子数量最多为 $2(n+m-2)$ ，复杂度可以做到 $O((n+m)m2^{2n})$
代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#define in(x,S) ((S>>x)&1)
#define chkmax(a,b) a=max(a,b)
#define chkmin(a,b) a=min(a,b)

using namespace std;
int n,m,ca,f[10][40][40],can[10],b[8][10],ans[40];
char cs[10];
bool o[40];
int popc(int x)
{
    int re=0;
    for(;x;x>>=1) re+=(x&1);
    return re;
}
int main()
{
    scanf("%d",&ca);
    for(int s=0;s<32;s++)
        o[s]=(popc(s)&1)^1;
    while(ca--)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++)
        {
            scanf("%s",cs+1);
            for(int j=1;j<=m;j++)
                b[i][j]=(cs[j]-'0');
        }
        memset(can,0,sizeof(can));
        for(int j=1;j<=m;j++)
            for(int i=0;i<n;i++)
                can[j]|=(b[i][j]<<i);       
        memset(f,0xc0,sizeof(f));
        f[0][0][0]=0;
        for(int j=0;j<m;j++)
            for(int s=0;s<(1<<n);s++)
                for(int k=0;k<=(n+m-2)*2;k++)
                    if(f[j][s][k]>=0)
                        for(int t=0;t<(1<<n);t++)
                            if(o[t]&&(t&can[j+1])==t)
                            {
                                int u=s^t,add=0,lst=-1,kk=k;
                                for(int i=0;i<n;i++)
                                    if(in(i,t))
                                    {
                                        if(lst==-1) lst=i;
                                        else add+=i-lst-1,lst=-1;       
                                        add++,kk++;
                                    }   
                                    else add+=in(i,u);
                                    chkmax(f[j+1][u][kk],f[j][s][k]+add);
                            }                   
        memset(ans,0x3f,sizeof(ans));
        for(int i=0;i<=(n+m-2)*2;i++)           
            chkmin(ans[i],4*n*m-f[m][0][i]*2);
        for(int i=1;i<=n*m;i++)
            chkmin(ans[i],ans[i-1]);
        for(int i=0;i<=n*m;i++)
            printf("%d ",ans[i]);
        puts("");                           
    }
    return 0;
}