当前位置:首页 > 资讯 > info5 > 正文

POJ 3185 The Water Bowls(高斯消元)

发表于: 2015-01-27   作者:u013013910   来源:转载   浏览:
摘要: 题目地址:POJ3185     醉了。,。高斯消元这种题样例都不过的时候怎么调试好。。。干瞪着看了半天代码终于发现是一个符合写错了。。    这题同样是高斯消元+自由元枚举。没什么好说的。代码如下:#include #include #include #include #include #include #include #include #include usingnamespacestd;

题目地址:POJ 3185

         醉了。,。高斯消元这种题样例都不过的时候怎么调试好。。。干瞪着看了半天代码终于发现是一个符合写错了。。

        这题同样是高斯消元+自由元枚举。没什么好说的。

代码如下:

#include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL __int64
#define pi acos(-1.0)
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const double eqs=1e-6;
int b[31], a[30][30];
int free_num, free_x[30], state[30];
int gauss()
{
        int i, j, k, h, max_r, tmp;
        free_num=0;
        for(i=0,j=0;i<20&&j<20;i++,j++){
                max_r=i;
                for(k=i+1;k<20;k++){
                        if(a[k][j]>a[max_r][j]) max_r=k;
                }
                if(max_r!=i){
                        for(k=j;k<=20;k++){
                                swap(a[i][k],a[max_r][k]);
                        }
                }
                if(!a[i][j]){
                        free_x[free_num++]=j;
                        i--;
                        continue ;
                }
                for(k=i+1;k<20;k++){
                        if(a[k][j]==0) continue ;
                        for(h=j;h<21;h++){
                                a[k][h]^=a[i][h];
                        }
                }
        }
        //printf("%d\n",free_num);
        tmp=i;
        for(i=tmp;i<20;i++){
                if(a[i][20])
                        return INF;
        }
        int tot=1<<free_num, cnt, ans=INF;
        //printf("%d\n",tot);
        for(i=0;i<tot;i++){
                cnt=0;
                for(j=0;j<free_num;j++){
                        if(i&(1<<j)){
                                state[free_x[j]]=1;
                                cnt++;
                        }
                        else{
                                state[free_x[j]]=0;
                        }
                }
                for(j=tmp-1;j>=0;j--){
                        int t=0;
                        while(a[j][t]==0) t++;
                        state[t]=a[j][20];
                        for(k=t+1;k<20;k++){
                                if(a[j][k]) state[t]^=state[k];
                        }
                        cnt+=state[t];
                }
                //printf("--%d\n",cnt);
                ans=min(ans,cnt);
        }
        return ans;
}
int main()
{
        int i;
        for(i=0;i<20;i++){
                scanf("%d",&b[i]);
        }
        memset(a,0,sizeof(a));
        for(i=0;i<20;i++){
                a[i][i]=1;
                if(i) a[i][i-1]=1;
                if(i!=19) a[i][i+1]=1;
                a[i][20]=b[i];
        }
        printf("%d\n",gauss());
        return 0;
}


POJ 3185 The Water Bowls(高斯消元)

  • 0

    开心

    开心

  • 0

    板砖

    板砖

  • 0

    感动

    感动

  • 0

    有用

    有用

  • 0

    疑问

    疑问

  • 0

    难过

    难过

  • 0

    无聊

    无聊

  • 0

    震惊

    震惊

编辑推荐
EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 8353 Accepted:
高斯消元法(理论+习题) 高斯消元主要用来解多元一次线性方程组,主要用到的线性代数的知识详见矩
[转:http://www.cnblogs.com/pegasus/archive/2011/07/31/2123195.html] 高斯消元法(Gaussian eli
Flip Game Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 37055 Accepted: 16125 De
Painter's Problem Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 5352 Accepted: 2
Language: Default EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissio
EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 6608 Accepted:
EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 7817 Accepted:
【POJ 1681】 Painter's Problem (高斯消元) Painter's Problem Time Limit: 1000MS Memory Limit
EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 8364 Accepted:
版权所有 IT知识库 CopyRight © 2009-2015 IT知识库 IT610.com , All Rights Reserved. 京ICP备09083238号