See	discussions,	stats,	and	author	profiles	for	this	publication	at:	https://www.researchgate.net/publication/266911500
A	Improved	Vogel's	Approximatio	Method	for
the	Transportation	Problem
Article		in		Mathematical	and	Computational	Applications	·	August	2011
DOI:	10.3390/mca16020370
CITATIONS
12
2	authors:
READS
660
Serdar	Korukoğlu
Ege	University
Serkan	Balli
Mugla	Üniversitesi
49	PUBLICATIONS			383	CITATIONS			
45	PUBLICATIONS			118	CITATIONS			
SEE	PROFILE
SEE	PROFILE
Some	of	the	authors	of	this	publication	are	also	working	on	these	related	projects:
Tabakalı	Kompozit	Plaklardaki	Hasar	Durumunun	Yapay	Zeka	Teknikleri	ile	Analiz	Edilmesi	View
project
All	content	following	this	page	was	uploaded	by	Serkan	Balli	on	30	October	2015.
The	user	has	requested	enhancement	of	the	downloaded	file.
 
 
 
 
 Mathematical and Computational Applications, Vol. 16, No. 2, pp. 370-381, 2011.  
 © Association for Scientific Research 
Aҭ IMPROVED VOGEL’S APPROXIMATIOҭ METHOD  
FOR THE TRAҭSPORTATIOҭ PROBLEM 
1Department of Computer Engineering, Ege University, 35000, Bornova, Đzmir, Turkey 
Serdar Korukoğlu1 and Serkan Ballı2 
 
2Department of Statistics, Muğla University, 48187, Muğla, Turkey 
serdar.korukoglu@ege.edu.tr 
serkan@mu.edu.tr 
 
Abstract- Determining efficient solutions for large scale transportation problems is an 
important  task  in  operations  research.  In  this  study,  Vogel’s  Approximation  Method 
(VAM)  which  is  one  of  well-known  transportation  methods  in  the  literature  was 
investigated to obtain more efficient initial solutions.  A variant of VAM was proposed 
by using total opportunity cost and regarding alternative allocation costs. Computational 
experiments were carried out to evaluate VAM and improved version of VAM (IVAM). 
It was seen that  IVAM  conspicuously obtains more efficient initial solutions for large 
scale transportation problems. Performance of IVAM over VAM was discussed in terms 
of iteration numbers and CPU times required to reach the optimal solutions. 
 
Keywords-  Transportation  Problem,  Integer  Programming,  Vogel’s  Approximation 
Method, Total Opportunity Cost, Simulation Experiments 
 
1. IҭTRODUCTIOҭ 
 
 
The  transportation  problem  is  a  special  kind  of  the  network  optimization 
problems.  It has the special data structure in solution characterized as a  transportation 
graph.  Transportation  models  play  an  important  role  in  logistics  and  supply  chains.   
The  problem  basically  deals  with  the  determination  of  a  cost  plan  for  transporting  a 
single  commodity  from  a  number  of  sources  to  a  number  of  destinations  [16].  The 
purpose is to minimize the cost of shipping goods from one location to another so that 
the  needs  of  each  arrival  area  are  met  and  every  shipping  location  operates  within  its 
capacity [10]. Network model of the transportation problem is shown in Figure 1 [17]. It 
aims to find the best way to fulfill the demand of n demand points using the capacities 
of m supply points. 
 
Cij
Supplies
available
S1
S2
S3
Sm
.
.
.
D1
D2
D3
Dn
.
.
Demand
requirements
 
Figure 1. Network model of the transportation problem 
 
 
S. Korukoğlu and S. Ballı 
371 
ijx is number of 
In Figure 1, S1-Sm are sources and D1-Dn are destinations. 
units shipped from supply point i to demand point j then the general linear programming 
representation of a transportation problem is : 
ijc is cost and 
m
n
    
min
∑ ∑  
c x
ij
ij
1
=
i
1
=
j
n
subject  to    
≤∑
   
S
x
ij
  (i=1,2,...,m)  
Supply constraints                                       (1) 
i
1
=
j
n
                   
≥∑
 
   
x D
ij
j
  (j=1,2,...,n)
   Demand constraints                                     (2) 
1
=
j
 
ijx ≥
0  and integer  for all  i,j 
                         
 
If  total  supply  equals  total  demand  then  the  problem  is  said  to  be  a  balanced 
transportation problem. The reader may refer to Wagner [17] and Taha [16] for detailed 
coverage of transportation problem. 
 
Transportation  problems  can  be  solved  by  using  general  simplex  based  integer 
programming  methods,  however  it  involves  time-consuming  computations.  There  are 
specialized algorithms for transportation problem that are much more efficient than the 
simplex algorithm [18].  The basic steps to solve transportation problem are: 
Step 1. Determination the initial feasible solution, 
Step 2. Determination optimal solution using the initial solution. 
 
In this study, basic idea is to get better initial solutions for the transportation problem. 
Therefore, study focused on Step 1 above. Several heuristic methods are available to get 
an  initial  basic  feasible  solution.  Although  some  heuristics  can  find  an  initial  feasible 
solution  very  quickly,  oftentimes  the  solution  they  find  is  not  very  good  in  terms  of 
minimizing  total  cost.  On  the  other  hand,  some  heuristics  may  not  find  an  initial 
solution as quickly, but the solution they find is often very good in terms of minimizing 
total  cost  [2].  Well-known  heuristics  methods  are  North  West  Corner  [4],  Best  Cell 
Method,  Vogel’s  Approximation  Method  (VAM)  [11],  Shimshak  et.  al.'s  version  of 
VAM [14], Goyal's version of VAM [6], Ramakrishnan's version of VAM [9] etc. Kirca 
and  Satir  [7]  developed  a  heuristic  to  obtain  efficient  initial  basic  feasible  solutions, 
called  Total  Opportunity-cost  Method  (TOM).  Balakrishnan  [3]  proposed  a  modified 
version  of  VAM  for  unbalanced  transportation  problems.  Gass  [5]  reviewed  various 
methods and discussed on solving the transportation problem. Sharma and Sharma [12] 
proposed  a  new  procedure  to  solve  the  dual  of  the  well-known  uncapacitated 
transportation  problem.  Sharma  and  Prasad  [13]  proposed  heuristic  gives  significantly 
better solutions than the well-known VAM. This is a best heuristic method than Vogel’s 
to  get  initial  solution  to  uncapacitated  transportation  problem.  Adlakha  and  Kowalski 
[1]  presented  a  simple  heuristic  algorithm  for  the  solution  of  small  fixed-charge 
transportation problems. Mathirajan and Meenakshi [8] were extended TOM using the 
VAM  procedure.  They  coupled  VAM  with  total  opportunity  cost  and  achieved  very 
efficient initial solutions. 
 
 
An Improved Vogel’s Approximation Method  
372 
 
In  this  paper,  VAM  was  improved  by  using  total  opportunity  cost  and  regarding 
alternative  allocation  costs.  Mathirajan  and  Meenakshi  [8]  applied  VAM  on  the  total 
opportunity cost matrix. In addition to this method, improved VAM (IVAM) considers 
highest  three  penalty  costs  and  calculates  alternative  allocation  costs  in  VAM 
procedure. Then it selects minimum one of them. 
 
Paper  is  organized  as  follows.  VAM  is  summarized  and  illustrated  with  solving  a 
sample transportation problem in the following section. IVAM is explained in the third 
section.  In  the  fourth  section,  simulation  experiments  are  given  and  performance  of 
IVAM over VAM is discussed. Results are clarified in fifth section. 
 
 
 
2. VOGEL’S APPROXIMATIOҭ METHOD (VAM) 
VAM  is  a  heuristic  and  usually  provides  a  better  starting  solution  than  other 
methods.  Application  of  VAM  to  a  given  problem  does  not  guarantee  that  an  optimal 
solution  will  result.  However,  a  very  good  solution  is  invariably  obtained  with 
comparatively little effort [15]. In fact, VAM generally  yields an optimum or close to 
optimum starting solution for small sized transportation problems [16].  
 
VAM is based on the concept of penalty cost or regret. A penalty cost is the difference 
between  the  largest  and  next  largest  cell  cost  in  a  row  or  column.  VAM  allocates  as 
much as possible to the minimum cost cell in the row or column with the largest penalty 
cost. Detailed processes of VAM are given below: 
 
 
Step 1:  Balance  the  given  transportation  problem  if  either  (total  supply>total 
demand) or (total supply
 
 
S. Korukoğlu and S. Ballı 
373 
Table 1. An example of 5x5 transportation problem 
From/To 
D1 
D2 
D3 
D4 
D5 
Supply 
S1 
S2 
S3 
S4 
S5 
46 
12 
35 
61 
85 
 
 
 
 
 
74 
75 
199 
81 
60 
 
 
 
 
 
9 
6 
4 
44 
14 
 
 
 
 
 
28 
36 
5 
88 
25 
 
 
 
 
 
99 
48 
71 
9 
79 
Demand 
278 
60 
461 
116 
1060 
 
 
 
 
 
461 
277 
356 
488 
393 
 
 
 
Initial basic solution for this problem was obtained using VAM and  given in Table 2. 
Using the values in Table 2, initial cost was calculated as 68804. 
 
Table 2. Initial solution tableau of VAM 
From/To 
D1 
D2 
D3 
D4 
D5 
Supply 
S1 
S2 
S3 
S4 
S5 
46 
12 
35 
61 
85 
1 
74 
60 
9 
68 
28 
 
99 
332 
461 
277 
 
 
 
75 
199 
81 
60 
 
 
 
 
6 
4 
44 
14 
 
 
 
393 
36 
5 
88 
25 
 
116 
 
 
48 
71 
9 
79 
 
240 
488 
 
277 
356 
488 
393 
Demand 
278 
60 
461 
116 
1060 
 
 
 
Optimal  solution  was  achieved  using  transportation  simplex  algorithm  [4]  within  five 
iterations and final cost was found as 59356. Optimal solution tableau is given in Table 
3.  Solution  of  proposed  method  (IVAM)  for  the  same  problem  is  illustrated  in  the 
following section. 
 
Table 3. Optimal solution tableau 
From/To 
D1 
D2 
D3 
D4 
D5 
Supply 
S1 
S2 
S3 
S4 
S5 
46 
12 
35 
61 
85 
 
277 
1 
 
 
74 
75 
199 
81 
60 
 
 
 
 
60 
9 
6 
4 
44 
14 
461 
 
 
 
 
28 
36 
5 
88 
25 
 
 
116 
 
 
99 
48 
71 
9 
79 
 
 
239 
488 
333 
461 
277 
356 
488 
393 
Demand 
278 
60 
461 
116 
1060 
 
 
 
 
 
 
 
An Improved Vogel’s Approximation Method  
374 
 
3. IMPROVED VAM (IVAM) 
 
VAM  was  improved  by  using  total  opportunity  cost  (TOC)  matrix  and  regarding 
alternative allocation costs. The TOC matrix is obtained by adding the "row opportunity 
cost matrix" (row opportunity cost matrix: for each row, the smallest cost of that row is 
subtracted  from  each  element  of  the  same  row)  and  the  "column  opportunity  cost 
matrix" (column opportunity cost matrix: for each column of the original transportation 
cost matrix the smallest cost of that column is subtracted from each element of the same 
column) [8]. Proposed algorithm is applied on the TOC matrix and it considers highest 
three penalty costs and calculates alternative allocation costs in VAM procedure. Then 
it selects minimum one of them. Detailed processes are given below: 
 
Step 1:  Balance the given transportation problem if either (total supply>total demand) or (total 
supply
 
 
S. Korukoğlu and S. Ballı 
375 
4. SIMULATIOҭ EXPERIMEҭTS 
 
For  evaluating  the  performance  of  the  VAM  and  its  variant  IVAM,  simulation 
experiments were carried out on a 2.13 GHz Intel Core 2 Duo machine with 4096 MB 
RAM. The main goal of the experiment was to evaluate the effectiveness of the initial 
solutions  obtained  by  VAM  and  IVAM  by  comparing  them  with  optimal  solutions. 
Effectiveness indicates closeness degree which is the lowest iteration number between 
initial solution and the optimal solution. Measures of effectiveness are explained below. 
 
4.1 Measure of effectiveness 
The performances of VAM and IVAM are compared using the following measures: 
Average Iteration (AI): Mean of iteration numbers to obtain optimal solutions using the 
initial solutions of VAM and IVAM over various sized problem instances. 
̱umber of best solutions (̱BS): A frequency which indicates the number of instances 
VAM and IVAM yielded optimal solution with lower iteration over the total of problem 
instances. NBS does not contain case of equal iteration between VAM and IVAM. 
Computation Time:  The CPU time is represented by three variables: T1, T2 and T3. T1 is 
the  time  to  reach  initial  solution.  T2  is  the  time  to  reach  optimal  solution  from  initial 
solution and T3 is the total time from the beginning that is sum of T1 and T2. 
 
4.2 Experimental design 
The transportation problems were randomly generated with twelve different sizes (row 
x  column):  5x5,  10x10,  10x20,  10x30,  10x40,  20x20,  10x60,  30x30,  10x100,  40x40, 
50x50  and  100x100,  respectively.  The  performance  of  the  VAM  and  IVAM  were 
compared over 1000 problem instances for each different sized problem. Total problem 
instances  were  12000.    Costs  were  generated  as  uniformly  discrete  in  the  range  of 
(0,1000)  and  supplies  and  demands  were  generated  in  the  range  of  (0,100)  uniformly 
discrete. All the 12000 problem instances were balanced. The experimental design was 
implemented using ANSI C. 
 
4.3 Comparison of VAM and IVAM 
The experiments and the analysis of the experimental data are presented in this section. 
For each problem instance, a linear programming model was implemented and solved. 
In  order  to  get  a  linear  programming  model  for  each  problem  instance,  a  matrix 
generator procedure and VAM and IVAM were implemented using ANSI C. For each 
problem  instance,  the  heuristic  solutions  were  obtained  using  VAM  and  IVAM.  The 
performance  of  the  VAM  and  IVAM  in  comparison  with  the  optimal  solution  is 
presented below. 
 
Performance  measure  -  AI:  AI  and  other  statistics  for  the  iteration  numbers  of  VAM 
and  IVAM over various sized problems are given in Table 5. For each different sized 
problem, statistics are calculated using 1000 problem instances. From Table 5, it is seen 
that  AI  of  VAM  is  better  than  IVAM  for  small  sized  transportation  problems  but  it 
deteriorates  when  problem  size  increases.  IVAM  yields  more  efficient  results  than 
VAM for large sized problems. 
 
 
 
An Improved Vogel’s Approximation Method  
376 
 
Table 5. Statistical indicators for the iteration numbers of VAM and IVAM 
Problem 
AI 
Standard Error  Median 
Range 
size  VAM  IVAM  VAM 
IVAM  VAM  IVAM VAM IVAM 
2.198  2.676  0.034 
5x5 
10x10  6.155  6.359  0.069 
10x20  10.063  9.913  0.094 
10x30  13.264  12.951  0.120 
10x40  17.761  17.495  0.141 
20x20  17.007  16.337  0.146 
10x60  22.869  22.118  0.173 
30x30  29.912  28.363  0.210 
10x100  32.561  31.139  0.221 
40x40  44.801  42.603  0.291 
50x50  60.505  57.651  0.346 
100x100 158.890 149.53  0.657 
0.035 
0.066 
0.092 
0.116 
0.144 
0.135 
0.166 
0.200 
0.220 
0.269 
0.328 
0.610 
2 
6 
10 
13 
17 
17 
23 
30 
32 
44 
60 
158 
5 
3 
14 
6 
18 
10 
24 
13 
28 
17 
31 
16 
36 
22 
56 
28 
40 
31 
55 
42 
74 
57 
149  137 
5 
14 
18 
27 
30 
31 
32 
39 
43 
59 
60 
112 
 
Both parametric and nonparametric statistical tests were performed using MINITAB-15 
Statistical  Package  for  comparing  iteration  numbers  of  VAM  and  IVAM.  Firstly, 
Student's  t-test  was  used  for  testing  the  mean  of  differences  of  the  iteration  numbers 
between VAM and IVAM based on 1000 samples. Secondly Wilcoxon test was used for 
testing the median of differences of the iteration numbers between VAM and IVAM on 
the  same  samples.    The  formal  test  is  Ho:  No  statistically  significant  difference  in  the 
amount  of  mean  (median)  to  complete  the  transportation  optimization  between  the 
VAM and IVAM initial methods. Table 6 gives a summary of the results of two-sided 
statistical tests.  
Table 6. Statistical tests for difference of iteration numbers between VAM and IVAM 
 
Problem 
size 
5x5 
10x10 
10x20 
10x30 
10x40 
20x20 
10x60 
30x30 
10x100 
40x40 
50x50 
100x100 
 Student's t-test 
Mean ± 
Standard Error 
Confidence 
Interval %95 
t 
P 
 Wilcoxon Test 
Wilcoxon 
Statistic 
Estimated 
Median 
P 
-0.478±0.039 
-0.204±0.074 
0.150±0.102 
0.313±0.125 
0.266±0.159 
0.670±0.157 
0.751±0.191 
1.549±0.229 
1.422±0.249 
2.198±0.313 
2.854±0.375 
9.360±0.725 
-0.555 ; -0.401  -12.190  0.000 
-0.348 ; -0.060 
-2.770  0.006 
1.470  0.143 
-0.051 ;  0.351 
2.510  0.012 
0.068  ;  0.558 
-0.047 ; 0.579 
1.670  0.096 
4.270  0.000 
0.362  ;  0.978 
3.930  0.000 
0.376  ;  1.126 
1.099  ;  1.999 
6.760  0.000 
5.720  0.000 
0.934  ;  1.910 
7.030  0.000 
1.584  ;  2.812 
2.118  ;  3.590 
7.600  0.000 
7.937  ;  10.783  12.900  0.000 
-0.500  52276.000  0.000 
0.000  148409.000 0.008 
0.000  203812.000 0.103 
0.500  218149.500 0.022 
0.500  229297.500 0.104 
0.500  245585.500 0.000 
1.000  248683.500 0.000 
1.500  281480.500 0.000 
1.500  276839.500 0.000 
2.500  290830.500 0.000 
3.000  293811.500 0.000 
9.000  349416.000 0.000