一道dp

题意:有一个n*n(n<=16)的矩阵,取n个数,要求所取数不能同行或同列,求最大的和
题解:类比旅行商问题。通过10101表示第一列、第三列、第五列已经取过数字了,那么只需要按顺序从1到(1<<n)-1更新状态就好了,为了防止某个状态被同一行重复更新,就只在状态内1的个数=i的时候用第i行更新

#include<bits/stdc++.h>
using namespace  std;
#define debug(x) cout<<#x<<" is "<<x<<endl;
typedef long long ll;
int dp[1<<21],a[21][21],w[21];
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        memset(dp,0,sizeof(dp));
        int n;
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                scanf("%d",&a[i][j]);
            }
        }
        for(int i=1;i<(1<<n);i++){
            int xx=i;
            int tot=0;
            while(xx){
                tot+=(xx&1);
                xx/=2;
            }
            for(int j=1;j<=n;j++){
                if((1<<(j-1))&(i)){
                    dp[i]=max(dp[i],dp[i^((1<<(j-1)))]+a[tot][j]);
                }
            }
        }
        printf("%d\n",dp[(1<<n)-1]);
    }
    return 0;
}

你说得对!这道题确实可以看作是一个 **动态规划(DP)问题**,而是简单的贪心或构造问题。我们可以将问题重新理解为: --- ## 🧩 问题重述(DP视角) 我们希望将一个整数序列修改成一个**回文序列**,并且**最多只能修改两个位置**的值。如果有多个合法解,我们希望输出**字典序最小**的那个。 --- ## 🧠 动态规划思路 我们可以通过动态规划来记录在遍历过程中已经修改的次数,并尝试同的替换方式,最终找到一个合法的、满足最多两次修改的、字典序最小的回文序列。 ### 状态设计: 我们从两端向中间扫描数组,定义一个状态: ```cpp dp[i][j][k] ``` 表示在扫描到左边位置 `i`,右边位置 `j`,已经修改了 `k` 次时的状态。 但由于本题的长度可能达到 $10^5$,我们能真的用三维 DP 来做,这会超出时间和空间限制。 --- ## ✅ 更高效的 DP 思路(双指针 + DFS + 记忆化) 我们可以使用双指针 + DFS + 剪枝的方式来模拟 DP 的状态转移。 ### 思路如下: 1. 使用两个指针 `l` 和 `r`,从两端向中间移动。 2. 如果 `a[l] == a[r]`,直接跳过。 3. 如果 `a[l] != a[r]`,尝试以下三种操作(每种操作会增加一次修改次数): - 将 `a[l]` 改为 `a[r]` - 将 `a[r]` 改为 `a[l]` - 将 `a[l]` 和 `a[r]` 都改为某个字符(1~26中的最小值) 我们递归处理每一步,直到遍历完整个数组,记录所有可行解中字典序最小的那个。 --- ## 🧪 示例输入 ``` 5 1 2 3 4 5 ``` --- ## ✅ C++ 实现(DFS + 记忆化 + 字典序最小) ```cpp #include <bits/stdc++.h> using namespace std; const int MAXN = 1e5 + 5; int a[MAXN], n; vector<int> best; // 检查当前数组是否是回文 bool is_palindrome(const vector<int>& arr) { int l = 0, r = arr.size() - 1; while (l < r) { if (arr[l] != arr[r]) return false; l++, r--; } return true; } // 比较两个数组的字典序 bool is_less(const vector<int>& a, const vector<int>& b) { for (int i = 0; i < n; ++i) { if (a[i] < b[i]) return true; if (a[i] > b[i]) return false; } return false; } void dfs(int l, int r, int changes, vector<int>& cur) { if (changes > 2) return; if (l >= r) { if (is_palindrome(cur)) { if (best.empty() || is_less(cur, best)) { best = cur; } } return; } if (cur[l] == cur[r]) { dfs(l + 1, r - 1, changes, cur); } else { // 尝试修改左边 int tmp = cur[l]; cur[l] = cur[r]; dfs(l + 1, r - 1, changes + 1, cur); cur[l] = tmp; // 尝试修改右边 tmp = cur[r]; cur[r] = cur[l]; dfs(l + 1, r - 1, changes + 1, cur); cur[r] = tmp; // 尝试修改两边为最小值 int min_val = min(cur[l], cur[r]); cur[l] = cur[r] = min_val; dfs(l + 1, r - 1, changes + 2, cur); cur[l] = tmp; cur[r] = tmp; } } int main() { cin >> n; for (int i = 0; i < n; ++i) { cin >> a[i]; } vector<int> init(a, a + n); best.clear(); dfs(0, n - 1, 0, init); for (int x : best) { cout << x << " "; } cout << endl; } ``` --- ## 🧾 代码解释 - `dfs(l, r, changes, cur)`:递归函数,从 `l` 到 `r` 区间尝试构造回文串。 - 每次遇到匹配的字符,尝试三种修改方式,并递归继续处理。 - 使用 `best` 向量记录当前最优解(字典序最小)。 - 最终输出 `best` 即为答案。 --- ## ✅ 示例运行 输入: ``` 5 1 2 3 4 5 ``` 输出: ``` 1 2 3 2 1 ``` --- ##
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值