題意：一個$n\times m$的數表，數值$\in[0,4)$，你可以任意次選擇一行或一列$+1,\text{mod }4$，要最小化所有數的和
因為$n\leq10$，所以數表可以看成$m$個$n$位$4$進制數$a_{1\cdots m}$，以下使用不進位加法
定義$f(x)=\min\limits_{i=0}^3\left(x+(i\cdots i)_4\right)$，如果確定了行的方案$s$，答案就是$\sum\limits_if(a_i+s)$
記$c_x$為$x$在$a_{1\cdots m}$中的出現次數，$\text{ans}(s)=\sum\limits_xc_xf(x+s)=\sum\limits_{a+b=s}c_{-a}f(b)$，是一個高維循環卷積
高維卷積需要用到高維FFT，根據wiki，只需在每一維分別dft即可，時間復雜度$O(n4^n)$
#include<stdio.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int mod=998244353,inf=2147483647;
int mul(int a,int b){return(ll)a*b%mod;}
int ad(int a,int b){return(a+=b)>=mod?a-mod:a;}
int pow(int a,int b){int s=1;while(b){if(b&1)s=mul(s,a);a=mul(a,a);b>>=1;}return s;
}
int*w[5],iN;
void pre(){int i,j,t;for(i=2;i<=4;i<<=1){w[i]=new int[i+1];t=pow(3,(mod-1)/i);w[i][0]=1;for(j=1;j<=i;j++)w[i][j]=mul(w[i][j-1],t);}iN=pow(4,mod-2);
}
void dft4(int*a,int sh,int on){#define a(i) a[(i)<<sh]const int N=4;int i,j,k,t;swap(a(1),a(2));for(i=2;i<=N;i<<=1){for(j=0;j<N;j+=i){for(k=0;k<i>>1;k++){t=mul(a(i/2+j+k),w[i][on==1?k:i-k]);a(i/2+j+k)=ad(a(j+k),mod-t);a(j+k)=ad(a(j+k),t);}}}if(on==-1){for(i=0;i<N;i++)a(i)=mul(a(i),iN);}#undef a
}
int n,N;
void dft(int*a,int on){int i,j;for(j=0;j<n;j++){for(i=0;i<N;i++){if(!(i>>j*2&3))dft4(a+i,j*2,on);}}
}
int a[10010];
int c[1048576],f[1048576];
int t[4];
int main(){int m,i,j,x,res;scanf("%d%d",&n,&m);N=1<<n*2;for(i=1;i<=n;i++){for(j=1;j<=m;j++){scanf("%d",&x);a[j]=a[j]<<2|x;}}for(i=1;i<=m;i++)c[a[i]]++;reverse(c,c+N);for(i=0;i<N;i++){t[0]=t[1]=t[2]=t[3]=0;for(j=0;j<n;j++)t[i>>j*2&3]++;f[i]=min(min(t[1]+2*t[2]+3*t[3],t[0]+2*t[1]+3*t[2]),min(t[3]+2*t[0]+3*t[1],t[2]+2*t[3]+3*t[0]));}pre();dft(c,1);dft(f,1);for(i=0;i<N;i++)c[i]=mul(c[i],f[i]);dft(c,-1);res=inf;for(i=0;i<N;i++)res=min(res,c[i]);printf("%d",res);
}
