#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
#include <iostream>
using namespace std;

#include "linear.h"

int Linear::recv(int *p)
{
  newp(p);
  return uselast();
}

void Linear::init()
{
  N=0;
  rank=0;
}

int Linear::uselast()
{
  cleanlast();
  if(rank<K)
    if(M[N-1][rank])
    {
      swap(rank,N-1);
      rank++;
      cleanall();
      gaussian();
    }
  return rank;
}

void Linear::gaussian()
{
  int i;
  if(!(rank<K))
    return;
  for(i=rank;i<N;i++)
    if(M[i][rank])
      break;
  if(i>=N)
    return;
  swap(rank,i);
  rank++;
  cleanall();
  gaussian();
}

void Linear::cleanall()
{
  for(int i=rank;i<N;i++)
    if(M[i][rank-1])
      doxor(i,rank-1);
}

void Linear::cleanlast()
{
  for(int i=0;i<rank;i++)
    if(M[N-1][i])
      doxor(N-1,i);
}

void Linear::doxor(int i, int j)
{
  for(int k=0;k<K;k++)
    M[i][k]^=M[j][k];
}

void Linear::swap(int i, int j)
{
  int *row;
  row=M[i];
  M[i]=M[j];
  M[j]=row;
}

void Linear::pr()
{
  printf("M[%d][%d]=\r\n",N,K);
  for(int i=0;i<N;i++)
  {
    printf("%3d:",i);
    for(int j=0;j<K;j++)
      printf(" %s",M[i][j]?"X":"-");
    printf("\r\n");
  }
}

void Linear::newp(int *p)
{
  N++;
  for(int i=0;i<K;i++)
    M[N-1][i]=p[i];
}

Linear::Linear(int K_, int maxN_)
{
  K=K_;
  maxN=maxN_;
  N=0;
  M=new (int*)[maxN];
  for(int i=0;i<maxN;i++)
    M[i]=new int[K];
  rank=0;
}

Linear::~Linear()
{
  for(int i=0;i<maxN;i++)
    delete[] M[i];
  delete[] M;
}