Startseite A polynomial algorithm for packing unit squares in a hypograph of a piecewise linear function
Artikel Open Access

A polynomial algorithm for packing unit squares in a hypograph of a piecewise linear function

  • Marat Arslanov EMAIL logo , Zhazira Amirgalieva und Chingiz Kenshimov
Veröffentlicht/Copyright: 29. Dezember 2017
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

We consider the problem of packing the maximal number of unit squares in a hypograph of a function. A polynomial time algorithm is described to solve this problem for a piecewise linear function.

1 Introduction

The problems of rectangular packing and cutting are difficult problems of discrete optimization (for references see [1]). For example, the computational complexity of the problem of packing in a rectangular sheet (A, B) the maximum number of equal small rectangles (c, d) (pallet loading problem (PLP)) is unknown. This problem has number 55 in the list of open problems “The Open Problems Project”. It is closely connected with another known problem, 56 in the list of open problems, about the complexity of deciding whether a given number of axis-parallel unit squares can be packed into a simple polygon (without holes). Appearance [2] conjecture the problem to be polynomially solvable. The problem is known to be NP-hard for polygons with holes [3], even if the polygon is an orthogonal polygon with all coordinates being multiples of 1/2. Recently, this version of the problem was shown to be in NP [6], making it NP-complete. This problem arises in wood industry (cutting a log into quadratic plates), in glass industry (cutting defecting pane glass into small equal squares) etc. Applications of analogous problems can be seen in [4, 5]. In this article we consider the problem of packing of unit squares in a hypograph of a function. The article is organized as follows. The second section contains the proof of a theorem about the structure of the optimal packing unit squares in a hypograph of a function. In the third section we propose a polynomial algorithm for the problem with piecewise linear function.

2 A theorem about optimal packing

It appears, that a two-dimensional problem about optimum packing unit squares in a hypograph of a function can be reduced to a one-dimensional problem of packing unit intervals in a one-dimensional set which have a trivial solution. If a compact set is represented as a separable union of intervals

M=iIAi,

then the maximum number of single intervals which are possible to pack in this set is given by the obvious formula

N=iI|Ai|.

Definition 1

A hypograph of a function f on the interval [a, b] is the set

H(f,a,b)={(x,y)|0yf(x),axb}.

Note that this definition is different from the classical definition of hypograph in [7]. Cleave this set by lines parallel to the x-axis with integer y-coordinates. We designate these lines by Li:

Li={(x,y)|y=iN={1,2,}}.

The intersection of the set Li and the hypograph H(f, a, b) is a compact set, which we designate by Hi(f, a, b). Let us suppose that this set consists of finite number of intervals. Let us suppose also without loss of generality that the function f is a bounded continuous function on the interval [a, b]. Then, if Hi(f, a, b) = ⋃jJAi,j, where Ai,j is an interval, the maximum number Ni of unit intervals in Hi(f, a, b) is given by

Ni=jJ|Ai,j|.

The following theorem reduces the packing of unit squares to the packing of unit intervals.

Theorem 1

The maximum number N of unit squares in a hypograph H(f, a, b) of a function f on the interval [a, b] is given by

N=iNNi.(1)

Proof

This theorem is almost obvious. We prove this claim by induction on the area of a hypograph. If the area of a hypograph is less than 1, then any unit square cannot be packed in such a hypograph. In the same way, any interval Ai,j = [ai,j, bi,j] has length less than 1, otherwise, a rectangle

Pi,j={(x,y)|ai,jxbi,j,i1yi},

lying inside a hypograph H(f, a, b), contains unit square. Let for all hypographs of area less then S, theorem is valid. Let us show its justice for hypographs of area greater then or equal to S. Let us consider a set H1(f, a, b) = ⋃jJA1,j. It can cross no more than

N1=jJ|A1,j|

unit squares. All remaining unit squares are above the line y = 1. Then these squares with decreased by 1 y-coordinate are in the hypograph of a function, f (x) – 1. By the inductive hypothesis, in the hypograph of a function f (x) – 1, it can be packed ∑iN\{1}Ni unit squares. Together with squares intersecting a line y = 1 it can be proved easily that it is impossible to locate in the hypograph H(f, a, b) more than N = ∑iNNi unit squares. But such a number of unit squares can be locate, if the first N1 squares we locate in a hypograph between lines y = 0 and y =1. The theorem is proved.

Let’s illustrate this theorem by a problem of optimal packing of unit squares in a semicircle of radius 5 (hypograph of a function 25x2). We have N1=96,N2=84,N3 = ⌊8⌋, N4 = ⌊6⌋. According to theorem, no more than 32=9+9+8+6 unit squares can be packed in a semicircle of radius 5 and optimal packing is presented in Figure 1.

Figure 1 Packing unit squares in a semicircle
Figure 1

Packing unit squares in a semicircle

As another application of this theorem, we consider a problem of packing unit squares in a hypograph H(f, a, b) of a function F(x)=nx.

If this problem can be solved in a polynomial time, then a known unsolved problem of factoring integers in a polynomial time has a positive answer. Let’s designate N(f, a, b) as the maximum number of unit squares, which can be packed into a hypograph of a function H(f, a, b), then

N(nx,a,b)N(n1x,a,b) is equal to the number of integral points on the graph of the function f(x)=nx,axb. As the points (1, n), (n, 1) lie on the graph of this function, an inequality N(nx,1,n)N(n1x,1,n)3 tests the primality of n. Then, to find prime factors of n, one can apply the method of half division of the interval [1, n]. Thus, the complexity of a problem of optimum packing of unit squares in a hypograph of a function H(nx,a,b) is limited from below by the complexity of factoring integers for which the best algorithms to date do not run in a polynomial time.

3 Polynomial algorithm for a piecewise-linear function

In this section we consider a problem of optimum packing of unit squares in a hypograph of a piecewise-linear function and show the existence of a polynomial algorithm for its solution.

Let us consider in the beginning, the case of trapezioidal function:

f(x)=c1+d1x,axa1c,a1xa2c2+d2x,a2xb.

Coefficients of this function satisfy the equalities

c1+d1a1=c=c2+d2a2,c1+d1a=c2+d2b=0,d1>0,d2<0.

The graph of this function is presented in Figure 2.

Figure 2 Trapezoidal function
Figure 2

Trapezoidal function

Then, by Theorem 1, a solution to the problem of packing of unit squares in a hypograph of a function H(f, a, b) is given by the formula

N(f,a,b)=i=1cNi,

where

Ni=bai(baa2+a1)c.

This expression is equal to the number of points of a twodimensional integral lattice without the points on axes of the abscissas in a trapezoid with vertices in points (0, 0), (0, c), (a2a1, c), (ba, 0). The formula for the number of integral points in a trapezoid in the vertices (0, 0), (0, a), (b, a), (c, 0) is given by

N=i=1aci(cb)a

Thus, the problem of optimum packing of unit squares in a hypograph of a trapezoidal function is reduced to the problem of counting the number of integral points in a trapezoid. This last problem in the general case of polygons is widely investigated and the existence of polynomial time algorithms is well known (see e.g. [8, 9]). So, for a solution of this reduced problem there exist effective polynomial algorithms. Thus, the problem of optimum packing of unit squares in the hypograph of trapezoidal function has effective, polynomial algorithms.

Let’s consider the case of a piecewise linear function

F(x)=c1+d1x,ifaxa1c2+d2x,ifa1xa2ifcn+dnx,ifan1xan=b.

Such a function can be represented as an ordered sequence of vertices of a hypograph considered as a polygon with vertices

(a1,b1),(a2,b2),(an,bn),a1<a2<<an.

It is easy to see that O(n) arithmetical operations are sufficient for transforming from one form of presentation to another. Moreover, a problem of optimum packing unit squares in a hypograph of a piecewise linear function can be considered by virtue of the theorem 1, as a problem of counting their numbers with the help of formulas (1), i.e. as counting the sum of integer part of lengths of intervals cutting from a hypograph of a piecewise linear function by the set of lines y = i, iN. Then a problem of counting the sum of integer part of lengths of intervals cutting from a hypograph of a piecewise linear function, considered as a sequence of vertices of a polygon

(a1,b1),(a2,b2),(an,bn),a1<a2<<an.

by the set of lines y = i, iN is easily decomposed into problems with smaller numbers of vertices, and the sum of such numbers of vertices will be no more than n. The main idea consists in cutting a trapezoidal part of a hypograph of a non-unimodal piecewise linear function. Before formulating an algorithm for general case, we shall consider the case of a unimodal piecewise linear function, for which a sequence bi is of the form:

b1b2bkbk+1bn1bn.

Such a polygon can be easily presented in a way in which

b1=bn1,b2=bn2,bk=bnk

Then, counting the sum of the integer part of lengths of intervals cutting from a hypograph of such a function can be reduced to the formula

N=i=1n/2Ni,

where Ni is equal to the sum of the integer part of lengths of intervals cutting from a hypograph of trapezoidal function as an appropriate trapezoid with vertices

(ai,bi),(ai+1,bi+1),(ani1,bni1),(ani,bni).

As we already have stated, for calculating Ni there is a polynomial algorithm based on counting the points of integral lattice in an appropriate trapezoid. Therefore, there is a polynomial algorithm for counting the sum of the integer part of lengths of intervals cutting from a hypograph of a unimodal piecewise linear function by lines y = i, iN.

Let’s consider the general case of a piecewise linear function, presented by the sorted list of vertices:

(a1,b1),(a2,b2),(an,bn),a1<a2<<an.

If this function is not unimodal on the interval [a1, an], there is a local minimum of this function bk:

bk1>bk<bk+1.
Figure 3 Piecewise linear function
Figure 3

Piecewise linear function

Select the least among such local minima. Let this minimum has an index k. Then by a line y = bk, we cut a hypograph into two parts. The bottom part corresponds to a unimodal function. The upper part corresponds to a function on a smaller interval, and this smaller interval can be divided into two intervals, on each of which it is possible to solve a problem of counting the sum of integer part of lengths of intervals cutting by lines y = i. Thus, the following decomposition has been derived (see figure 2):

N=N1+N2+N3,

where N1 corresponds to a solution of a problem for unimodal function, N2 is a solution of a problem for function on an interval [a′, ak], N3 is a solution of a problem for function on an interval [ak, b′], a′, b′ are the left and right intersection of a line Y = bk with the graph of considered function.

Then the time of solution T(n) of considered problem satisfies the recurrent ratio

T(n)=P(n)+T(k)+T(nk),

where P(n) is the time of the solution of a problem for trapezoidal component of a function plus elementary operations on calculation of k, a′, b′. It is easy to show that a solution of the last recurrent ratio requires an easily proved polynomial evaluation T(n). Thus, taking into account that a problem of calculating the sum of integer parts of intervals cutting by lines y = i in a hypograph of function is equivalent to a problem of optimum packing of unit squares in a hypograph of this function, the following theorem is proved.

Theorem 2

For a problem of optimum packing of unit squares in a hypograph H(f, a, b) of a piecewise linear function f on the interval [a, b] there exists a polynomial algorithm.

Let’s remark, that the crude estimates allow to evaluate a computing time of the given algorithm as quadratic number of arithmetical operations. That is for a hypograph of function, consisting n vertices, the number of arithmetical operations required for calculating the maximum number of unit squares, which can be packed into the hypograph, is evaluated by O(n2).

4 Conclusions

The proved theorems make probable a hypothesis about an existence of polynomial algorithm for a general problem of optimum packing the unit squares in a polygon without holes. It is natural to consider at first the following weaker problems.

Problem 1

To prove the existence of a polynomial algorithm for problems of packing of unit squares in a convex polygon.

Problem 2

To prove the existence of a polynomial algorithm for the problem about packing of squares 2 × 2 in a rectilinear polygon without holes, whose all vertices have integer coordinates.

References

[1] H. Dyckhoff, A typology of cutting and packing problems, European Journal of Operational Research 44(1990) 145–159.10.1016/0377-2217(90)90350-KSuche in Google Scholar

[2] C. Baur and S. P. Fekete, Approximation of geometric dispersion problems, Algorithmica, 30:450-470, 2001.10.1007/BFb0053964Suche in Google Scholar

[3] R. J. Fowler, M. S. Paterson, and S. L. Tanimoto. Optimal packing and covering in the plane are NP-complete. Inform. Process. Lett., 12(3):133-137, 1981.10.1016/0020-0190(81)90111-3Suche in Google Scholar

[4] C. Turan and E.N. Amirgaliev. A robust leaky-LMS algorithm for sparse system identification, Proceeding of the 9th international Conference on Discrete Optimization and Operations Research, Vladivostok, Russia, 2016. PP 538-546.10.1007/978-3-319-44914-2_42Suche in Google Scholar

[5] Rysbaiyly B., Yunicheva N. An iterative method to calculate the thermal characteristics of the rock mass with inaccurate initial data. Open Engineering. 2016. PP 627-636.10.1515/eng-2016-0089Suche in Google Scholar

[6] Muriel Dulieu, Dania El-Khechen, John Iacono, and Nikolaj van Omme. Packing 2 × 2 unit squares into grid polygons is NP-complete. In Proceedings of the 21st Canadian Conference on Computational Geometry, Vancouver, Canada, August 2009.Suche in Google Scholar

[7] Rockafellar, Ralph Tyrell (1996), Convex Analysis, Princeton University Press, Princeton, NJ.Suche in Google Scholar

[8] Alexander I. Barvinok. A polynomial time algorithm for counting integral points in polyhedra when the dimension is fixed. Mathematics of Operations Research, 19(4):769-779, 1994.10.1287/moor.19.4.769Suche in Google Scholar

[9] M. Beck and S. Robins. Explicit and efficient formulas for the lattice point count in rational polygons using Dedekind-Rademacher sums. Discrete and Computational Geometry, 27(4):443-459, 2002.10.1007/s00454-001-0082-3Suche in Google Scholar

Received: 2017-6-13
Accepted: 2017-11-19
Published Online: 2017-12-29

© 2017 M. Z. Arslanov et al.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Artikel in diesem Heft

  1. Regular Articles
  2. The Differential Pressure Signal De-noised by Domain Transform Combined with Wavelet Threshold
  3. Regular Articles
  4. Robot-operated quality control station based on the UTT method
  5. Regular Articles
  6. Regression Models and Fuzzy Logic Prediction of TBM Penetration Rate
  7. Regular Articles
  8. Numerical study of chemically reacting unsteady Casson fluid flow past a stretching surface with cross diffusion and thermal radiation
  9. Regular Articles
  10. Experimental comparison between R409A and R437A performance in a heat pump unit
  11. Regular Articles
  12. Rapid prediction of damage on a struck ship accounting for side impact scenario models
  13. Regular Articles
  14. Implementation of Non-Destructive Evaluation and Process Monitoring in DLP-based Additive Manufacturing
  15. Regular Articles
  16. Air purification in industrial plants producing automotive rubber components in terms of energy efficiency
  17. Regular Articles
  18. On cyclic yield strength in definition of limits for characterisation of fatigue and creep behaviour
  19. Regular Articles
  20. Development of an operation strategy for hydrogen production using solar PV energy based on fluid dynamic aspects
  21. Regular Articles
  22. An exponential-related function for decision-making in engineering and management
  23. Regular Articles
  24. Usability Prediction & Ranking of SDLC Models Using Fuzzy Hierarchical Usability Model
  25. Regular Articles
  26. Exact Soliton and Kink Solutions for New (3+1)-Dimensional Nonlinear Modified Equations of Wave Propagation
  27. Regular Articles
  28. Entropy generation analysis and effects of slip conditions on micropolar fluid flow due to a rotating disk
  29. Regular Articles
  30. Application of the mode-shape expansion based on model order reduction methods to a composite structure
  31. Regular Articles
  32. A Combinatory Index based Optimal Reallocation of Generators in the presence of SVC using Krill Herd Algorithm
  33. Regular Articles
  34. Quality assessment of compost prepared with municipal solid waste
  35. Regular Articles
  36. Influence of polymer fibers on rheological properties of cement mortars
  37. Regular Articles
  38. Degradation of flood embankments – Results of observation of the destruction mechanism and comparison with a numerical model
  39. Regular Articles
  40. Mechanical Design of Innovative Electromagnetic Linear Actuators for Marine Applications
  41. Regular Articles
  42. Influence of addition of calcium sulfate dihydrate on drying of autoclaved aerated concrete
  43. Regular Articles
  44. Analysis of Microstrip Line Fed Patch Antenna for Wireless Communications
  45. Regular Articles
  46. PEMFC for aeronautic applications: A review on the durability aspects
  47. Regular Articles
  48. Laser marking as environment technology
  49. Regular Articles
  50. Influence of grain size distribution on dynamic shear modulus of sands
  51. Regular Articles
  52. Field evaluation of reflective insulation in south east Asia
  53. Regular Articles
  54. Effects of different production technologies on mechanical and metallurgical properties of precious metal denture alloys
  55. Regular Articles
  56. Mathematical description of tooth flank surface of globoidal worm gear with straight axial tooth profile
  57. Regular Articles
  58. Earth-based construction material field tests characterization in the Alto Douro Wine Region
  59. Regular Articles
  60. Experimental and Mathematical Modeling for Prediction of Tool Wear on the Machining of Aluminium 6061 Alloy by High Speed Steel Tools
  61. Special Issue on Current Topics, Trends and Applications in Logistics
  62. 10.1515/eng-2017-0001
  63. Special Issue on Current Topics, Trends and Applications in Logistics
  64. The Methodology of Selecting the Transport Mode for Companies on the Slovak Transport Market
  65. Special Issue on Current Topics, Trends and Applications in Logistics
  66. Determinants of Distribution Logistics in the Construction Industry
  67. Special Issue on Current Topics, Trends and Applications in Logistics
  68. Management of Customer Service in Terms of Logistics Information Systems
  69. Special Issue on Current Topics, Trends and Applications in Logistics
  70. The Use of Simulation Models in Solving the Problems of Merging two Plants of the Company
  71. Special Issue on Current Topics, Trends and Applications in Logistics
  72. Applying the Heuristic to the Risk Assessment within the Automotive Industry Supply Chain
  73. Special Issue on Current Topics, Trends and Applications in Logistics
  74. Modeling the Supply Process Using the Application of Selected Methods of Operational Analysis
  75. Special Issue on Current Topics, Trends and Applications in Logistics
  76. Possibilities of Using Transport Terminals in South Bohemian Region
  77. Special Issue on Current Topics, Trends and Applications in Logistics
  78. Comparison of the Temperature Conditions in the Transport of Perishable Foodstuff
  79. Special Issue on Current Topics, Trends and Applications in Logistics
  80. E-commerce and its Impact on Logistics Requirements
  81. Topical Issue Modern Manufacturing Technologies
  82. Wear-dependent specific coefficients in a mechanistic model for turning of nickel-based superalloy with ceramic tools
  83. Topical Issue Modern Manufacturing Technologies
  84. Effects of cutting parameters on machinability characteristics of Ni-based superalloys: a review
  85. Topical Issue Desktop Grids for High Performance Computing
  86. Task Scheduling in Desktop Grids: Open Problems
  87. Topical Issue Desktop Grids for High Performance Computing
  88. A Volunteer Computing Project for Solving Geoacoustic Inversion Problems
  89. Topical Issue Desktop Grids for High Performance Computing
  90. Improving “tail” computations in a BOINC-based Desktop Grid
  91. Topical Issue Desktop Grids for High Performance Computing
  92. LHC@Home: a BOINC-based volunteer computing infrastructure for physics studies at CERN
  93. Topical Issue Desktop Grids for High Performance Computing
  94. Comparison of Decisions Quality of Heuristic Methods with Limited Depth-First Search Techniques in the Graph Shortest Path Problem
  95. Topical Issue Desktop Grids for High Performance Computing
  96. Using Volunteer Computing to Study Some Features of Diagonal Latin Squares
  97. Topical Issue on Mathematical Modelling in Applied Sciences, II
  98. A polynomial algorithm for packing unit squares in a hypograph of a piecewise linear function
  99. Topical Issue on Mathematical Modelling in Applied Sciences, II
  100. Numerical Validation of Chemical Compositional Model for Wettability Alteration Processes
  101. Topical Issue on Mathematical Modelling in Applied Sciences, II
  102. Innovative intelligent technology of distance learning for visually impaired people
  103. Topical Issue on Mathematical Modelling in Applied Sciences, II
  104. Implementation and verification of global optimization benchmark problems
  105. Topical Issue on Mathematical Modelling in Applied Sciences, II
  106. On a program manifold’s stability of one contour automatic control systems
  107. Topical Issue on Mathematical Modelling in Applied Sciences, II
  108. Multi-agent grid system Agent-GRID with dynamic load balancing of cluster nodes
Heruntergeladen am 14.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/eng-2017-0045/html
Button zum nach oben scrollen