1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93
| #include<cstdio> #include<cstring> #include<algorithm> #include<iostream> #include<queue> #include<map> #include<stack> #include<cmath> #include<set> #include<bitset> #include<complex> #include<cstdlib> #include<assert.h> #define inc(i,l,r) for(int i=l;i<=r;i++) #define dec(i,l,r) for(int i=l;i>=r;i--) #define link(x) for(edge *j=h[x];j;j=j->next) #define mem(a) memset(a,0,sizeof(a)) #define ll long long #define eps 1e-8 #define succ(x) (1<<x) #define lowbit(x) (x&(-x)) #define mid (x+y>>1) #define sqr(x) ((x)*(x)) #define NM 1100005 #define nm 400005 using namespace std; const double pi=acos(-1); const ll inf=998244353; ll read(){ ll x=0,f=1;char ch=getchar(); while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();} while(isdigit(ch))x=x*10+ch-'0',ch=getchar(); return f*x; }
int n,m,tot; ll a[NM],b[NM],c[NM],ans; char _s[NM];
void fwt(ll*a,int f=0){ for(int k=1;k<n;k<<=1) for(int i=0;i<n;i+=k<<1) for(int j=0;j<k;j++){ ll x=a[i+j],y=a[i+j+k]; a[i+j]=x+y; a[i+j+k]=x-y; } if(f)inc(i,0,n-1)a[i]/=n; }
int main(){ n=read();m=read(); inc(i,1,n){ scanf("%s",_s+1); inc(j,1,m)if(_s[j]-'0')a[j]|=succ(i-1); } inc(i,1,m)c[a[i]]++; tot=succ(n); inc(i,0,tot-1)b[i]=__builtin_popcount(i),b[i]=min(b[i],n-b[i]); n=tot; fwt(c);fwt(b); inc(i,0,n-1)b[i]=b[i]*c[i]; fwt(b,1); ans=inf; inc(i,0,n-1)ans=min(ans,b[i]); printf("%lld\n",ans); return 0; }
|