爪哇- 522 434 228 213个字符
通过直接乘法系统地检查所有可能的整数n元组,直到找到一个可行的整数,然后进行求解。
函数将增强矩阵A,试验解向量x和维度n作为输入-输出解向量x。请注意,向量x实际上比维度大一,以帮助逐步解决可能的问题。(如果我将变量A,x,n,j,k,s声明为实例变量,则该函数将短31个字符-共计182个字符,但这感觉像是将规则弯曲得太远了。)
int[]Z(int[][]A,int[]x,int n){int j,k,s;for(;;){for(j=0;j<n;j++){for(k=s=0;k<n;s+=A[j][k]*x[k++]);if(s!=A[j][n])j+=n;}if(j==n)return x;for(j=0;j<=n;j++)if(x[j]!=x[n]||j==n){x[j]++;for(k=0;k<j;x[k++]=-x[n]);j=n;}}}
测试程序(有点不实用):
import java.util.*;
class MatrixSolver{
public MatrixSolver() {
Scanner p=new Scanner(System.in); //initialize everything from stdin
int j,k,n=p.nextInt(),A[][]=new int[n][n+1],x[]=new int[n+1];
for(j=0;j<n;j++)for(k=0;k<=n;A[j][k++]=p.nextInt());
x=Z(A,x,n); //call the magic function
for(j=0;j<n;j++) System.out.print(x[j]+" "); //print the output
}
public static void main(String[]args){
new MatrixSolver();
}
int[]Z(int[][]A,int[]x,int n){
int j,k,s;
for(;;){
for(j=0;j<n;j++){ //multiply each row of matrix by trial solution and check to see if it is correct
for(k=s=0;k<n;s+=A[j][k]*x[k++]);
if(s!=A[j][n])j+=n;
}
if(j==n)return x; //if it is correct return the trial solution
for(j=0;j<=n;j++)if(x[j]!=x[n]||j==n){//calculate the next trial solution
x[j]++;
for(k=0;k<j;x[k++]=-x[n]);
j=n;
}
}
}
}
程序将来自stdin的输入作为空格分隔的整数,如下所示:首先,问题的维数,其次,是逐行扩展矩阵的条目。
样品运行:
$java -jar MatrixSolver.jar
3 2 1 -1 8 -3 -1 2 -11 -2 1 2 -3
2 3 -1
我遵循Victor关于循环和“公共”的建议,剃掉了几个角色,将RHS而不是单独存储在增强矩阵中,并在我的试用解决方案中添加了一个额外的条目以简化每个新试用解决方案的生成。OP还表示,一个功能就足够了-无需计算整个程序。