1799. Maximize Score after N Operations - Explanation

Prerequisites

Before attempting this problem, you should be comfortable with:

GCD (Greatest Common Divisor) - Computing GCD using Euclidean algorithm and understanding its properties
Backtracking - Exploring all possible pairings through recursive enumeration with state restoration
Bitmask Dynamic Programming - Using bitmasks to represent subsets of elements and memoizing states efficiently

1. Brute Force (Backtracking)

Intuition

We need to perform exactly n operations on an array of 2n elements, where each operation picks two elements, computes their GCD, and multiplies it by the operation number. Since we want to maximize the total score, we should try all possible ways to pair up elements and track the best result.

The backtracking approach explores every possible pairing. At each step, we pick any two unvisited elements, mark them as used, compute the score for this operation, then recursively solve for the remaining elements. After exploring that path, we backtrack by unmarking those elements and trying different pairs.

Algorithm

Create a visit array to track which elements have been used.
Define a recursive dfs(n) function where n is the current operation number (starting from 1).
Base case: if n > N/2, all operations are done, return 0.
For each pair of unvisited indices (i, j):
- Mark both as visited.
- Compute n * gcd(nums[i], nums[j]) plus the recursive result for the next operation.
- Track the maximum score.
- Backtrack by unmarking both indices.
Return the maximum score found.

class Solution:
    def maxScore(self, nums: List[int]) -> int:
        N = len(nums)
        visit = [False] * N

        def dfs(n):
            if n > (N // 2):
                return 0

            res = 0
            for i in range(N):
                if visit[i]:
                    continue
                visit[i] = True
                for j in range(i + 1, N):
                    if visit[j]:
                        continue
                    visit[j] = True
                    g = gcd(nums[i], nums[j])
                    res = max(res, n * g + dfs(n + 1))
                    visit[j] = False
                visit[i] = False

            return res

        return dfs(1)

public class Solution {
    public int maxScore(int[] nums) {
        int N = nums.length;
        boolean[] visit = new boolean[N];
        return dfs(nums, visit, 1, N);
    }

    private int dfs(int[] nums, boolean[] visit, int n, int N) {
        if (n > N / 2) {
            return 0;
        }

        int res = 0;
        for (int i = 0; i < N; i++) {
            if (visit[i]) continue;
            visit[i] = true;
            for (int j = i + 1; j < N; j++) {
                if (visit[j]) continue;
                visit[j] = true;
                int g = gcd(nums[i], nums[j]);
                res = Math.max(res, n * g + dfs(nums, visit, n + 1, N));
                visit[j] = false;
            }
            visit[i] = false;
        }

        return res;
    }

    private int gcd(int a, int b) {
        return b == 0 ? a : gcd(b, a % b);
    }
}

class Solution {
public:
    int maxScore(vector<int>& nums) {
        int N = nums.size();
        vector<bool> visit(N, false);
        return dfs(nums, visit, 1, N);
    }

private:
    int dfs(vector<int>& nums, vector<bool>& visit, int n, int N) {
        if (n > N / 2) {
            return 0;
        }

        int res = 0;
        for (int i = 0; i < N; i++) {
            if (visit[i]) continue;
            visit[i] = true;
            for (int j = i + 1; j < N; j++) {
                if (visit[j]) continue;
                visit[j] = true;
                int g = gcd(nums[i], nums[j]);
                res = max(res, n * g + dfs(nums, visit, n + 1, N));
                visit[j] = false;
            }
            visit[i] = false;
        }

        return res;
    }

    int gcd(int a, int b) {
        return b == 0 ? a : gcd(b, a % b);
    }
};

class Solution {
    /**
     * @param {number[]} nums
     * @return {number}
     */
    maxScore(nums) {
        const N = nums.length;
        const visit = new Array(N).fill(false);

        const gcd = (a, b) => {
            return b === 0 ? a : gcd(b, a % b);
        };
        const dfs = (n) => {
            if (n > N / 2) {
                return 0;
            }

            let res = 0;
            for (let i = 0; i < N; i++) {
                if (visit[i]) continue;
                visit[i] = true;
                for (let j = i + 1; j < N; j++) {
                    if (visit[j]) continue;
                    visit[j] = true;
                    let g = gcd(nums[i], nums[j]);
                    res = Math.max(res, n * g + dfs(n + 1));
                    visit[j] = false;
                }
                visit[i] = false;
            }

            return res;
        };

        return dfs(1);
    }
}

public class Solution {
    public int MaxScore(int[] nums) {
        int N = nums.Length;
        bool[] visit = new bool[N];
        return Dfs(nums, visit, 1, N);
    }

    private int Dfs(int[] nums, bool[] visit, int n, int N) {
        if (n > N / 2) {
            return 0;
        }

        int res = 0;
        for (int i = 0; i < N; i++) {
            if (visit[i]) continue;
            visit[i] = true;
            for (int j = i + 1; j < N; j++) {
                if (visit[j]) continue;
                visit[j] = true;
                int g = Gcd(nums[i], nums[j]);
                res = Math.Max(res, n * g + Dfs(nums, visit, n + 1, N));
                visit[j] = false;
            }
            visit[i] = false;
        }

        return res;
    }

    private int Gcd(int a, int b) {
        return b == 0 ? a : Gcd(b, a % b);
    }
}

func maxScore(nums []int) int {
    N := len(nums)
    visit := make([]bool, N)

    var gcd func(a, b int) int
    gcd = func(a, b int) int {
        if b == 0 {
            return a
        }
        return gcd(b, a%b)
    }

    var dfs func(n int) int
    dfs = func(n int) int {
        if n > N/2 {
            return 0
        }

        res := 0
        for i := 0; i < N; i++ {
            if visit[i] {
                continue
            }
            visit[i] = true
            for j := i + 1; j < N; j++ {
                if visit[j] {
                    continue
                }
                visit[j] = true
                g := gcd(nums[i], nums[j])
                score := n*g + dfs(n+1)
                if score > res {
                    res = score
                }
                visit[j] = false
            }
            visit[i] = false
        }

        return res
    }

    return dfs(1)
}

class Solution {
    fun maxScore(nums: IntArray): Int {
        val N = nums.size
        val visit = BooleanArray(N)

        fun gcd(a: Int, b: Int): Int = if (b == 0) a else gcd(b, a % b)

        fun dfs(n: Int): Int {
            if (n > N / 2) return 0

            var res = 0
            for (i in 0 until N) {
                if (visit[i]) continue
                visit[i] = true
                for (j in i + 1 until N) {
                    if (visit[j]) continue
                    visit[j] = true
                    val g = gcd(nums[i], nums[j])
                    res = maxOf(res, n * g + dfs(n + 1))
                    visit[j] = false
                }
                visit[i] = false
            }

            return res
        }

        return dfs(1)
    }
}

class Solution {
    func maxScore(_ nums: [Int]) -> Int {
        let N = nums.count
        var visit = [Bool](repeating: false, count: N)

        func gcd(_ a: Int, _ b: Int) -> Int {
            return b == 0 ? a : gcd(b, a % b)
        }

        func dfs(_ n: Int) -> Int {
            if n > N / 2 {
                return 0
            }

            var res = 0
            for i in 0..<N {
                if visit[i] { continue }
                visit[i] = true
                for j in (i + 1)..<N {
                    if visit[j] { continue }
                    visit[j] = true
                    let g = gcd(nums[i], nums[j])
                    res = max(res, n * g + dfs(n + 1))
                    visit[j] = false
                }
                visit[i] = false
            }

            return res
        }

        return dfs(1)
    }
}

Time & Space Complexity

Time complexity: $O(n ^ n * \log m)$
Space complexity: $O(n)$

Where $n$ is the size of the array $nums$ and $m$ is the maximum element in the array.

2. Bitmask DP (Top-Down) - I

Intuition

The brute force approach has redundant computation because the same subset of remaining elements can be reached through different pairing sequences. We can optimize by caching results based on which elements have been used.

A bitmask elegantly represents which elements are taken. Each bit position corresponds to an array index: bit i being 1 means nums[i] is used. Since the operation number can be derived from counting set bits (every 2 used elements means one completed operation), the mask alone uniquely defines the state.

Algorithm

Use a hash map cache to store the maximum score achievable from each bitmask state.
Define dfs(mask, op) where mask tracks used elements and op is the current operation number.
If mask is already cached, return the stored value.
For each pair of indices (i, j) where neither bit is set:
- Create newMask by setting bits i and j.
- Compute op * gcd(nums[i], nums[j]) plus recursive result.
- Update the maximum score for this mask.
Cache and return the result.

class Solution:
    def maxScore(self, nums: List[int]) -> int:
        cache = collections.defaultdict(int)

        def dfs(mask, op):
            if mask in cache:
                return cache[mask]

            for i in range(len(nums)):
                for j in range(i + 1, len(nums)):
                    if (1 << i) & mask or (1 << j) & mask:
                        continue

                    newMask = mask | (1 << i) | (1 << j)
                    score = op * math.gcd(nums[i], nums[j])
                    cache[mask] = max(
                        cache[mask],
                        score + dfs(newMask, op + 1)
                    )

            return cache[mask]

        return dfs(0, 1)

public class Solution {
    private Map<Integer, Integer> cache;

    public int maxScore(int[] nums) {
        cache = new HashMap<>();
        return dfs(0, 1, nums);
    }

    private int dfs(int mask, int op, int[] nums) {
        if (cache.containsKey(mask)) {
            return cache.get(mask);
        }

        int maxScore = 0;
        int n = nums.length;
        for (int i = 0; i < n; i++) {
            if ((mask & (1 << i)) != 0) continue;
            for (int j = i + 1; j < n; j++) {
                if ((mask & (1 << j)) != 0) continue;
                int newMask = mask | (1 << i) | (1 << j);
                int score = op * gcd(nums[i], nums[j]) + dfs(newMask, op + 1, nums);
                maxScore = Math.max(maxScore, score);
            }
        }

        cache.put(mask, maxScore);
        return maxScore;
    }

    private int gcd(int a, int b) {
        return b == 0 ? a : gcd(b, a % b);
    }
}

class Solution {
public:
    int maxScore(vector<int>& nums) {
        return dfs(0, 1, nums);
    }

private:
    unordered_map<int, int> cache;

    int dfs(int mask, int op, vector<int>& nums) {
        if (cache.count(mask)) {
            return cache[mask];
        }

        int maxScore = 0;
        int n = nums.size();
        for (int i = 0; i < n; i++) {
            if ((mask & (1 << i)) != 0) continue;
            for (int j = i + 1; j < n; j++) {
                if ((mask & (1 << j)) != 0) continue;
                int newMask = mask | (1 << i) | (1 << j);
                int score = op * gcd(nums[i], nums[j]) + dfs(newMask, op + 1, nums);
                maxScore = max(maxScore, score);
            }
        }

        return cache[mask] = maxScore;
    }
};

class Solution {
    /**
     * @param {number[]} nums
     * @return {number}
     */
    maxScore(nums) {
        const cache = new Map();

        const gcd = (a, b) => (b === 0 ? a : gcd(b, a % b));

        const dfs = (mask, op) => {
            if (cache.has(mask)) {
                return cache.get(mask);
            }

            let maxScore = 0;
            const n = nums.length;
            for (let i = 0; i < n; i++) {
                if ((mask & (1 << i)) !== 0) continue;
                for (let j = i + 1; j < n; j++) {
                    if ((mask & (1 << j)) !== 0) continue;
                    let newMask = mask | (1 << i) | (1 << j);
                    let score =
                        op * gcd(nums[i], nums[j]) + dfs(newMask, op + 1);
                    maxScore = Math.max(maxScore, score);
                }
            }

            cache.set(mask, maxScore);
            return maxScore;
        };

        return dfs(0, 1);
    }
}

public class Solution {
    private Dictionary<int, int> cache;

    public int MaxScore(int[] nums) {
        cache = new Dictionary<int, int>();
        return Dfs(0, 1, nums);
    }

    private int Dfs(int mask, int op, int[] nums) {
        if (cache.ContainsKey(mask)) {
            return cache[mask];
        }

        int maxScore = 0;
        int n = nums.Length;
        for (int i = 0; i < n; i++) {
            if ((mask & (1 << i)) != 0) continue;
            for (int j = i + 1; j < n; j++) {
                if ((mask & (1 << j)) != 0) continue;
                int newMask = mask | (1 << i) | (1 << j);
                int score = op * Gcd(nums[i], nums[j]) + Dfs(newMask, op + 1, nums);
                maxScore = Math.Max(maxScore, score);
            }
        }

        cache[mask] = maxScore;
        return maxScore;
    }

    private int Gcd(int a, int b) {
        return b == 0 ? a : Gcd(b, a % b);
    }
}

func maxScore(nums []int) int {
    cache := make(map[int]int)

    var gcd func(a, b int) int
    gcd = func(a, b int) int {
        if b == 0 {
            return a
        }
        return gcd(b, a%b)
    }

    var dfs func(mask, op int) int
    dfs = func(mask, op int) int {
        if val, ok := cache[mask]; ok {
            return val
        }

        maxScore := 0
        n := len(nums)
        for i := 0; i < n; i++ {
            if (mask & (1 << i)) != 0 {
                continue
            }
            for j := i + 1; j < n; j++ {
                if (mask & (1 << j)) != 0 {
                    continue
                }
                newMask := mask | (1 << i) | (1 << j)
                score := op*gcd(nums[i], nums[j]) + dfs(newMask, op+1)
                if score > maxScore {
                    maxScore = score
                }
            }
        }

        cache[mask] = maxScore
        return maxScore
    }

    return dfs(0, 1)
}

class Solution {
    private val cache = HashMap<Int, Int>()

    fun maxScore(nums: IntArray): Int {
        return dfs(0, 1, nums)
    }

    private fun dfs(mask: Int, op: Int, nums: IntArray): Int {
        if (cache.containsKey(mask)) {
            return cache[mask]!!
        }

        var maxScore = 0
        val n = nums.size
        for (i in 0 until n) {
            if ((mask and (1 shl i)) != 0) continue
            for (j in i + 1 until n) {
                if ((mask and (1 shl j)) != 0) continue
                val newMask = mask or (1 shl i) or (1 shl j)
                val score = op * gcd(nums[i], nums[j]) + dfs(newMask, op + 1, nums)
                maxScore = maxOf(maxScore, score)
            }
        }

        cache[mask] = maxScore
        return maxScore
    }

    private fun gcd(a: Int, b: Int): Int = if (b == 0) a else gcd(b, a % b)
}

class Solution {
    private var cache = [Int: Int]()

    func maxScore(_ nums: [Int]) -> Int {
        cache = [:]
        return dfs(0, 1, nums)
    }

    private func dfs(_ mask: Int, _ op: Int, _ nums: [Int]) -> Int {
        if let val = cache[mask] {
            return val
        }

        var maxScore = 0
        let n = nums.count
        for i in 0..<n {
            if (mask & (1 << i)) != 0 { continue }
            for j in (i + 1)..<n {
                if (mask & (1 << j)) != 0 { continue }
                let newMask = mask | (1 << i) | (1 << j)
                let score = op * gcd(nums[i], nums[j]) + dfs(newMask, op + 1, nums)
                maxScore = max(maxScore, score)
            }
        }

        cache[mask] = maxScore
        return maxScore
    }

    private func gcd(_ a: Int, _ b: Int) -> Int {
        return b == 0 ? a : gcd(b, a % b)
    }
}

Time & Space Complexity

Time complexity: $O(n ^ 2 * 2 ^ n * \log m)$
Space complexity: $O(2 ^ n)$

Where $n$ is the size of the array $nums$ and $m$ is the maximum element in the array.

3. Bitmask DP (Top-Down) - II

Intuition

We can further optimize by precomputing all pairwise GCD values. In the previous approach, we recompute gcd(nums[i], nums[j]) every time we consider that pair. Since the same pair appears in many different states, precomputing these values into a 2D table saves repeated GCD calculations.

Additionally, using an array instead of a hash map for memoization provides faster access since bitmask states are contiguous integers from 0 to 2^n - 1.

Algorithm

Precompute a GCD[i][j] table storing gcd(nums[i], nums[j]) for all pairs where i < j.
Create a dp array of size 2^n initialized to -1 (unvisited).
Define dfs(mask, op) similar to before, but look up GCD values from the precomputed table.
Use the dp array for memoization instead of a hash map.
Return dfs(0, 1) as the answer.

class Solution:
    def maxScore(self, nums: List[int]) -> int:
        n = len(nums)
        GCD = [[0] * n for _ in range(n)]
        for i in range(n):
            for j in range(i + 1, n):
                GCD[i][j] = gcd(nums[i], nums[j])

        dp = [-1] * (1 << n)
        def dfs(mask, op):
            if dp[mask] != -1:
                return dp[mask]

            max_score = 0
            for i in range(n):
                if mask & (1 << i):
                    continue
                for j in range(i + 1, n):
                    if mask & (1 << j):
                        continue
                    new_mask = mask | (1 << i) | (1 << j)
                    max_score = max(
                        max_score,
                        op * GCD[i][j] + dfs(new_mask, op + 1)
                    )

            dp[mask] = max_score
            return max_score

        return dfs(0, 1)

public class Solution {
    private int[][] GCD;
    private int[] dp;

    public int maxScore(int[] nums) {
        int n = nums.length;
        GCD = new int[n][n];
        dp = new int[1 << n];
        Arrays.fill(dp, -1);

        for (int i = 0; i < n; i++) {
            for (int j = i + 1; j < n; j++) {
                GCD[i][j] = gcd(nums[i], nums[j]);
            }
        }

        return (int) dfs(0, 1, nums);
    }

    private int dfs(int mask, int op, int[] nums) {
        if (dp[mask] != -1) return dp[mask];

        int maxScore = 0;
        for (int i = 0; i < nums.length; i++) {
            if ((mask & (1 << i)) != 0) continue;
            for (int j = i + 1; j < nums.length; j++) {
                if ((mask & (1 << j)) != 0) continue;
                int newMask = mask | (1 << i) | (1 << j);
                maxScore = Math.max(
                    maxScore,
                    op * GCD[i][j] + dfs(newMask, op + 1, nums)
                );
            }
        }
        return dp[mask] = maxScore;
    }

    private int gcd(int a, int b) {
        return b == 0 ? a : gcd(b, a % b);
    }
}

class Solution {
    /**
     * @param {number[]} nums
     * @return {number}
     */
    maxScore(nums) {
        const n = nums.length;
        const GCD = Array.from({ length: n }, () => Array(n).fill(0));
        const dp = Array(1 << n).fill(-1);

        const gcd = (a, b) => (b === 0 ? a : gcd(b, a % b));
        for (let i = 0; i < n; i++) {
            for (let j = i + 1; j < n; j++) {
                GCD[i][j] = gcd(nums[i], nums[j]);
            }
        }

        const dfs = (mask, op) => {
            if (dp[mask] !== -1) return dp[mask];

            let maxScore = 0;
            for (let i = 0; i < n; i++) {
                if (mask & (1 << i)) continue;
                for (let j = i + 1; j < n; j++) {
                    if (mask & (1 << j)) continue;
                    const newMask = mask | (1 << i) | (1 << j);
                    maxScore = Math.max(
                        maxScore,
                        op * GCD[i][j] + dfs(newMask, op + 1),
                    );
                }
            }
            return (dp[mask] = maxScore);
        };

        return dfs(0, 1);
    }
}

public class Solution {
    private int[][] GCD;
    private int[] dp;

    public int MaxScore(int[] nums) {
        int n = nums.Length;
        GCD = new int[n][];
        dp = new int[1 << n];
        Array.Fill(dp, -1);

        for (int i = 0; i < n; i++) {
            GCD[i] = new int[n];
            for (int j = i + 1; j < n; j++) {
                GCD[i][j] = Gcd(nums[i], nums[j]);
            }
        }

        return Dfs(0, 1, nums);
    }

    private int Dfs(int mask, int op, int[] nums) {
        if (dp[mask] != -1) return dp[mask];

        int maxScore = 0;
        for (int i = 0; i < nums.Length; i++) {
            if ((mask & (1 << i)) != 0) continue;
            for (int j = i + 1; j < nums.Length; j++) {
                if ((mask & (1 << j)) != 0) continue;
                int newMask = mask | (1 << i) | (1 << j);
                maxScore = Math.Max(
                    maxScore,
                    op * GCD[i][j] + Dfs(newMask, op + 1, nums)
                );
            }
        }
        return dp[mask] = maxScore;
    }

    private int Gcd(int a, int b) {
        return b == 0 ? a : Gcd(b, a % b);
    }
}

class Solution {
    private lateinit var GCD: Array<IntArray>
    private lateinit var dp: IntArray

    fun maxScore(nums: IntArray): Int {
        val n = nums.size
        GCD = Array(n) { IntArray(n) }
        dp = IntArray(1 shl n) { -1 }

        for (i in 0 until n) {
            for (j in i + 1 until n) {
                GCD[i][j] = gcd(nums[i], nums[j])
            }
        }

        return dfs(0, 1, nums)
    }

    private fun dfs(mask: Int, op: Int, nums: IntArray): Int {
        if (dp[mask] != -1) return dp[mask]

        var maxScore = 0
        for (i in nums.indices) {
            if ((mask and (1 shl i)) != 0) continue
            for (j in i + 1 until nums.size) {
                if ((mask and (1 shl j)) != 0) continue
                val newMask = mask or (1 shl i) or (1 shl j)
                maxScore = maxOf(
                    maxScore,
                    op * GCD[i][j] + dfs(newMask, op + 1, nums)
                )
            }
        }
        dp[mask] = maxScore
        return maxScore
    }

    private fun gcd(a: Int, b: Int): Int = if (b == 0) a else gcd(b, a % b)
}

class Solution {
    private var GCD: [[Int]] = []
    private var dp: [Int] = []

    func maxScore(_ nums: [Int]) -> Int {
        let n = nums.count
        GCD = Array(repeating: Array(repeating: 0, count: n), count: n)
        dp = Array(repeating: -1, count: 1 << n)

        for i in 0..<n {
            for j in (i + 1)..<n {
                GCD[i][j] = gcd(nums[i], nums[j])
            }
        }

        return dfs(0, 1, nums)
    }

    private func dfs(_ mask: Int, _ op: Int, _ nums: [Int]) -> Int {
        if dp[mask] != -1 { return dp[mask] }

        var maxScore = 0
        for i in 0..<nums.count {
            if (mask & (1 << i)) != 0 { continue }
            for j in (i + 1)..<nums.count {
                if (mask & (1 << j)) != 0 { continue }
                let newMask = mask | (1 << i) | (1 << j)
                maxScore = max(
                    maxScore,
                    op * GCD[i][j] + dfs(newMask, op + 1, nums)
                )
            }
        }
        dp[mask] = maxScore
        return maxScore
    }

    private func gcd(_ a: Int, _ b: Int) -> Int {
        return b == 0 ? a : gcd(b, a % b)
    }
}

Time & Space Complexity

Time complexity: $O(n ^ 2 * (2 ^ n + \log m))$
Space complexity: $O(n ^ 2 + 2 ^ n)$

Where $n$ is the size of the array $nums$ and $m$ is the maximum element in the array.

4. Bitmask DP (Bottom-Up)

Intuition

Instead of recursion with memoization, we can fill the DP table iteratively. The key observation is that a state with more bits set depends on states with fewer bits set. By iterating from higher masks (more elements used) to lower masks (fewer elements used), we ensure that when we process a state, all states it transitions to are already computed.

The operation number for a given mask is determined by counting how many bits are set and dividing by 2, then adding 1 for the next operation.

Algorithm

Precompute the GCD[i][j] table for all pairs.
Create a dp array of size 2^n initialized to 0.
Iterate mask from 2^n - 1 down to 0:
- Count the number of set bits. Skip if it is odd (invalid state).
- Calculate op = bits / 2 + 1 (the next operation number).
- For each pair (i, j) not in the current mask:
  - Compute new_mask = mask | (1 << i) | (1 << j).
  - Update dp[mask] = max(dp[mask], op * GCD[i][j] + dp[new_mask]).
Return dp[0].

class Solution:
    def maxScore(self, nums: List[int]) -> int:
        n = len(nums)
        N = 1 << n
        GCD = [[0] * n for _ in range(n)]
        for i in range(n):
            for j in range(i + 1, n):
                GCD[i][j] = gcd(nums[i], nums[j])

        dp = [0] * N
        for mask in range(N - 1, -1, -1):
            bits = bin(mask).count('1')
            if bits % 2 == 1:
                continue
            op = bits // 2 + 1

            for i in range(n):
                if mask & (1 << i):
                    continue
                for j in range(i + 1, n):
                    if mask & (1 << j):
                        continue
                    new_mask = mask | (1 << i) | (1 << j)
                    dp[mask] = max(dp[mask], op * GCD[i][j] + dp[new_mask])

        return dp[0]

class Solution {
public:
    int maxScore(vector<int>& nums) {
        int n = nums.size();
        int N = 1 << n;
        vector<vector<int>> GCD(n, vector<int>(n, 0));

        for (int i = 0; i < n; i++) {
            for (int j = i + 1; j < n; j++) {
                GCD[i][j] = __gcd(nums[i], nums[j]);
            }
        }

        vector<int> dp(N, 0);
        for (int mask = N - 1; mask >= 0; mask--) {
            int bits = __builtin_popcount(mask);
            if (bits % 2 == 1) continue;
            int op = bits / 2 + 1;

            for (int i = 0; i < n; i++) {
                if (mask & (1 << i)) continue;
                for (int j = i + 1; j < n; j++) {
                    if (mask & (1 << j)) continue;
                    int newMask = mask | (1 << i) | (1 << j);
                    dp[mask] = max(dp[mask], op * GCD[i][j] + dp[newMask]);
                }
            }
        }
        return dp[0];
    }
};

class Solution {
    /**
     * @param {number[]} nums
     * @return {number}
     */
    maxScore(nums) {
        const n = nums.length;
        const N = 1 << n;
        const GCD = Array.from({ length: n }, () => Array(n).fill(0));

        const gcd = (a, b) => (b === 0 ? a : gcd(b, a % b));
        for (let i = 0; i < n; i++) {
            for (let j = i + 1; j < n; j++) {
                GCD[i][j] = gcd(nums[i], nums[j]);
            }
        }

        const dp = Array(N).fill(0);
        for (let mask = N - 1; mask >= 0; mask--) {
            let bits = mask.toString(2).split('0').join('').length;
            if (bits % 2 === 1) continue;
            let op = bits / 2 + 1;

            for (let i = 0; i < n; i++) {
                if ((mask & (1 << i)) !== 0) continue;
                for (let j = i + 1; j < n; j++) {
                    if ((mask & (1 << j)) !== 0) continue;
                    let newMask = mask | (1 << i) | (1 << j);
                    dp[mask] = Math.max(dp[mask], op * GCD[i][j] + dp[newMask]);
                }
            }
        }
        return dp[0];
    }
}

public class Solution {
    public int MaxScore(int[] nums) {
        int n = nums.Length;
        int N = 1 << n;
        int[][] GCD = new int[n][];

        for (int i = 0; i < n; i++) {
            GCD[i] = new int[n];
            for (int j = i + 1; j < n; j++) {
                GCD[i][j] = Gcd(nums[i], nums[j]);
            }
        }

        int[] dp = new int[N];
        for (int mask = N - 1; mask >= 0; mask--) {
            int bits = BitCount(mask);
            if (bits % 2 == 1) continue;
            int op = bits / 2 + 1;

            for (int i = 0; i < n; i++) {
                if ((mask & (1 << i)) != 0) continue;
                for (int j = i + 1; j < n; j++) {
                    if ((mask & (1 << j)) != 0) continue;
                    int newMask = mask | (1 << i) | (1 << j);
                    dp[mask] = Math.Max(dp[mask], op * GCD[i][j] + dp[newMask]);
                }
            }
        }
        return dp[0];
    }

    private int Gcd(int a, int b) {
        return b == 0 ? a : Gcd(b, a % b);
    }

    private int BitCount(int n) {
        int count = 0;
        while (n != 0) {
            count += n & 1;
            n >>= 1;
        }
        return count;
    }
}

func maxScore(nums []int) int {
    n := len(nums)
    N := 1 << n
    GCD := make([][]int, n)
    for i := range GCD {
        GCD[i] = make([]int, n)
    }

    var gcd func(a, b int) int
    gcd = func(a, b int) int {
        if b == 0 {
            return a
        }
        return gcd(b, a%b)
    }

    for i := 0; i < n; i++ {
        for j := i + 1; j < n; j++ {
            GCD[i][j] = gcd(nums[i], nums[j])
        }
    }

    dp := make([]int, N)
    for mask := N - 1; mask >= 0; mask-- {
        bits := popcount(mask)
        if bits%2 == 1 {
            continue
        }
        op := bits/2 + 1

        for i := 0; i < n; i++ {
            if (mask & (1 << i)) != 0 {
                continue
            }
            for j := i + 1; j < n; j++ {
                if (mask & (1 << j)) != 0 {
                    continue
                }
                newMask := mask | (1 << i) | (1 << j)
                score := op*GCD[i][j] + dp[newMask]
                if score > dp[mask] {
                    dp[mask] = score
                }
            }
        }
    }
    return dp[0]
}

func popcount(x int) int {
    count := 0
    for x != 0 {
        count += x & 1
        x >>= 1
    }
    return count
}

class Solution {
    fun maxScore(nums: IntArray): Int {
        val n = nums.size
        val N = 1 shl n
        val GCD = Array(n) { IntArray(n) }

        for (i in 0 until n) {
            for (j in i + 1 until n) {
                GCD[i][j] = gcd(nums[i], nums[j])
            }
        }

        val dp = IntArray(N)
        for (mask in N - 1 downTo 0) {
            val bits = Integer.bitCount(mask)
            if (bits % 2 == 1) continue
            val op = bits / 2 + 1

            for (i in 0 until n) {
                if ((mask and (1 shl i)) != 0) continue
                for (j in i + 1 until n) {
                    if ((mask and (1 shl j)) != 0) continue
                    val newMask = mask or (1 shl i) or (1 shl j)
                    dp[mask] = maxOf(dp[mask], op * GCD[i][j] + dp[newMask])
                }
            }
        }
        return dp[0]
    }

    private fun gcd(a: Int, b: Int): Int = if (b == 0) a else gcd(b, a % b)
}

class Solution {
    func maxScore(_ nums: [Int]) -> Int {
        let n = nums.count
        let N = 1 << n
        var GCD = Array(repeating: Array(repeating: 0, count: n), count: n)

        func gcd(_ a: Int, _ b: Int) -> Int {
            return b == 0 ? a : gcd(b, a % b)
        }

        for i in 0..<n {
            for j in (i + 1)..<n {
                GCD[i][j] = gcd(nums[i], nums[j])
            }
        }

        var dp = Array(repeating: 0, count: N)
        for mask in stride(from: N - 1, through: 0, by: -1) {
            let bits = mask.nonzeroBitCount
            if bits % 2 == 1 { continue }
            let op = bits / 2 + 1

            for i in 0..<n {
                if (mask & (1 << i)) != 0 { continue }
                for j in (i + 1)..<n {
                    if (mask & (1 << j)) != 0 { continue }
                    let newMask = mask | (1 << i) | (1 << j)
                    dp[mask] = max(dp[mask], op * GCD[i][j] + dp[newMask])
                }
            }
        }
        return dp[0]
    }
}

Time & Space Complexity

Time complexity: $O(n ^ 2 * (2 ^ n + \log m))$
Space complexity: $O(n ^ 2 + 2 ^ n)$

Where $n$ is the size of the array $nums$ and $m$ is the maximum element in the array.

Common Pitfalls

Forgetting to Multiply by Operation Number

Each operation's score is i * gcd(nums[a], nums[b]) where i is the operation number (1-indexed). Forgetting the multiplier or using 0-indexed operation numbers produces incorrect scores.

Using Wrong Bit Count for Operation Number

In bitmask DP, the operation number is derived from the number of set bits divided by 2, plus 1 (since each operation uses 2 elements). Off-by-one errors in this calculation cascade through all subsequent operations.

Not Precomputing GCD Values

Recomputing gcd(nums[i], nums[j]) inside the DP loop for every state transition leads to TLE on larger inputs. Precomputing all pairwise GCD values into a 2D table provides significant speedup.

Skipping Invalid Bitmask States

In bottom-up DP, states with an odd number of set bits are invalid (you cannot pair up an odd number of elements). Failing to skip these states wastes computation and can cause incorrect transitions.

Incorrect Bitmask Transition

When marking elements as used, the new mask should be mask | (1 << i) | (1 << j). Using XOR or other operations instead of OR can incorrectly toggle bits that are already set, corrupting the state.