暑期实践日志（四）

标签：

Description

There is a square wall which is made of n*n small square bricks. Some bricks are white while some bricks are yellow. Bob is a painter and he wants to paint all the bricks yellow. But there is something wrong with Bob‘s brush. Once he uses this brush to paint brick (i, j), the bricks at (i, j), (i-1, j), (i+1, j), (i, j-1) and (i, j+1) all change their color. Your task is to find the minimum number of bricks Bob should paint in order to make all the bricks yellow. 

Input

Output

Sample Input

2
3
yyy
yyy
yyy
5
wwwww
wwwww
wwwww
wwwww
wwwww

Sample Output

0
15
这道题目也是高斯消元法的经典题目，但是和之前的题目想比有些不同。
之前的开关问题在经过高斯消元后，只需要统计其中自有变量的个数，在进行移位运算即可；
在这道题中，需要解决的是在高斯消元后，若只有唯一解则，统计变化的旗子个数，若有多个解即有自由变量
那么要枚举自由变量的所有可能，取最小值。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std;
const int maxn=300;
int a[maxn][maxn];
int x[maxn];
int free_x[maxn];
int free_num;
int equ,var;
int gauss()
{
    int i,j,k,col,max_r;
    free_num=0;
    for(k=0,col=0;k<equ&&col<var;k++,col++)
    {
        max_r=k;
        for(int i=k+1;i<equ;i++)
        {
            if(abs(a[i][col])>abs(a[max_r][col]))
                max_r=i;
        }
        if(a[max_r][col]==0)
        {
            k--;
            free_x[free_num++]=col;
            continue;
        }
        if(max_r!=k)
        {
            for(int j=col;j<var+1;j++)
                swap(a[k][j],a[max_r][j]);
        }
        for(int i=k+1;i<equ;i++)
        {
            if(a[i][col]!=0)
            {
                for(int j=col;j<var+1;j++)
                    a[i][j]^=a[k][j];
            }
        }
    }
    for(int i=k;i<equ;i++)
    {
        if(a[i][col]!=0)
            return -1;
    }
    if(k<var)
        return var-k;
    for(int i=var-1;i>=0;i--)
    {
        x[i]=a[i][var];
        for(int j=i+1;j<var;j++)
            x[i]^=(a[i][j]&&x[j]);
    }
    return 0;
}
void init(int r)
{
    memset(a,0,sizeof(a));
    memset(x,0,sizeof(x));
    for(int i=0;i<r;i++)
    {
        for(int j=0;j<r;j++)
        {
            int t=i*r+j;
            a[t][t]=1;
            if(i>0)
                a[(i-1)*r+j][t]=1;
            if(i<r-1)
                a[(i+1)*r+j][t]=1;
            if(j>0)
                a[i*r+j-1][t]=1;
            if(j<r-1)
                a[i*r+j+1][t]=1;
        }
    }
}
void solve()
{
    int f=gauss();
    if(f==-1)
    {
        printf("inf\n");
        return;
    }
    else if(f==0)
    {
        int ans=0;
        for(int i=0;i<equ;i++)
            ans+=x[i];
        printf("%d\n",ans);
        return;
    }
    else
    {
        int ans=0x3f3f3f3f;
        int tot=(1<<f);
        for(int i=0;i<tot;i++)
        {
            int cnt=0;
            for(int j=0;j<f;j++)
            {
                if(i&(1<<j))
                {
                    x[free_x[j]]=1;
                    cnt++;
                }
                else
                    x[free_x[j]]=0;
            }
            for(int j=var-f-1;j>=0;j--)
            {
                int idx;
                for(idx=j;idx<var;idx++)
                    if(a[j][idx])
                        break;
                x[idx]=a[j][var];
                for(int l=idx+1;l<var;l++)
                    if(a[j][l])
                        x[idx]^=x[l];
                cnt+=x[idx];
            }
            ans=min(ans,cnt);
        }
        printf("%d\n",ans);
    }
}
char p[30][30];
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        memset(p,‘\0‘,sizeof(p));
        int n;
        scanf("%d",&n);
        getchar();
        equ=n*n;
        var=n*n;
        init(n);
        for(int i=0;i<n;i++)
        {
            scanf("%s",p[i]);
            for(int j=0;j<n;j++)
            {
                if(p[i][j]==‘y‘)
                    a[i*n+j][n*n]=0;
                else
                    a[i*n+j][n*n]=1;
            }
        }
        solve();
    }
    return 0;
}

暑期实践日志（四）

标签：

原文地址：http://www.cnblogs.com/hbutACMER/p/4652741.html