经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C » 查看文章
cf868F. Yet Another Minimization Problem(决策单调性 分治dp)
来源:cnblogs  作者:自为风月马前卒  时间:2018/10/10 9:09:25  对本文有异议

题意

题目链接

给定一个长度为\(n\)的序列。你需要将它分为\(m\)段,每一段的代价为这一段内相同的数的对数,最小化代价总和。

\(n<=10^5,m<=20\)

Sol

看完题解之后的感受:

首先列出裸的dp方程,\(f[i][j]\)表示前\(i\)个位置,切了\(j\)次,转移的时候枚举上一次且在了哪儿

\(f[i][j] = max(f[k][j - 1] + w(k, i))\)

\(w(k, i)\)表示\([k, i]\)内相同的数的对数。。

然后sb的我以为拿个单调队列维护一下就完了结果发现转移是\(O(n)\)的??

标算真神仙Orz

因为转移不是\(O(n)\)的,所以我们可以分治的去做

假设当前要求的区间为\([l, r]\),可以从\([L, R]\)转移而来,\(mid = (l + r)/ 2\)的决策点为\(p\)

那么\([l, mid - 1]\)的转移区间一定在\([L, p]\)\([mid+1, r]\)的转移区间一定是\([p, R]\),递归的做即可

画出图来应该是这样的

求解区间:\(|\gets\)预处理\(\to |\) \(l\frac{\qquad\qquad\qquad\downarrow^{mid}\qquad\qquad\qquad}{}r\)

决策区间:\(L\frac{\qquad\qquad\qquad\downarrow^{p}\qquad\qquad\qquad}{}R\)

转移的时候需要记录每个数的出现次数,每次转移时\(O(1)\)

总的时间复杂度为:\(O(nlognk)\)

推荐一篇写的不错的blog

  1. #include<bits/stdc++.h>
  2. #define LL long long
  3. using namespace std;
  4. const int MAXN = 1e5 + 10;
  5. inline int read() {
  6. int x = 0, f = 1; char c = getchar();
  7. while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
  8. while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
  9. return x * f;
  10. }
  11. int N, K, a[MAXN], b[MAXN];
  12. LL f[MAXN], g[MAXN];
  13. void solve(int l, int r, int L, int R, LL w) {
  14. if(l > r) return ;
  15. int mid = l + r >> 1, p = min(mid, R), k = 0;
  16. for(int i = l; i <= mid; i++) w += b[a[i]]++;
  17. for(int i = L; i <= p; i++) w -= --b[a[i]], f[mid] > g[i] + w ? f[mid] = g[i] + w, k = i : 0;
  18. for(int i = l; i <= mid; i++) w -= --b[a[i]];
  19. for(int i = L; i <= p; i++) w += b[a[i]]++;
  20. solve(l, mid - 1, L, k, w);
  21. for(int i = L; i < k; i++) w -= --b[a[i]];
  22. for(int i = l; i <= mid; i++) w += b[a[i]]++;
  23. solve(mid + 1, r, k, R, w);
  24. for(int i = l; i <= mid; i++) --b[a[i]];
  25. for(int i = L; i < k; i++) ++b[a[i]];
  26. }
  27. main() {
  28. N = read(); K = read();
  29. for(int i = 1; i <= N; i++) a[i] = read(), g[i] = g[i - 1] + b[a[i]]++; memset(b, 0, sizeof(b));
  30. while(K--) {memset(f, 0x7f, sizeof(f)); solve(1, N, 1, N, 0); swap(f, g);}
  31. cout << f[N];
  32. }
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号