hihocoder 1195&1196

1195
1196

存一发高斯消元模板。（填坑）
1195是纯模板，1196是解异或方程组。
代码中有注释。

【代码】

• 1195

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <queue>
#include <cmath>
#define N 505
#define M 20205
#define INF 0x7fffffff
using namespace std;
typedef long long ll;
const double EPS=1e-10;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m;
double s[N][N];
double a[N][N];
double ans[N];

int Gauss_Elimination()
{
    int k=1,col=1;
    for(;col<=n&&k<=m;col++,k++)
    {
        int r=k;
        for(int i=k+1;i<=m;i++) if(fabs(a[i][col])>fabs(a[r][col])) r=i;
        if(fabs(a[r][col])<EPS){k--;continue;}//自由元
        if(r!=k) swap(a[r],a[k]);
        for(int i=k+1;i<=m;i++)
        {
            if(fabs(a[i][col])>EPS)
            {
                double t=a[i][col]/a[k][col];
                for(int j=col;j<=n+1;j++) a[i][j]-=a[k][j]*t;
                //当前行为k,当前列为col,要化成上三角矩阵，所以这里将后面的式子中当前列全部变为0，做法就是将第k行化为第col列与当前行系数相同的式子，然后当前行依次减掉它。
            }
        }
    }
    for(int i=k;i<=m;i++)
        if(fabs(a[i][n+1])>EPS) return -1;//无解
    if(k<=n) return n-k+1;//返回自由元个数（有多解）
    for(int i=n;i;i--)
    {
        for(int j=i+1;j<=n;j++)
            a[i][n+1]-=a[i][j]*a[j][n+1];
        a[i][n+1]/=a[i][i];
    }//回带过程
    return 0;
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=m;i++) 
    for(int j=1;j<=n+1;j++) a[i][j]=read();
    int t=Gauss_Elimination();
    if (t<0) puts("No solutions");
    else if (t>0) puts("Many solutions");
    else 
        for(int i=1;i<=n;i++) printf("%d\n",int(a[i][n+1]+0.5));
    return 0;
}

• 1196

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <queue>
#include <cmath>
#define N 505
#define M 20205
#define INF 0x7fffffff
using namespace std;
typedef long long ll;
const double EPS=1e-10;
ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n;
int s[7][7],a[35][35]; 

int id(int x,int y){
    return (x-1)*6+y;
}

void Get(int X){
    int x=X%6;if(!x) x=6;
    int y=(X-x)/6+1;
    printf("%d %d\n",y,x);
}

void Gauss_Elimination()
{
    int k=1,col=1;
    for(;k<=n&&col<=n;k++,col++)
    {
        int r=k;
        for(int i=k;i<=n;i++) if(a[i][col]==1) {
            r=i;break;
        }
        if(!a[r][col]) {k--;continue;}
        if(r!=k) swap(a[r],a[k]);
        for(int i=k+1;i<=n;i++)
        {
            if(!a[i][col]) continue;
            for(int j=col;j<=n+1;j++) a[i][j]^=a[k][j];
            //和解线性方程组有些不同，但都是化成上三角形式，想想发现其实就是将两个式子合并，将当前列化为0的新等式。
        }
    }
    for(int i=n;i;i--)
    {
        for(int j=i+1;j<=n;j++)
            a[i][n+1]^=(a[j][n+1]*a[i][j]);//只有两个都为1才算。
        //这里并不需要管a[i][i]，输出看看发现全是1.
    }
    int sum=0;
    for(int i=1;i<=n;i++) if(a[i][n+1]) sum++;
    printf("%d\n",sum);
    for(int i=1;i<=n;i++) if(a[i][n+1]) Get(i);
}

int main()
{
    n=30;
    for(int i=1;i<=5;i++)
    {
        char ch[7];scanf("%s",ch);
        for(int j=1;j<=6;j++) 
        {
            s[i][j]=ch[j-1]-'0';
            a[id(i,j)][31]=s[i][j]^1;
        }
    }
    for(int i=1;i<=30;i++)
    {
        a[i][i]=1;
        if(i%6!=1) a[i-1][i]=1;
        if(i%6!=0) a[i+1][i]=1;
        if(i>6) a[i-6][i]=1;
        if(i<25) a[i+6][i]=1;
    }
    Gauss_Elimination();
    return 0;
}