Codeforces 453B B. Little Pony and Harmony Chest(dp+数论)

题解 ： 首先发现这个题目中的 b 数组中每一个元素都不会超过 58 为什么呢 ？ 自己想去 ….
然后我们发现，两个数互质就是没有相同的素因子，然后 1 －57 只有 18 个素数，先将这 18 个数打进一个表然后，用
state[i] 存下每个数包含的所有的素因子。
剩下的就是转移了
dp[i][j] 表示前 i 个 用了 j 这些素数 (j 是一个二维的压缩状态)。
dp[i + 1][j | state[k]] = min (dp[i][j] + abs (k - a[i]),dp[i + 1][j | k]) ;
由于这个题目要求保留路径，那么就还需要开一个数组去储存路径就好了

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
using namespace std;

const int INF = 1e9 + 7;
const int maxn = (1 << 18);
int prime[] = {2,3,5,7,11,13,17,19,23,29,31,37,39,41,43,47,53,57};
int a[105] = {0};
int b[105] = {0};
int state[60] = {0};
int dp[105][maxn] = {0};
int pos[105][maxn] = {0};
int main() {
    ios_base :: sync_with_stdio(false);
    int n;
    cin >> n;
    for (int i = 1;i <= n; ++ i) cin >> a[i];
    for (int i = 2;i < 60; ++ i) {
        for (int j = 0;j < 18; ++ j) {
            if (i % prime[j] == 0) {
                state[i] |= (1 << j);
            }
        }
    }
    for (int i = 0;i < 105; ++ i) {
        for (int j = 0; j < maxn; ++ j) dp[i][j] = INF;
    }
    dp[0][0] = 0;
    for (int i = 0;i < n; ++ i) {
        for (int j = 0;j < maxn; ++ j) {
            if (dp[i][j] == INF) continue;
            if (dp[i + 1][j] > dp[i][j] + a[i + 1] - 1) {
                dp[i + 1][j] = dp[i][j] + a[i + 1] - 1;
                pos[i + 1][j] = 1;
            }
            for (int k = 2;k < 60; ++ k) {
                if (state[k] & j) continue;
                else {
                    int u = abs (a[i + 1] - k);
                    if (u + dp[i][j] < dp[i + 1][j | state[k]]) {
                        dp[i + 1][j | state[k]] = dp[i][j] + u;
                        pos[i + 1][j | state[k]] = k;
                    }
                }
            }
        }
    }
    int ans = INF;
    int tem = 0;
    for (int i = 0;i < maxn; ++ i) {
        if (dp[n][i] < ans) {
            ans = dp[n][i];
            tem = i;
        }
    }
    int cnt = 0;
    for (int i = n;i > 0;-- i) {
        b[cnt ++] = pos[i][tem];
        tem -= state[pos[i][tem]];
    }
    for (int i = cnt - 1;i >= 0 ; -- i) {
        cout << b[i] << ' ';
    }
    return 0;
}