一个数学编程题目

题目：
　　有十个开关等间距排成一线，每个开关对应其上方的一盏灯(十盏灯也排成一线)。每按动一下开关，可以使对应的灯改变状态(原来亮着的将熄灭，原来熄灭的将被点亮)。

　　但是，由于开关之间的距离很小，每次按动开关时，相邻的一个开关也将被按动。例如：按动第5个开关，则实际上第4、5、6个开关都被按动。而按动靠边的第1个开关时，第1、2个开关都被按动。并且，无法只按动最靠边的一个开关。

　　现在给出十盏灯的初始的状态和目标状态，要求计算：从初始状态改变到目标状态所需要的最少操作次数。

　　函数接口：

　　int MinChange(const int Start[],const int End[]);

　　其中：Start表示了初始状态，End表示了目标状态。表示状态的数组(Start和End)中，若某元素为0表示对应的灯亮着，否则表示对应的灯没有亮。调用函数时保证Start和End数组长度均为10，并保证有解。

　　看了很多人的解法都是用循环遍历来判断是否达到最后要求，但是如果和线形代数结合的话，就有一种很快速的解法。

　　约定：以下所用的‘+’号都是‘异或’的运算。

　　先简化一下，假设有四个灯，初始状态s0～s3，目标状态是e0～e3，转换一次状态就是和1进行异或运算一次，所以状态转移矩阵为：

　　(s0,s1,s2,s3)+k0*(1,1,0,0)+k1*(1,1,1,0)+k2*(0,1,1,1)+k3*(0,0,1,1)=(e0,e1,e2,e3);

　　其中k(n)表示第n个开关所翻动的次数。并且，注意异或运算中a+b+b=a，所以，某个开关翻动偶数次的效果相当于没有翻动，翻动奇数次的效果相当于翻动一次;又由于异或运算满足交换律，所以翻动的顺序没有影响。综上每个开关翻动的次数只有1次或0次就足够了。

　　设m(n)=s(n)+e(n)，注意异或运算中的'-'也就是'+'，所以解线性方程组：

　　k0+k1 =m1;

　　k0+k1+k2 =m2;

　　k1+k2+k3=m3;

　　k2+k3=m4;

　　假设解存在，就可以算出通解(k0,k1,k2,k3)，再统计出通解中1的个数，就是所需要翻动的次数了。并且还可以知道哪些开关需要拨动，比如算出解是(1,0,1,0)就是第0和2个开关需要拨动一次。

　　因此针对本题目的10个灯泡，本人已算出这10元线性方程组的通解：

　　k0=m0+m2+m3+m5+m6+m8+m9;

　　k1=m2+m3+m5+m6+m8+m9;

　　k2=m0+m1;

　　k3=m3+m0+m1+m5+m6+m8+m9;

　　k4=m5+m6+m8+m9;

　　k5=m4+m3+m0+m1;

　　k6=m6+m4+m3+m0+m1+m8+m9;

　　k7=m8+m9;

　　k8=m7+m6+m4+m3+m0+m1;

　　k9=m9+m7+m6+m4+m3+m0+m1;

　　和上面一样，m(n)为开始状态与目标状态的每位异或。至于是否存在解，本人已将次系数矩阵化简为对角矩阵，可以看到系数矩阵的秩(Rank)与未知数的个数相等，所以无论是任何的输入和输出变换都能找到唯一解。

　　本人程序如下：

　　int MinChange(const int Start[],const int End[]){

　　int m[10];

　　int k[10];

　　int res=0;

　　int i,j,n;

　　for(i=0;i<10;i++){

　　m=Start^End; /* m[]=Start[] XOR End[] */

　　}

　　/* calculate roots */

　　k[0]=m[0]^m[2]^m[3]^m[5]^m[6]^m[8]^m[9];

　　k[1]=m[2]^m[3]^m[5]^m[6]^m[8]^m[9];

　　k[2]=m[0]^m[1];

　　k[3]=m[3]^m[0]^m[1]^m[5]^m[6]^m[8]^m[9];

　　k[4]=m[5]^m[6]^m[8]^m[9];

　　k[5]=m[4]^m[3]^m[0]^m[1];

　　k[6]=m[6]^m[4]^m[3]^m[0]^m[1]^m[8]^m[9];

　　k[7]=m[8]^m[9];

　　k[8]=m[7]^m[6]^m[4]^m[3]^m[0]^m[1];

　　k[9]=m[9]^m[7]^m[6]^m[4]^m[3]^m[0]^m[1];

　　/* count for switch times */

　　for(j=0;j<10;j++){

　　if(k[j]) res++;

　　}

　　/* display k(n); */

　　for(n=0;n<10;n++) printf("%d,",k[n]);

　　return res;

　　}

　　测试主程序：

　　main(){

　　int A[10]={1,1,1,0,0,1,0,1,1,0};

　　int B[10]={0,0,1,1,0,0,1,1,1,1};

　　int C;

　　C=MinChange(A,B);

　　printf("**%d**",C);

　　}

　　显示结果为：

　　1，0，0，0，1，1，1，1，0，1，**6**

　　就是如果要把状态A转为状态B，需要把第0，4，5，6，7，9号开关翻动一次，共6次。

　　测试验证结果正确：)

　　当然，此做法也有一个缺点，就是当灯的个数改变时，就要重新设定线形方程组的特解形式。