\(\mathcal{Descriptoin}\)

  Link.

  你有一个容量为 \(k\) 的空书架,现在共有 \(n\) 个请求,每个请求给定一本书 \(a_i\)。如果你的书架里没有这本书,你就必须以 \(c_{a_i}\) 的价格购买这本书放入书架。当然,你可以在任何时候丢掉书架里的某本书。请求出完成这 \(n\) 个请求所需要的最少价钱。

  \(n,k\le80\)。

\(\mathcal{Solution}\)

  网络瘤嘛……

  费用流,考虑先全部买入,再抵消花费。具体地,假设 \(i<j\),第 \(i\) 天的书和第 \(j\) 天的书相同,就可以一直保留第 \(i\) 天的书到第 \(j\) 天,减少一次花费。脑洞一番之后,建图如下:

  • \(S\) 向 \(v_i~(i=1,2,\dots n)\) 连边,容量为 \(1\),费用为 \(c_{a_i}\),表示买入。
  • \(v_i\) 向 \(v_{i+1}~(i=1,2,\cdots,n-1)\) 连边,容量为 \(k-1\),费用为 \(0\),表示保留至多 \(k-1\) 本,剩下一本给 \(a_{i+1}\) 留位置。
  • \(v_i\) 向 \(v_i'~(i=1,2,\cdots,n)\) 连边,容量为 \(1\),费用为 \(0\),表示将这本书出手(丢掉或卖掉)。
  • \(v_{i-1}\) 向上一次 \(a_i\) 出现的位置 \(j\) 所对应的 \(v_j'\) 连边,容量为 \(1\),费用为 \(-c_{a_i}\),表示上次的“出手”是卖掉,以抵消 本次 \(a_i\) 的花费。
  • \(v_i'\) 向 \(T\) 连边,容量为 \(1\),费用为 \(0\)。

  费用流求出的最小费用就是答案。

\(\mathcal{Code}\)

#include <queue>
#include <cstdio> typedef std::pair<int, int> pii; const int MAXN = 80, MAXND = MAXN * 2 + 2, MAXED = 5 * MAXN, INF = 0x3f3f3f3f;
int n, K, S, T, ecnt = 1, a[MAXN + 5], c[MAXN + 5], las[MAXN + 5];
int head[MAXND + 5], curh[MAXND + 5], d[MAXND + 5];
bool vis[MAXND + 5]; struct Edge { int to, flow, cost, nxt; } graph[MAXED * 2 + 5]; inline void link ( const int s, const int t, const int f, const int c ) {
graph[++ ecnt] = { t, f, c, head[s] };
head[s] = ecnt;
} inline void addEdge ( const int s, const int t, const int f, const int c ) {
link ( s, t, f, c ), link ( t, s, 0, -c );
} inline pii DFS ( const int u, int iflow ) {
vis[u] = true;
if ( u == T ) return { iflow, 0 };
int oflow = 0, ocost = 0;
for ( int& i = curh[u], v; i; i = graph[i].nxt ) {
if ( ! vis[v = graph[i].to] && d[v] == d[u] + graph[i].cost && graph[i].flow ) {
pii of ( DFS ( v, std::min ( iflow, graph[i].flow ) ) );
oflow += of.first, ocost += of.first * graph[i].cost + of.second;
graph[i].flow -= of.first, graph[i ^ 1].flow += of.first;
if ( ! ( iflow -= of.first ) ) break;
}
}
if ( ! oflow ) d[u] = INF;
return { oflow, ocost };
} inline bool SPFA () {
static std::queue<int> que;
static bool inq[MAXND + 5];
for ( int i = 0; i <= T; ++ i ) d[i] = INF, inq[i] = false;
que.push ( S ), d[S] = 0, inq[S] = true;
for ( int u; ! que.empty (); ) {
inq[u = que.front ()] = false, que.pop ();
for ( int i = head[u], v; i; i = graph[i].nxt ) {
if ( d[v = graph[i].to] > d[u] + graph[i].cost && graph[i].flow ) {
d[v] = d[u] + graph[i].cost;
if ( ! inq[v] ) que.push ( v ), inq[v] = true;
}
}
}
return d[T] ^ INF;
} inline int Dinic () {
int ret = 0;
for ( ; SPFA (); ret += DFS ( S, INF ).second ) {
for ( int i = 0; i <= T; ++ i ) {
vis[i] = false;
curh[i] = head[i];
}
}
return ret;
} int main () {
scanf ( "%d %d", &n, &K );
for ( int i = 1; i <= n; ++ i ) scanf ( "%d", &a[i] );
for ( int i = 1; i <= n; ++ i ) scanf ( "%d", &c[i] );
S = 0, T = n << 1 | 1;
for ( int i = 1; i <= n; ++ i ) {
addEdge ( S, i, 1, c[a[i]] ), addEdge ( i, i + n, 1, 0 );
if ( las[a[i]] ) addEdge ( i - 1, las[a[i]] + n, 1, -c[a[i]] );
if ( i < n ) addEdge ( i, i + 1, K - 1, 0 );
addEdge ( i + n, T, 1, 0 ), las[a[i]] = i;
}
printf ( "%d\n", Dinic () );
return 0;
}

Solution -「CF 802C」Heidi and Library (hard)的更多相关文章

  1. Solution -「CF 1342E」Placing Rooks

    \(\mathcal{Description}\)   Link.   在一个 \(n\times n\) 的国际象棋棋盘上摆 \(n\) 个车,求满足: 所有格子都可以被攻击到. 恰好存在 \(k\ ...

  2. Solution -「CF 1622F」Quadratic Set

    \(\mathscr{Description}\)   Link.   求 \(S\subseteq\{1,2,\dots,n\}\),使得 \(\prod_{i\in S}i\) 是完全平方数,并最 ...

  3. Solution -「CF 923F」Public Service

    \(\mathscr{Description}\)   Link.   给定两棵含 \(n\) 个结点的树 \(T_1=(V_1,E_1),T_2=(V_2,E_2)\),求一个双射 \(\varph ...

  4. Solution -「CF 923E」Perpetual Subtraction

    \(\mathcal{Description}\)   Link.   有一个整数 \(x\in[0,n]\),初始时以 \(p_i\) 的概率取值 \(i\).进行 \(m\) 轮变换,每次均匀随机 ...

  5. Solution -「CF 1586F」Defender of Childhood Dreams

    \(\mathcal{Description}\)   Link.   定义有向图 \(G=(V,E)\),\(|V|=n\),\(\lang u,v\rang \in E \Leftrightarr ...

  6. Solution -「CF 1237E」Balanced Binary Search Trees

    \(\mathcal{Description}\)   Link.   定义棵点权为 \(1\sim n\) 的二叉搜索树 \(T\) 是 好树,当且仅当: 除去最深的所有叶子后,\(T\) 是满的: ...

  7. Solution -「CF 623E」Transforming Sequence

    题目 题意简述   link.   有一个 \(n\) 个元素的集合,你需要进行 \(m\) 次操作.每次操作选择集合的一个非空子集,要求该集合不是已选集合的并的子集.求操作的方案数,对 \(10^9 ...

  8. Solution -「CF 1023F」Mobile Phone Network

    \(\mathcal{Description}\)   Link.   有一个 \(n\) 个结点的图,并给定 \(m_1\) 条无向带权黑边,\(m_2\) 条无向无权白边.你需要为每条白边指定边权 ...

  9. Solution -「CF 599E」Sandy and Nuts

    \(\mathcal{Description}\)   Link.   指定一棵大小为 \(n\),以 \(1\) 为根的有根树的 \(m\) 对邻接关系与 \(q\) 组 \(\text{LCA}\ ...

随机推荐

  1. BIO、NIO、AIO --- 个人理解

    1.前言 什么是 BIO.NIO.AIO  ,不难看出,都是共同的字符IO , IO的意思是input output  ,即输入输出 , 那么 B . N .A 分别指不同的io模型 ,而io又分为 ...

  2. iframe页面二次登录问题

    原文链接:iframe页面二次登录问题 生产问题 问题背景 由于历史原因,公司内部系统有一些页面是基于iframe嵌入的其他系统的页面,之前一直运行正常,最近不知什么原因接连出现访问所有iframe页 ...

  3. 如何向内核提交补丁?——FirstKernelPatch

    参考 https://kernelnewbies.org/FirstKernelPatch

  4. BitMap算法知识笔记以及在大数据方向的使用

    概述 所谓的BitMap算法就是位图算法,简单说就是用一个bit位来标记某个元素所对应的value,而key即是该元素,由于BitMap使用了bit位来存储数据,因此可以大大节省存储空间,这是很常用的 ...

  5. 《剑指offer》面试题26. 树的子结构

    问题描述 输入两棵二叉树A和B,判断B是不是A的子结构.(约定空树不是任意一个树的子结构) B是A的子结构, 即 A中有出现和B相同的结构和节点值. 例如: 给定的树 A:      3     / ...

  6. SQLServer触发器调用JavaWeb接口

    这几天接到一个需求需要吧不同系统的数据库进行同步,需要我做一个中间平台进行连接,瞬间就想到了触发器调用接口然后通过API进行传递再写入另一个数据库. sqlServer触发器调用JavaWeb接口 1 ...

  7. Javascript——ES6( ECMAScript 6.0)语法

    ES6( ECMAScript 6.0)语法 一.let/const与var的区别 var 会进行预解析,let/const不会 var可以声明两个重名的变量,let/const不能 var没有块级作 ...

  8. 【小记录】android下opencv的cv::dft()函数,CPU版本与opencl版本的性能相差16倍

            cv::dft  相差15.9倍         cpu版本  单次调用  0.029448 毫秒         opencl版本  单次调用  0.468688  毫秒   差别仅 ...

  9. HashMap和TreeMap的内部结构

    一.HashMap 1.基于哈希表的 Map 接口的实现.此实现提供所有可选的映射操作,并允许使用 null 值和 null 键.(除了非同步和允许使用 null 之外,HashMap 类与 Hash ...

  10. 【自写信息搜集工具】ThunderSearch开发原理解析

    前段时间结合zoomeye的开发文档做了个简易的信息搜集工具ThunderSearch[项目地址 / 博客地址],这次来讲讲具体的实现原理和开发思路 首先要能看懂开发文档,https://www.zo ...