hihoCoder 1196 高斯消元·二

Description

一个黑白网格,点一次会改变这个以及与其连通的其他方格的颜色,求最少点击次数使得所有全部变成黑色.

Sol

高斯消元解异或方程组.

先建立一个方程组.

\(x_i\) 表示这个点是否被用过.

因为第二次使用同一个点,这个点的贡献就被消除了,所以每个点只会被用 0/1 次.

\(a_{ij}\) 表示 \(j\) 点对 \(i\) 是否有影响,有影响为 1 否则为 0.

最后的一位表示最后的状态^最初的状态.

这样就列出来了 \(n*m\) 个方程组,一共 \(n*m\) 个未知数.

然后求解,就跟高斯消元一样,用异或的方式消掉就可以了.

Code

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std; #define H(i,j) ((i-1)*m+j)
const int N = 35; int n,m;
int a[N][N]; void print(){
cout<<"----------------------------"<<endl;
for(int i=1;i<=30;i++){
for(int j=1;j<=31;j++) cout<<a[i][j];
cout<<endl;
}
} void init(){
n=5,m=6;
for(int i=1;i<=n;i++){
char s[N];memset(s,0,sizeof(s));
scanf("%s",s+1);
for(int j=1;j<=m;j++) a[H(i,j)][n*m+1]=1^(s[j]-'0');
// cout<<a[H(i,j)][n*m+1];
// cout<<endl<<s+1<<endl<<endl;
} for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) {
a[H(i,j)][H(i,j)]=1;
if(i+1<=n) a[H(i,j)][H(i+1,j)]=1;
if(i-1>=1) a[H(i,j)][H(i-1,j)]=1;
if(j+1<=m) a[H(i,j)][H(i,j+1)]=1;
if(j-1>=1) a[H(i,j)][H(i,j-1)]=1;
}
// print();
} void gauss(int n){
for(int i=1,r,j;i<n;i++){
for(r=j=i;j<n;j++) if(a[j][i]) { r=j;break; }
// cout<<i<<" "<<r<<endl;
if(r!=i) for(int k=i;k<=n;k++) swap(a[i][k],a[r][k]);
for(int k=i+1;k<n;k++) if(a[k][i]) {
for(j=i;j<=n;j++) a[k][j]^=a[i][j];
}
// print();
}
// print();
for(int i=n-1;i;i--){
for(int j=i+1;j<n;j++) a[i][n]^=a[j][n]*a[i][j];
} } void out(){
int ans=0;
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(a[H(i,j)][n*m+1]) ans++;
cout<<ans<<endl;
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(a[H(i,j)][n*m+1]) cout<<i<<" "<<j<<endl;
} int main(){
// freopen("in.in","r",stdin);
// freopen("log.out","w",stdout);
init();
gauss(31);
out(); return 0;
}

  

上一篇:iOS UICollectionView(转一) XIB+纯代码创建:cell,头脚视图 cell间距


下一篇:TP5.x——打印SQL语句