转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove
惨痛的记忆,HDU第4次热身赛的某题,很快有了算法,完成代码。比赛期间花了3个小时在这题调试,晚上又花了两个多小时对照官方数据调试。终于找到了问题。
还是经典的开关问题,受最近的影响,直接高斯消元解异或方程组,由于需要最少次数,枚举自由变元即可,之前做的POJ 3185,类似的题目。可是POJ的数据弱了,当时没有暴露出问题,结果在今天的比赛中坑爹了,无限的WA让我十分郁闷。
之前的做法是找出自由变元的个数,对矩阵化简,然后搜索枚举自由变元,结果发现自由变元所在列并不在最后几列,这使得枚举自由变元的时候出现了问题,囧,每次把当前全为0的列后移到最后,也就是自由变元。之后再枚举
/*
ID:cxlove
PROB:poj 4200
DATA:2012.4.3
HINT:高斯消元法
*/
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,d;
int a[105][105],cnt;
int ta[105][105],x[105],ans[105],var0;
void debug(){
for(int i=1;i<=n;i++){
for(int j=1;j<=n+1;j++)
printf("%d ",a[i][j]);
printf("\n");
}
}
void dfs(int v){
if(v>n){
int temp=0;
for(int i=1;i<=n;i++)
x[i]=ans[i];
for(int i=1;i<=n;i++)
for(int j=1;j<=n+1;j++)
ta[i][j]=a[i][j];
for(int i=var0-1;i>=1;i--){
for(int j=i+1;j<=n;j++)
ta[i][n+1]^=(x[j]&&ta[i][j]);
x[i]=ta[i][n+1];
}
for(int i=1;i<=n;i++)
if(x[i])
temp++;
cnt=min(temp,cnt);
}
ans[v]=0;
dfs(v+1);
ans[v]=1;
dfs(v+1);
}
void gauss(){
int i,j,jj;
for(i=1,j=1,jj=1;i<=n&&jj<=n;jj++){
int k=i;
for(;k<=n;k++)
if(a[k][j])
break;
if(a[k][j]){
for(int r=1;r<=n+1;r++)
swap(a[i][r],a[k][r]);
for(k=1;k<=n;k++)
if(k!=i&&a[k][j]){
for(int r=1;r<=n+1;r++)
a[k][r]^=a[i][r];
}
i++; j++;
}
else{
int val=n-jj+j;
for(int cow=1;cow<=n;cow++)
swap(a[cow][j],a[cow][val]);
}
}
if(i==n+1){
cnt=0;
for(int r=1;r<=n;r++)
if(a[r][n+1])
cnt++;
printf("%d\n",cnt);
return ;
}
if(i<=n){
for(int r=i;r<=n;r++)
if(a[r][n+1]){
printf("impossible\n");
return ;
}
}
var0=i;
cnt=1<<30;
dfs(i);
printf("%d\n",cnt);
}
int main(){
int t;
// freopen("D.in","r",stdin);
// freopen("test.out","w",stdout);
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&d);
memset(a,0,sizeof(a));
for(int i=1;i<=n;i++)
scanf("%d",&a[i][n+1]);
for(int i=1;i<=n;i++)
for(int j=max(1,i-d);j<=min(n,i+d);j++)
a[i][j]=1;
gauss();
}
return 0;
}