<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AJOR</journal-id><journal-title-group><journal-title>American Journal of Operations Research</journal-title></journal-title-group><issn pub-type="epub">2160-8830</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajor.2016.61001</article-id><article-id pub-id-type="publisher-id">AJOR-62642</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Solving the Binary Linear Programming Model in Polynomial Time
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>lias</surname><given-names>Munapo</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>School of Accounting, Economics and Decision Sciences, North West University, Mafeking, South Africa</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>emunapo@gmail.com</email></corresp></author-notes><pub-date pub-type="epub"><day>11</day><month>01</month><year>2016</year></pub-date><volume>06</volume><issue>01</issue><fpage>1</fpage><lpage>7</lpage><history><date date-type="received"><day>26</day>	<month>August</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>5</month>	<year>January</year>	</date><date date-type="accepted"><day>11</day>	<month>January</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   The paper presents a technique for solving the binary linear programming model in polynomial time. The general binary linear programming problem is transformed into a convex quadratic programming problem. The convex quadratic programming problem is then solved by interior point algorithms. This settles one of the open problems of whether P = NP or not. The worst case complexity of interior point algorithms for the convex quadratic problem is polynomial. It can also be shown that every liner integer problem can be converted into binary linear problem. 
 
</p></abstract><kwd-group><kwd>NP-Complete</kwd><kwd> Binary Linear Programming</kwd><kwd> Convex Function</kwd><kwd> Convex Quadratic Programming Problem</kwd><kwd> Interior Point Algorithm and Polynomial Time</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The binary linear programming (BLP) model is NP-complete and up to now we have not been aware of any polynomial algorithm for this model. See for example Fortnow [<xref ref-type="bibr" rid="scirp.62642-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.62642-ref2">2</xref>] for more on complexity. In this paper we present a technique for transforming the BLP model into a convex quadratic programming (QP) problem. The optimal solution of the resultant convex QP is also the optimal solution of the original problem BLP. This solves one of the famous open problems of whether P = NP or not.</p></sec><sec id="s2"><title>2. The BLP Model</title><p>Let any BLP model be represented by</p><disp-formula id="scirp.62642-formula1"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x6.png"  xlink:type="simple"/></disp-formula><p>Any minimization BLP can be converted into maximization form and vice versa. There are several strategies for solving mixed 0 - 1 integer problems that are presented in Adams and Sherali [<xref ref-type="bibr" rid="scirp.62642-ref3">3</xref>] .</p></sec><sec id="s3"><title>3. Convex Quadratic Programming Model</title><p>Let a quadratic programming problem be represented by (2).</p><disp-formula id="scirp.62642-formula2"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x7.png"  xlink:type="simple"/></disp-formula><p>We assume that:</p><p>1) matrix Q is symmetric and positive definite,</p><p>2) function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x8.png" xlink:type="simple"/></inline-formula> is strictly convex,</p><p>3) since constraints are linear then the solution space is convex,</p><p>4) any maximization quadratic problem can be changed into a minimization and vice versa.</p><p>When the function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x9.png" xlink:type="simple"/></inline-formula> is strictly convex for all points in the convex region then the quadratic problem has a unique local minimum which is also the global minimum [<xref ref-type="bibr" rid="scirp.62642-ref4">4</xref>] -[<xref ref-type="bibr" rid="scirp.62642-ref6">6</xref>] .</p></sec><sec id="s4"><title>4. Transforming BLP into a Convex/Concave Quadratic Programming Problem</title><p>Our problem is to transform problem (1) into (2) and once that is done then (2) can be solved in polynomial time implying P = NP. Interior point algorithms can solve the convex/concave QP problem in polynomial time.</p><sec id="s4_1"><title>4.1. Rules with Binary Variables</title><p>Binary variables have certain special features that we can capitalize on when solving.</p><sec id="s4_1_1"><title>4.1.1. Rule 1</title><p>Given any binary variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x10.png" xlink:type="simple"/></inline-formula> then slack variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x11.png" xlink:type="simple"/></inline-formula> is also binary in the optimal solution.</p><p>Proof</p><disp-formula id="scirp.62642-formula3"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x12.png"  xlink:type="simple"/></disp-formula><p>Case 1: When <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x13.png" xlink:type="simple"/></inline-formula> then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x14.png" xlink:type="simple"/></inline-formula></p><p>Case 2: When <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x15.png" xlink:type="simple"/></inline-formula> then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x16.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s4_1_2"><title>4.1.2. Rule 2</title><p>For any binary variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x17.png" xlink:type="simple"/></inline-formula> and slack variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x18.png" xlink:type="simple"/></inline-formula> the following must hold at optimality for BLPs.</p><disp-formula id="scirp.62642-formula4"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x19.png"  xlink:type="simple"/></disp-formula><p>The proof is the same as the one given in 4.1.1. Note that it is only binary variables that can satisfy (4). Even though none binary values such as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x20.png" xlink:type="simple"/></inline-formula> then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x21.png" xlink:type="simple"/></inline-formula> can satisfy (3) the same values cannot satisfy (4), i.e., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x22.png" xlink:type="simple"/></inline-formula>The binary variable slack relationship given in (4) is the backbone of this paper.</p></sec></sec><sec id="s4_2"><title>4.2. Forcing Variables to Assume Binary Variables</title><p>The main weakness of the objective function given in (1) is that it does not force variables to assume binary values. In this paper we alleviate this challenge by adding a nonlinear extension to the objective function as given in (5).</p><disp-formula id="scirp.62642-formula5"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x23.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x24.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x25.png" xlink:type="simple"/></inline-formula> is a very large constant. The constant <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x26.png" xlink:type="simple"/></inline-formula> is very large in terms of its size compared to any of the coefficients in the objective function. This large value can be approximated as:</p><disp-formula id="scirp.62642-formula6"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x27.png"  xlink:type="simple"/></disp-formula><p>Proof</p><disp-formula id="scirp.62642-formula7"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x28.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62642-formula8"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x29.png"  xlink:type="simple"/></disp-formula><p>Since from Rule 2, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x30.png" xlink:type="simple"/></inline-formula>then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x31.png" xlink:type="simple"/></inline-formula> is minimized when</p><disp-formula id="scirp.62642-formula9"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x32.png"  xlink:type="simple"/></disp-formula><p>In other words <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x33.png" xlink:type="simple"/></inline-formula> is maximized when variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x34.png" xlink:type="simple"/></inline-formula> and slack variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x35.png" xlink:type="simple"/></inline-formula> are integers. In this paper the nonlinear extension <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x36.png" xlink:type="simple"/></inline-formula> is called an enforcer. An enforcer is a function, a set of constraint(s) or combination of both added to a problem to force an optimal solution with desired features such integrality.</p></sec><sec id="s4_3"><title>4.3. Convexity of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x37.png" xlink:type="simple"/></inline-formula></title><p>A function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x38.png" xlink:type="simple"/></inline-formula> is convex if and only if it has second-order partial derivatives for each point <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x39.png" xlink:type="simple"/></inline-formula> and for each <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x40.png" xlink:type="simple"/></inline-formula> all principal minors of the Hessian matrix are none negative.</p><p>Proof</p><p>In this case</p><disp-formula id="scirp.62642-formula10"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x41.png"  xlink:type="simple"/></disp-formula><p>This has continuous second order partial derivatives and the 2n by 2n Hessian matrix is given by</p><disp-formula id="scirp.62642-formula11"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x42.png"  xlink:type="simple"/></disp-formula><p>Since all principal minors of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x43.png" xlink:type="simple"/></inline-formula> are nonnegative then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x44.png" xlink:type="simple"/></inline-formula> is convex. See Winston [<xref ref-type="bibr" rid="scirp.62642-ref7">7</xref>] for more on convex functions.</p></sec><sec id="s4_4"><title>4.4. Convex Quadratic Programming Form</title><p>The function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x45.png" xlink:type="simple"/></inline-formula> can be expressed in the convex quadratic programming form</p><disp-formula id="scirp.62642-formula12"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x46.png"  xlink:type="simple"/></disp-formula><p>where matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x47.png" xlink:type="simple"/></inline-formula> is of dimension 2n by 2n, symmetric and positive definite as given in (13).</p><disp-formula id="scirp.62642-formula13"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x48.png"  xlink:type="simple"/></disp-formula><p>Thus matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x49.png" xlink:type="simple"/></inline-formula> is symmetric and positive definite. Note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x50.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s4_5"><title>4.5. Complexity of Convex Quadratic Programming</title><p>The main reason for converting a BLP into a convex quadratic programming model is to take advantage of the availability of interior point algorithms which can solve convex QPs in polynomial time [<xref ref-type="bibr" rid="scirp.62642-ref8">8</xref>] . If any BLP can be converted into a convex quadratic problem, then any BLP can be solved in polynomial time.</p></sec><sec id="s4_6"><title>4.6. Proof of Optimality</title><p>The proof is easily shown by reducing the convex quadratic objective function to the original linear form given in (1). The proposed objective function of the convex QP is reduced as follows:</p><disp-formula id="scirp.62642-formula14"><graphic  xlink:href="http://html.scirp.org/file/1-1040416x51.png"  xlink:type="simple"/></disp-formula><p>Since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x52.png" xlink:type="simple"/></inline-formula> then,</p><disp-formula id="scirp.62642-formula15"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x53.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62642-formula16"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x54.png"  xlink:type="simple"/></disp-formula><p>In other words <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x55.png" xlink:type="simple"/></inline-formula> is a constant and the objective function is the same as:</p><disp-formula id="scirp.62642-formula17"><graphic  xlink:href="http://html.scirp.org/file/1-1040416x56.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x57.png" xlink:type="simple"/></inline-formula> is binary for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x58.png" xlink:type="simple"/></inline-formula> this is the original form given in (1).</p></sec><sec id="s4_7"><title>4.7. Infeasible Binary Integer Solution Space</title><p>In this case the solution of the convex OP will not be integer. The objective,</p><disp-formula id="scirp.62642-formula18"><graphic  xlink:href="http://html.scirp.org/file/1-1040416x59.png"  xlink:type="simple"/></disp-formula><p>forces variables to binary or integral values if an integer point exists in the solution space. If an integer point does not exists in the solution space the large constant <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x60.png" xlink:type="simple"/></inline-formula> in the objective forces variables to assume values whose sum of squares are near one and not necessarily one. In other words the variables will assume values <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x61.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x62.png" xlink:type="simple"/></inline-formula> such that</p><disp-formula id="scirp.62642-formula19"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x63.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62642-formula20"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x64.png"  xlink:type="simple"/></disp-formula></sec><sec id="s4_8"><title>4.8. Mixed BLP Models</title><p>In some BLP problems that occur in real life, a fraction of some of the variables may not be restricted to integer values. In this case the enforcer <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x65.png" xlink:type="simple"/></inline-formula> is composed of only those variables that are supposed be binary and integer.</p></sec><sec id="s4_9"><title>4.9. Interior Point Algorithm for Convex QP</title><p>Any maximization BLP problem can be converted into a minimization BLP and vice versa. This can be done by the substitution given in (18).</p><disp-formula id="scirp.62642-formula21"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x66.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x67.png" xlink:type="simple"/></inline-formula> is also a binary variable.</p><p>Suppose the primal-dual pair of the convex QP is given by (19) and (20).</p><p>Primal:</p><disp-formula id="scirp.62642-formula22"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x68.png"  xlink:type="simple"/></disp-formula><p>Dual:</p><disp-formula id="scirp.62642-formula23"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x69.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x70.png" xlink:type="simple"/></inline-formula> is free, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x71.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x72.png" xlink:type="simple"/></inline-formula> is a diagonal matrix.</p><p>The first order optimality conditions for (19) and (20) are given by (21)</p><disp-formula id="scirp.62642-formula24"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x73.png"  xlink:type="simple"/></disp-formula><p>where e is a vector of ones. The primal-dual central path method can be used to solve the convex QP. Detailed information on this interior point algorithm and other variants can be obtained in Gondzio [<xref ref-type="bibr" rid="scirp.62642-ref8">8</xref>] .</p></sec></sec><sec id="s5"><title>5. BLP and Convex QP Relationship</title><disp-formula id="scirp.62642-formula25"><graphic  xlink:href="http://html.scirp.org/file/1-1040416x74.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62642-formula26"><graphic  xlink:href="http://html.scirp.org/file/1-1040416x75.png"  xlink:type="simple"/></disp-formula><p>From the two versions of the same problem</p><disp-formula id="scirp.62642-formula27"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x76.png"  xlink:type="simple"/></disp-formula></sec><sec id="s6"><title>6. Numerical Illustration</title><p>The following numerical illustration shows how a BLP problem is transformed into convex quadratic programming model and then solved.</p><sec id="s6_1"><title>6.1. Pure Binary Linear Programming</title><disp-formula id="scirp.62642-formula28"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x77.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x78.png" xlink:type="simple"/></inline-formula> are binary variables.</p><p>Transforming into a convex quadratic programming problem becomes (24)</p><disp-formula id="scirp.62642-formula29"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x79.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x80.png" xlink:type="simple"/></inline-formula> are also binary variables.</p><p>The solution to the convex quadratic problem is given in (25).</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x81.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x82.png" xlink:type="simple"/></inline-formula> (25)</p></sec><sec id="s6_2"><title>6.2. Mixed Binary Linear Programming Problem</title><p>In the case of a mixed binary linear programming problem, only the binary integer variables occupy the enforcer. In other words, if only the r binary variables <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x83.png" xlink:type="simple"/></inline-formula> are integer then use</p><disp-formula id="scirp.62642-formula30"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x84.png"  xlink:type="simple"/></disp-formula><p>Suppose in 5.1, the variables <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x85.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x86.png" xlink:type="simple"/></inline-formula> are not restricted to integer but both variables are less than 1.</p><p>The transformation becomes as shown in (27).</p><disp-formula id="scirp.62642-formula31"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x87.png"  xlink:type="simple"/></disp-formula><p>The solution to the convex quadratic problem is:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x88.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x89.png" xlink:type="simple"/></inline-formula> (28)</p></sec></sec><sec id="s7"><title>7. From Mixed Integer Problem to BLP</title><p>The problems that occur in real life do not have binary variables only. These practical problems occur as general mixed integer problem (MIP) where variables assume integer values greater than 1. There are methods that can be used to solve these problems but we are not aware of any method that can solve these mixed integer problems in polynomial time up now. The obvious strategy is to expand the general mixed integer variable into binary ones.</p><sec id="s7_1"><title>7.1. Converting MIP into BLP</title><p>Any MIP variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x90.png" xlink:type="simple"/></inline-formula> can be expanded into binary variables as given in (29).</p><disp-formula id="scirp.62642-formula32"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x91.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x92.png" xlink:type="simple"/></inline-formula> is a binary variable for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x93.png" xlink:type="simple"/></inline-formula> This procedure is explained in Owen and Mehrotra [<xref ref-type="bibr" rid="scirp.62642-ref9">9</xref>] .</p></sec><sec id="s7_2"><title>7.2. Numerical Illustration</title><p>Convert the following MIP into a BLP.</p><disp-formula id="scirp.62642-formula33"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x94.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x95.png" xlink:type="simple"/></inline-formula> are integers.</p><p>The following substitutions change the problem into a BLP.</p><disp-formula id="scirp.62642-formula34"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1040416x96.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x97.png" xlink:type="simple"/></inline-formula> is a binary variable for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x98.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1040416x99.png" xlink:type="simple"/></inline-formula></p></sec></sec><sec id="s8"><title>8. Conclusion</title><p>The general BLP problem has been given so much attention by researchers all over the world for over half a century without a breakthrough. A difficult category of BLP models includes the traveling salesman, generalized assignment, quadratic assignment and set covering problems. The paper presented a technique to solve BLP problems by first transforming them into convex QPs and then applying interior point algorithms to solve them in polynomial time. We also showed that the proposed technique worked for both pure and mixed BLPs and also for the general linear integer model where variables were expanded into BLPs. We hope the proposed approach will give more clues to researchers in the hunt for efficient solutions to the general difficult integer programming problem.</p></sec><sec id="s9"><title>Acknowledgements</title><p>The author is thankful to the referees for their helpful and constructive comments.</p></sec><sec id="s10"><title>Cite this paper</title><p>EliasMunapo, (2016) Solving the Binary Linear Programming Model in Polynomial Time. American Journal of Operations Research,06,1-7. doi: 10.4236/ajor.2016.61001</p></sec></body><back><ref-list><title>References</title><ref id="scirp.62642-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Fortnow, F. (2009) The Status of the P versus NP Problem. Communications of the ACM, 52, 78-86. &lt;/br&gt;http://dx.doi.org/10.1145/1562164.1562186</mixed-citation></ref><ref id="scirp.62642-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Fortnow, F. (2013) The Golden Ticket: P, NP, and the Search for the Impossible. Princeton University Press, Princeton. &lt;/br&gt;http://dx.doi.org/10.1515/9781400846610</mixed-citation></ref><ref id="scirp.62642-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Adams, W.P. and Sherali, H.D. (1990) Linearization Strategies for a Class of Zero-One Mixed Integer Programming Problems. Operations Research, 38, 217-226. &lt;/br&gt;http://dx.doi.org/10.1287/opre.38.2.217</mixed-citation></ref><ref id="scirp.62642-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Freund, R.M. (2002) Solution Methods for Quadratic Optimization: Lecture Notes. Massachusetts Institute of Technology, Cambridge, MA.</mixed-citation></ref><ref id="scirp.62642-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Jensen, P.A and Bard, J.F. (2012) Operations Research Models and Methods. John Wiley &amp;Sons, Inc., Hoboken, NJ.</mixed-citation></ref><ref id="scirp.62642-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Taha, H.A. (2004) Operations Research: An Introduction. 7th Edition, Pearson Educators, New Delhi.</mixed-citation></ref><ref id="scirp.62642-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Winston, W.L. (2004) Operations Research Applications and Algorithms. 4th Edition, Duxbury Press, Ontario.</mixed-citation></ref><ref id="scirp.62642-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Gondzio, J. (2012) Interior Point Methods 25 Years Later. European Journal of Operational Research, 218, 587-601. &lt;/br&gt;http://dx.doi.org/10.1016/j.ejor.2011.09.017</mixed-citation></ref><ref id="scirp.62642-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Owen, J.H. and Mehrotra, S. (2002) On the Value of Binary Expansions for General Mixed-Integer Linear Programs. Operations Research, 50, 810-819. &lt;/br&gt;http://dx.doi.org/10.1287/opre.50.5.810.370</mixed-citation></ref></ref-list></back></article>