This paper addresses the single-item capacitated lot sizing problem with a 1-breakpoint all-units quantity discount in a monotonic setting where the purchase prices are non-increasing over the planning horizon. For this case, we establish several novel properties of the optimal solution and develop a hybrid dynamic programming approach that maintains a compact representation of the solution space by storing only essential information about the states and using linear equations for intermediate values. Our algorithm runs in \(O(n\log n)\) time, where \(n\) denotes the number of periods. Our result is an improvement over the previous state-of-the-art algorithm, which has an \(O(n^2)\) time complexity.
An O(nlogn) Algorithm for Single-Item Capacitated Lot Sizing with a One-Breakpoint All-Units Discount and Non-Increasing Prices
- Paper ID: 2510.11368
- Title: An O(nlogn) Algorithm for Single-Item Capacitated Lot Sizing with a One-Breakpoint All-Units Discount and Non-Increasing Prices
- Author: Kleitos Papadopoulos
- Classification: cs.DS (Data Structures and Algorithms)
- Publication Date: October 2025 (arXiv preprint)
- Paper Link: https://arxiv.org/abs/2510.11368
This paper investigates the single-item capacitated lot sizing problem with a one-breakpoint all-units discount under monotonically decreasing prices. The author establishes several novel structural properties of optimal solutions and develops a hybrid dynamic programming approach that maintains a compact representation of the solution space by storing only critical state information and computing intermediate values using linear equations. The algorithm achieves a time complexity of O(nlogn), where n denotes the number of time periods, representing a significant improvement over the previous best O(n2) algorithm.
Lot sizing problems occupy a central position in manufacturing and supply chain management, where enterprises must minimize total costs—including procurement and storage costs—while satisfying demand. Variants of this problem are ubiquitous in practical business applications, particularly when quantity discounts are present.
The related work comparison table provided in the paper reveals that existing algorithms have high time complexity:
- Federgruen and Lee (1990): O(n3)
- Atamtürk and Hochbaum (2001): O(n3) to O(n5)
- Malekian et al. (2021): O(n4)
- Down et al. (2021): O(n2) (current best)
The author employs a "gas station analogy" to intuitively explain the problem: time periods correspond to gas stations, inventory corresponds to fuel in the tank, demand corresponds to distances between stations, and item costs correspond to fuel prices. This analogy facilitates understanding the intuition behind algorithm design.
- Established structural properties of optimal solutions: Proved that optimal solutions possess specific mathematical properties under non-decreasing prices and single-breakpoint discounts
- Proposed a hybrid dynamic programming algorithm: Developed a compact solution space representation method based on segment representation
- Achieved O(nlogn) time complexity: A significant improvement over the existing best algorithm of O(n2)
- Designed efficient data structures: Utilized augmented balanced binary search trees to maintain segment information and MV thresholds
Input:
- n time periods, each period t with demand dt
- Storage capacity B(t) and unit holding cost ht
- Hierarchical pricing function: pt(x)={p1,txp2,txif x<Qif x≥Q
- Non-increasing prices: p1,t≥p1,t+1, p2,t≥p2,t+1
Output: Procurement and storage strategy that minimizes total cost
Objective Function:
minx,I∑t=1n(pt(xt)+htIt)
Constraints:
- It=It−1+xt−dt
- 0≤It≤B(t)
- I0=0
The author introduces the concept of "state segments," where each segment contains:
- Explicit states: (f,dp(i,f)), where f is the inventory level and dp(i,f) is the minimum cost to reach that state
- Linear equations: Used to generate implicit states within the segment's coverage range
- Auxiliary information: p(i,f) (price at which p2 fuel was last obtained), r(i,f) (available additional fuel quantity)
Lemma 1 (2Q Bound): At any station i, the optimal solution set Sb(i) containing the first 2Q states includes all states necessary to construct the optimal solution set Sb(i+1).
Lemma 2 (Price Monotonicity): For any two states dp(i,f) and dp(i,w) in Sb(i) with f<w, we have p(i,f)≥p(i,w).
Lemma 3 (Continuity): If state dp(i,f) is used to generate an optimal state in S(i), then the generated states form a continuous interval.
To efficiently handle dominance relationships between segments, the author designed the MV (Maximum Value) threshold:
Regular MV (boundary checkpoint):
MV(S1)=g−fdp(i,g)−dp(i,f)
Terminal MV (when S2 uses p2,i on the right):
MVterm(S1)=L−adp(i,g)+Tp2,i−dp(i,f)−ap(i,f)
Processing at each station includes:
- Pruning with p1,i: Remove dominated adjacent segments
- Forming dp(i+1,0): Compare costs of direct arrival and refueling arrival
- Splitting at d(i,i+1): Partition segments into reachable and unreachable categories
- Batch Updates: Add Q units of p2,i fuel to unreachable segments
- Linear Merging: Merge batch-updated segments with existing segments in the [Q,2Q] interval
- Final Pruning: Perform final dominance checks with p1,i
Traditional dynamic programming requires explicit storage of every possible state, whereas the segment representation in this paper only needs to store critical state points and linear equations, significantly reducing space complexity.
By precomputing MV thresholds and maintaining them in a balanced binary search tree, the algorithm can determine dominance relationships between segments in O(logn) time.
Batch updates and holding cost updates employ lazy labels, avoiding unnecessary computations while maintaining algorithm efficiency.
The paper primarily provides theoretical analysis rather than experimental validation, focusing on proving:
- Correctness: Through induction, the algorithm produces correct 2Q-approximate solution sets at each station
- Time Complexity:
- At most 2 new segments created per station
- Total segment count mi≤2i
- Each operation has time complexity O(logn)
- Total time complexity is O(nlogn)
Segment Count Bound (Lemma 7): The optimal Q-approximate solution set Sb(i) contains at most 2i segments.
Operation Complexity:
- Individual updates: Each removal of dominated segments requires O(logn)
- Batch updates: Lazy label application is O(1)
- Splitting/merging: Standard BST operations are O(logn)
The paper provides rigorous theoretical guarantees:
Theorem 1 (Correctness): Under non-increasing unit prices, a single all-units discount threshold, and linear holding costs, Algorithm 1 correctly computes the 2Q-approximate solution set S(i) and optimal boundary state dp(i+1,0) for each station i.
Theorem 2 (Time Complexity): Under segment bounds mi≤2i and augmented balanced tree primitives, the total time complexity for constructing S(i) from Sb(i) is O(nlogn), with space complexity O(n).
The paper also provides two practical extensions:
- Handling B(i)>2Q cases: Addressed through transient in-place expansion for large capacity queries
- Decision tracking: Mechanism supporting recovery of optimal procurement plans
Research on lot sizing problems dates back decades, with major development directions including:
- Cost function extensions: From linear to piecewise linear and concave functions
- Constraint enrichment: Capacity limits, buyback options, subcontracting, etc.
- Algorithm improvements: Progressive refinement from O(n5) to O(n2)
Building upon Down et al. (2021)'s O(n2) algorithm, this paper achieves a breakthrough improvement to O(nlogn) through the introduction of segment representation and MV threshold mechanisms.
- Algorithm efficiency: Achieved time complexity improvement from O(n2) to O(nlogn)
- Theoretical contribution: Established new structural properties of lot sizing problems under monotonic price environments
- Practical value: The algorithm applies equally to integer and non-integer quantities, with broad applicability
- Price assumptions: Requires non-increasing prices, limiting applicability scope
- Single-breakpoint restriction: Only handles cases with a single discount threshold
- Missing experimental validation: The paper primarily provides theoretical analysis without empirical verification on real data
The author proposes research directions including:
- Investigating broader classes of cost and holding functions that preserve lemma validity
- Extension to multi-breakpoint discount scenarios
- Handling non-monotonic price environments
- Strong theoretical innovation: Segment representation and MV threshold mechanisms are original contributions
- Mathematical rigor: Provides complete correctness and complexity proofs
- High practical value: Significant time complexity improvement has important practical implications
- Clear exposition: The gas station analogy makes complex problems accessible
- Insufficient experimental validation: Lacks practical performance comparisons with existing algorithms
- Limited applicability scope: Non-increasing price assumptions may not always hold in practice
- Implementation complexity: Augmented BST implementation is complex and may impact practical application
- Academic contribution: Provides new theoretical tools and analytical frameworks for lot sizing problems
- Practical value: O(nlogn) complexity makes the algorithm suitable for large-scale problems
- Reproducibility: Paper provides detailed algorithm descriptions and data structure implementations
The algorithm is particularly suitable for:
- Large-scale production planning: Long-term planning problems with many time periods
- Decreasing price environments: Such as technology products and seasonal goods
- Capacity-constrained inventory management: Supply chain management with limited warehouse capacity
The paper cites important works in the lot sizing domain, including:
- Chung and Lin (1988): Early O(n2) algorithm
- Federgruen and Lee (1990): O(n3) method introducing all-units discounts
- Atamtürk and Hochbaum (2001): Handling concave cost functions
- Down et al. (2021): Current best O(n2) algorithm
Overall Assessment: This is a high-quality theoretical computer science paper that achieves significant breakthroughs in lot sizing problems. While lacking experimental validation, its theoretical contributions and algorithmic innovations possess important academic value and practical significance.