Codeforces 1242C Sum Balance(状压DP)

题目链接
题意

给你 n n n 行每行若干整数(所有整数互不相同),每行拿出任意一个整数,再给每行放回一个拿出的整数,求这样操作一次后所有行的整数和相同的操作方案

思路

第一次打DIV1,没想到代码能力太弱想到思路却没敲出来,题还是刷太少。

先验证所有整数和除以n是否存在余数判断

之后对每个点求其所在的环,每个环长度不超过n,所以这里复杂度是 O ( n ∗ 5000 ∗ n ) O(n*5000*n) O(n5000n)
将每个环上以环包含那几行(二进制)为key,将元素下标存入vector。
之后求一个包含所有行的集合,其中各个环不能相交的方案。这个使用状压dp解决

最后根据 d p [ ( 1 < < n ) − 1 ] dp[(1<<n)-1] dp[(1<<n)1] 判断是否有解,和输出路径。

整体复杂度 O ( n ∗ n ∗ 5000 + 3 n ) O(n*n*5000+3^n) O(nn5000+3n)

代码
#include <bits/stdc++.h>
using namespace std;

#define ll long long
#define _p pair<ll,ll>
#define fi first
#define se second

ll n, sum, up[20], e[20][5005], rsum[20];
_p to[20][5005];
map<ll,_p> mp;
vector<_p> path[1<<15];
ll vis[1<<15], dp[1<<15], pre[1<<15];
ll book[20];
_p ans[20];

int main() {
    scanf("%lld",&n);
    for(ll i = 1; i <= n; ++i) {
        scanf("%lld",&up[i]);
        for(ll j = 0; j < up[i]; ++j) {
            scanf("%lld",&e[i][j]);
            sum += e[i][j];
            rsum[i] += e[i][j];
            mp[e[i][j]] = {i,j};
        }
    }
    if(sum%n) return !printf("No\n");
    // build
    for(ll i = 1; i <= n; ++i) {
        for(ll j = 0; j < up[i]; ++j) {
            to[i][j] = mp.count(sum/n-(rsum[i]-e[i][j])) ? mp[sum/n-(rsum[i]-e[i][j])] : make_pair(0ll,0ll);
        }
    }
    // getcircle
    for(ll i = 1; i <= n; ++i) {
        for(ll j = 0; j < up[i]; ++j) {
            _p now = {i,j};
            vector<_p> tmp;
            memset(book,0,sizeof(book));
            ll sta = 0, f = 0;
            while(1) {
                sta += (1<<(now.fi-1));
                book[now.fi] = 1;
                tmp.push_back(now);
                now = to[now.fi][now.se];
                if(now == make_pair(0ll,0ll)) break;
                if(book[now.fi]) {
                    if(now == make_pair(i,j)) f = 1;
                    break;
                }
            }
            if(f && path[sta].size() == 0) {
                path[sta] = tmp;
                vis[sta] = 1;
            }
        }
    }
    // getans
    dp[0] = 1;
    for(ll i = 1; i < (1<<n); ++i) {
        for(ll j = i; j; j = i&(j-1)) {
            if(vis[j] && dp[i^j]) {
                dp[i] = 1;
                pre[i] = j;
            }
        }
    }
    if(dp[(1<<n)-1] == 0) return !printf("No\n");
    printf("Yes\n");
    ll now = (1<<n)-1;
    while(now) {
        _p la;
        ll f = 0;
//        printf("%lld -\n",pre[now]);
        for(auto v : path[pre[now]]) {
//            printf("%lld -- %lld\n",v.fi,v.se);
            if(f) ans[v.fi] = make_pair(e[v.fi][v.se],la.fi);
            ++f;
            la = v;
        }
        ans[path[pre[now]][0].fi] = make_pair(e[path[pre[now]][0].fi][path[pre[now]][0].se], la.fi);
        now ^= pre[now];
    }
    for(ll i = 1; i <= n; ++i) printf("%lld %lld\n",ans[i].fi,ans[i].se);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
引用\[1\]中提到了一种树形动态规划的方法来解决CodeForces - 982C问题。在这个问题中,subtree指的是子连通块,而不是子树。为了使cnt_white - cnt_black尽可能大,可以使用两次树形动态规划来求解。第一次是自底向上的过程,维护一个dp数组,表示以每个节点为根的子树中的最大连通块。第二次是自顶向下的过程,处理自底向上过程中无法包含的树链所代表的子树。在第二次遍历中,需要维护一个sum变量,用于存储树链所代表的子树的贡献。根据ans\[u\]的正负,决定是否能对相邻的子节点做出贡献。如果ans\[u\]为正,则减去dp\[v\]就是树链所代表的子树的权值。最终,ans\[u\]代表包含节点u在内的子连通块的最大权值。\[1\] 问题: CodeForces - 982C 树形DP是什么问题?如何解决? 回答: CodeForces - 982C是一个树形动态规划问题。在这个问题中,需要求解子连通块的最大权值和,使得cnt_white - cnt_black尽可能大。解决这个问题的方法是使用两次树形动态规划。第一次是自底向上的过程,维护一个dp数组,表示以每个节点为根的子树中的最大连通块。第二次是自顶向下的过程,处理自底向上过程中无法包含的树链所代表的子树。在第二次遍历中,需要维护一个sum变量,用于存储树链所代表的子树的贡献。根据ans\[u\]的正负,决定是否能对相邻的子节点做出贡献。最终,ans\[u\]代表包含节点u在内的子连通块的最大权值。\[1\] #### 引用[.reference_title] - *1* *2* [CodeForces - 1324F Maximum White Subtree(树形dp)](https://blog.csdn.net/qq_45458915/article/details/104831678)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^koosearch_v1,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值