BZOJ4676 Xor-Mul棋盘/BZOJ3254 Xor-Mul Chessboard

首先位运算具有独立性，那么我们可以每一位分别考虑，那么我们只需要确定这一位是0还是1

考虑f[i][j]表示第i行，状态为j的最小花费

复杂度20nm4^n

预处理两个状态之间的转移，复杂度20m4^n

我写了一发，但是还是T了

我们可以考虑逐格转移，f[d]表示当前轮廓线状态是d的最小花费，然后逐格转移，每次把轮廓线移动一格，枚举状态，对新格子上面原来是0还是1两种状态计算新格子的影响后取min即可

感觉逐格比逐行好写-_-

rank1 get

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<bitset>
#include<set>
#include<map>
using namespace std;
#define MAXN 6
#define MAXM 10010
#define MAXD 33
#define INF 1000000000
#define eps 1e-8
#define MOD 1000000007
#define ll long long
char xch,xB[1<<15],*xS=xB,*xTT=xB;
#define getc() (xS==xTT&&(xTT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xTT)?0:*xS++)
inline int read()
{
int x=0,f=1;char ch=getc();
while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getc();}
return x*f;
}
int n,m,N;
int a[MAXN][MAXM],b[MAXN][MAXM],e1[MAXN][MAXM],e2[MAXN][MAXM];
ll f[MAXD];
ll g[MAXD];
ll ans;
int w(int x,int y){
if(!y){
return 0;
}
return (x>>(y-1))&1;
}
int main(){
int I,i,j,d;
n=read();
m=read();
for(i=1;i<=n;i++){
for(j=1;j<=m;j++){
a[i][j]=read();
}
}
for(i=1;i<=n;i++){
for(j=1;j<=m;j++){
b[i][j]=read();
}
}
for(i=1;i<=n;i++){
for(j=2;j<=m;j++){
e1[i][j]=read();
}
}
for(i=1;i<=n;i++){
for(j=1;j<=m;j++){
e2[i][j]=read();
}
}
N=(1<<n);
for(I=0;I<=19;I++){
memset(f,0,sizeof(f));
for(j=1;j<=m;j++){
for(i=1;i<n;i++){
for(d=0;d<N;d++){
g[d]=min(f[d],f[d^(1<<(i-1))]+e1[i][j])+(w(d,i)^w(d,i-1))*e2[i-1][j]+(w(d,i)^w(a[i][j],I+1))*b[i][j];
}
memcpy(f,g,sizeof(f));
}
for(d=0;d<N;d++){
g[d]=min(f[d],f[d^(1<<(i-1))]+e1[i][j])+(w(d,i)^w(d,i-1))*e2[i-1][j]+(w(d,i)^w(d,1))*e2[i][j]+(w(d,n)^w(a[i][j],I+1))*b[i][j];
}
memcpy(f,g,sizeof(f));
}
ll tmp=1000000000000000000ll;
for(d=0;d<N;d++){
tmp=min(tmp,f[d]);
}
ans+=tmp*(1<<I);
}
printf("%lld\n",ans);
return 0;
}

/*
2 2
3 6
7 3
9 9
1 8
9 5
3 9
4 3

*/