M_4的一个算法

设a_1,a_2,a_3,a_4是正整数,(a_1,a_2,a_3,a_4)=1,线性型a_1x_1+a_2x_2+a_3x_3+a_4x_4 ,x_i≥0,i=1,2,3,4.不能表出的最大整数记为M_4.而线性型ax+by+cz+dw,x≥1,y≥1,z≥1,w≥1,不能表出的最大整数记为N_4~′.其中,a_1=a·(r_1r_3r_4),a_2=b·(r_1r_2r_4),a_3=c·(r_1r_2r_3),a_4=d·(r_2r_3r_4);r_1=(a_1,a_2,a_3),r_2=(a_2,a_3,a_4),r_3=(a_3,a_4,a_1),r_4=(a_4,a_1,a_2).通过范式组:ak_a=bx_a+cy_a+dz_a,x_a≥0,y_a≥0,z_a≥0,bk_b=cx_b+dy_b+az_b,x_b≥0,y_b≥0,z_b≥0,ck_c=dx_c+ay_c+bz_c,x_c≥0,y_c≥0,z_c≥0,dk_d=ax_d+by_d+cz_d,x_d≥0,y_d≥0,z_d≥0.算出N_4~′,则M_4=(r_1r_2r_3r_4)N_4~′-a_1-a_2-a_3-a_4.