ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

[CF1215E] Marbles - 状压dp

2021-02-12 16:32:24  阅读:150  来源: 互联网

标签:20 Marbles int 状压 珠子 序列 颜色 多少 dp


[CF1215E] Marbles

Description

有 n \((n \le 4 \times 10^5)\) 个珠子,第 \(i\) 个珠子颜色是 \(c_i (c_i \le 20)\),每次操作把相邻的两个珠子交换。现在要把相同颜色的珠子排列在相连的一段,问至少要多少次操作。

Solution

设 \(f[s]\) 表示处理了颜色集合为 s 的珠子(排在目标序列的开头),总操作次数最小是多少

转移无非就是枚举下一次选什么颜色的珠子

代价计算的方法有两种

一种是把目标位置和初始位置求和做差,这样挺麻烦的

另一种是对于已经排上的所有 q,现在新拍上的 p,计算所有 x 为 p 色,y 为 q 色,原始序列中 x 在 y 前面的情况

这可以进一步描述为一个数组 \(inv[i][j]\),表示 i 色的排在 j 色的前面的有多少个,可以 O(nmm) 计算出

具体地,枚举 i,j,然后扫一遍,边扫边统计 i 色已有的个数,然后碰到一个 j 色的就给他加上去

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

#define int long long 

const int N = 400005;
const int M = 25;

int g[M][M], f[N * 6], n, c[N];
//设 $f[s]$ 表示处理了颜色集合为 s 的珠子(排在目标序列的开头),总操作次数最小是多少
// 表示 i 色的排在 j 色的前面的有多少个
//对于已经排上的所有 q,现在新拍上的 p,计算所有 x 为 p 色,y 为 q 色,原始序列中 x 在 y 前面的情况

signed main()
{
    ios::sync_with_stdio(false);
    cin >> n;
    for (int i = 1; i <= n; i++)
    {
        cin >> c[i];
        --c[i];
    }
    for (int i = 0; i < 20; i++)
        for (int j = 0; j < 20; j++)
        {
            int sum = 0;
            for (int k = 1; k <= n; k++)
            {
                if (c[k] == i)
                    ++sum;
                if (c[k] == j)
                    g[i][j] += sum;
            }
        }
    memset(f, 0x3f, sizeof f);
    f[0] = 0;
    for (int s = 0; s < 1 << 20; s++)
    {
        for (int i = 0; i < 20; i++)
            if ((~s) & (1 << i))
            {
                int t = s | (1 << i);
                int delta = 0;
                for (int j = 0; j < 20; j++)
                    if (s & (1 << j))
                    {
                        delta += g[i][j];
                    }
                f[t] = min(f[t], f[s] + delta);
            }
    }

    cout << f[(1 << 20) - 1] << endl;
}

标签:20,Marbles,int,状压,珠子,序列,颜色,多少,dp
来源: https://www.cnblogs.com/mollnn/p/14398930.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有