#include
#include
#include
//const?int?LENGTH=5000;
void?main（）
{??
???ifstream?fin（“putin.txt“）;
???ofstream?fout（“putout.txt“）;
???int?n_interiorm_interiorm_iijkiter_k=0;
???float?（*a_matrix）[5000];
???a_matrix=new?float[5000][5000];
???float?*d_matrix;
???d_matrix=new?float[5000];
???float?*c_matrix;
???c_matrix=new?float[5000];
???float?*b_matrix;
???b_matrix=new?float[5000];
???float?*x0_matrix;
???x0_matrix=new?float[5000];
???float?*y0_matrix;
???y0_matrix=new?float[5000];
???float?*z0_matrix;
???z0_matrix=new?float[5000];
???float?*ep_node;
???ep_node=new?float[5000];
???float?*ep_node1;
???ep_node1=new?float[5000];
???float?*x_node;
???x_node=new?float[5000];
???float?*z_node;
???z_node=new?float[5000];
???float?*dx_node;
???dx_node=new?float[5000];
???float?*uu_node;
???uu_node=new?float[5000];
???float?*rd_node;
???rd_node=new?float[5000];
???float?*rp_node;
???rp_node=new?float[5000];
???float?（*ada_matrix）[5000];
???ada_matrix=new?float[5000][5000];
???float?*bda_matrix;
???bda_matrix=new?float[5000];


???float?*BB;
???BB=new?float[5000];
???float?*U;
???U=new?float[500000];
???int?*JU;
???JU=new?int[500000];
???int?*MD;
???MD=new?int[5000];
???float?*D;
???D=new?float[5000];
???float?*DJ;
???DJ=new?float[5000];
???

??//?float?a_matrix[5000][5000]d_matrix[5000]c_matrix[5000]b_matrix[5000];
??//???float?x0_matrix[5000]y0_matrix[5000]z0_matrix[5000];
??//???float?ep_node[5000]ep_node1[5000]x_node[5000]z_node[5000];
??//???float?dx_node[5000]uu_node[5000]rd_node[5000]rp_node[5000];
??//???float?m_iada_matrix[5000][5000]bda_matrix[5000];
???float?primaldualu_node;
???float?gap=0.00;
???float?ruta;
?????ruta=0.01;///非常重要的中心參數
//???int?BB[5000]U[500000]JU[500000]MD[5000]D[5000]DJ[5000];
???/////讀入數據
??

???fin>>n_interior>>m_interior;
?????????for（i=1;i<=m_interior;i++）
??????	???for（j=1;j<=n_interior;j++）
??????		???fin>>a_matrix[i][j];
?????????for（i=1;i<=n_interior;i++）
??????	???fin>>c_matrix[i];
?????????for（i=1;i<=m_interior;i++）
??????	???fin>>b_matrix[i];
?????????
?????????//////////////賦初值
?????????for（i=1;i<=n_interior;i++）
?????????{??x0_matrix[i]=1;
????????????z0_matrix[i]=1;
?????????}
?????????for（i=1;i<=m_interior;i++）
??????	???y0_matrix[i]=0.0;
?????????for（i=1;i<=n_interior;i++）
??????	???gap=gap+z0_matrix[i]*x0_matrix[i];
?????????u_node=float（fabs（gap））/n_interior;?//對數障礙因子
?????????///step?1
?????????aa:
??????	??for（i=1;i<=m_interior;i++）??//計算A1-6
?????????{??rp_node[i]=b_matrix[i];
????????????for（j=1;j<=n_interior;j++）
??????		??if?（a_matrix[i][j]!=0.0）?
??????			??rp_node[i]=rp_node[i]-a_matrix[i][j]*x0_matrix[j];
?????????}
?????????for（i=1;i<=n_interior;i++）?????????//計算A1-7
?????????{??rd_node[i]=c_matrix[i]-z0_matrix[i];
????????????for（j=1;j<=m_interior;j++）
??????		??if?（a_matrix[j][i]!=0.0）?
??????			??rd_node[i]=rd_node[i]-a_matrix[j][i]*y0_matr

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-04-05?08:07??內點法機組組合\
?????目錄???????????0??2011-04-05?08:07??內點法機組組合\Debug\
?????文件??????258145??2008-09-04?15:27??內點法機組組合\Debug\interior?method.exe
?????文件???????21276??2008-09-04?15:27??內點法機組組合\Debug\interior?method.obj
?????文件??????590848??2008-09-04?15:27??內點法機組組合\Debug\interior?method.pdb
?????文件???????26095??2008-05-23?21:44??內點法機組組合\Debug\interior?method.sbr
?????文件???????61440??2008-09-04?15:27??內點法機組組合\Debug\vc60.pdb
?????文件???????10143??2008-09-04?15:27??內點法機組組合\interior?method.cpp
?????文件????????4388??2008-03-23?20:42??內點法機組組合\interior?method.dsp
?????文件?????????555??2008-03-23?20:08??內點法機組組合\interior?method.dsw
?????文件??????237648??2008-05-23?18:08??內點法機組組合\interior?method.exe
?????文件???????50176??2010-01-28?12:07??內點法機組組合\interior?method.ncb
?????文件???????20652??2008-05-23?18:08??內點法機組組合\interior?method.obj
?????文件???????69632??2010-01-28?12:07??內點法機組組合\interior?method.opt
?????文件??????574464??2008-05-23?18:08??內點法機組組合\interior?method.pdb
?????文件????????1816??2008-09-04?15:27??內點法機組組合\interior?method.plg
?????文件????????1543??2008-11-06?14:30??內點法機組組合\putin.txt
?????文件?????2380366??2008-08-29?10:09??內點法機組組合\putin000.txt
?????文件?????????115??2010-01-28?12:04??內點法機組組合\putout.txt
?????文件???????61440??2008-05-23?18:08??內點法機組組合\vc60.pdb