请教伽罗华域矩阵求逆

2011-03-31 05:07:25 749 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

使用GDScript(godot)生成仿杀戮尖塔地图使用GDScript(godot)生成仿杀戮尖塔地图 extends Node2D export (PackedScene) var Ball var 被连接数 = {} # 每个怪物被下层怪物连接的数量最多为2 var 连接1 = {} # 第一次连接的上层怪物位置 var 连接2 = {} # 第二次连接的上层怪物位置 var 怪物数量 = [0,0,0,0,0,0,0,0,0,0,0,0,0,0] # 每层怪物的数量 var 总层数 =8 # 总层数 var 怪物坐标行 = {} # 阅读详情

最近在学Raptor码，里面有一个矩阵求逆的过程，是在伽罗华域里进行的，做纠错码的朋友肯定很熟悉，伽罗华域的元操作是异或，我将清华大学出版《C语言常用算法程序集》第二版这本书里的矩阵求逆改了一下，这本书有随书的源代码，我改为如下：



int brinv(int **a,int n)  //a对应矩阵二维矩阵A，n是矩阵的维数L ，计算A的逆

{

 int *is,*js,i,j,k;

 int d,p;

 is=malloc(n*sizeof(int));

 js=malloc(n*sizeof(int));

 for (k=0; k<=n-1; k++)

 { 

  d=0;

     for (i=k; i<=n-1; i++)

  for (j=k; j<=n-1; j++)

  {

   p=a[i][j];

        if (p>d) { d=p; is[k]=i; js[k]=j;}

     }

    if (d+1==1)

     { 

     free(is); free(js); printf("err**not inv\n");

        return(1);

     }

 if (is[k]!=k)

  for (j=0; j<=n-1; j++)

  { 

       p=a[k][j]; a[k][j]=a[is[k]][j]; a[is[k]][j]=p;

     }

 if (js[k]!=k)

  for (i=0; i<=n-1; i++)

  { 

         p=a[i][k]; a[i][k]=a[i][js[k]]; a[i][js[k]]=p;

     }

 a[k][k]=a[k][k];

 for (j=0; j<=n-1; j++)

  if (j!=k)

  {  a[k][j]=a[k][j]*a[k][k];}

  for (i=0; i<=n-1; i++)

   if (i!=k)

    for (j=0; j<=n-1; j++)

     if (j!=k)

     { 

                a[i][j]=a[i][j]^(a[i][k]*a[k][j]);

              }

  for (i=0; i<=n-1; i++)

   if (i!=k)

   {  a[i][k]=a[i][k]*a[k][k];}

 }

 for (k=n-1; k>=0; k--)

 { if (js[k]!=k)

 for (j=0; j<=n-1; j++)

 {

    p=a[k][j]; a[k][j]=a[js[k]][j]; a[js[k]][j]=p;

 }

 if (is[k]!=k)

  for (i=0; i<=n-1; i++)

  { 

       p=a[i][k]; a[i][k]=a[i][is[k]]; a[i][is[k]]=p;

     }

 }

 free(is); free(js);

 return 0;

}