题目描述
LYK在学习dp,有一天它看到了一道关于dp的题目。
这个题目是这个样子的:一开始有n个数,一段区间的价值为这段区间相同的数的对数。我们想把这n个数切成恰好k段区间。之后这n个数的价值为这k段区间的价值和。我们想让最终这n个数的价值和尽可能少。
例如6个数1,1,2,2,3,3要切成3段,一个好方法是切成[1],[1,2],[2,3,3],这样只有第三个区间有1的价值。因此这6个数的价值为1。
LYK并不会做,丢给了你。

输入格式(dp.in)
第一行两个数n,k。
接下来一行n个数ai表示这n个数。

输出格式(dp.out)
一个数表示答案。

输入样例
10 2
1 2 1 2 1 2 1 2 1 2

输出样例
8

数据范围
对于30%的数据n<=10。
对于60%的数据n<=1000。
对于100%的数据1<=n<=100000,1<=k<=min(n,20),1<=ai<=n。
其中有30%的数据满足ai完全相同均匀分布在所有数据中。

分析:很容易想到dp的做法,设f[i][j]表示前i个数中分成了k个区间的最小价值和.

f[i][j] = min{f[k][j-1] + sum[k+1][i]}.这样dp只有60分,需要优化.注意到当j不变时,i增长,k是一定不减的.因为如果f[i+1][j]能从k前面转移过来,那么f[i][j]一定也能从k前面转移过来,这与f[i][j]从k转移过来是自相矛盾的,所以k一定不减.那么可以用1D1D优化.

类似分治,如果f[mid][j]从f[pos][j-1]转移而来,那么f[l~mid-1][j]一定从f[l'~pos][j-1]转移而来,f[mid + 1~r][j]一定从f[pos~r'][j-1]转移而来,递归下去,就能更新完所有答案.因为涉及到多次的区间求相同数对的个数,可以像莫队一样用两个指针来维护.

    复杂度一般就是把朴素dp的一个n变成log,这道题的复杂度就是O(n*logn*k).

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm> using namespace std;
typedef long long ll;
const ll inf = 1LL << ;
ll f[], n, k, a[], c[], g[], L = , R = , tot, pos, mx; void move(int l, int r)
{
while (L < l)
{
c[a[L]]--;
tot -= c[a[L]];
L++;
}
while (L > l)
{
L--;
tot += c[a[L]];
c[a[L]]++;
}
while (R < r)
{
R++;
tot += c[a[R]];
c[a[R]]++;
}
while (R > r)
{
c[a[R]]--;
tot -= c[a[R]];
R--;
}
} void solve(ll l, ll r, ll x, ll y)
{
if (x > y)
return;
ll mid = (x + y) >> ;
mx = inf;
for (ll i = l; i <= r; i++)
if (i < mid) //mid只能从比mid小的地方转移过来
{
move(i + , mid);
if (f[i] + tot < mx)
{
mx = f[i] + tot;
pos = i;
}
}
g[mid] = mx;
solve(l, pos, x, mid - );
solve(pos, r, mid + , y);
} int main()
{
scanf("%lld%lld", &n, &k);
for (int i = ; i <= n; i++)
scanf("%lld", &a[i]);
for (int i = ; i <= n; i++)
f[i] = inf;
while (k--)
{
memset(c, , sizeof(c));
L = ;
R = ;
tot = ;
solve(, n - , , n);
memcpy(f, g, sizeof(g));
memset(g, , sizeof(g));
}
printf("%lld\n", f[n]); return ;
}

noip模拟赛 动态规划的更多相关文章

  1. contesthunter暑假NOIP模拟赛第一场题解

    contesthunter暑假NOIP模拟赛#1题解: 第一题:杯具大派送 水题.枚举A,B的公约数即可. #include <algorithm> #include <cmath& ...

  2. 2016-06-19 NOIP模拟赛

          2016-06-19 NOIP模拟赛 by coolyangzc 共3道题目,时间3小时 题目名 高级打字机 不等数列 经营与开发 源文件 type.cpp/c/pas num.cpp/c ...

  3. NOIP模拟赛20161022

    NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...

  4. NOIP模拟赛 by hzwer

    2015年10月04日NOIP模拟赛 by hzwer    (这是小奇=> 小奇挖矿2(mining) [题目背景] 小奇飞船的钻头开启了无限耐久+精准采集模式!这次它要将原矿运到泛光之源的矿 ...

  5. 大家AK杯 灰天飞雁NOIP模拟赛题解/数据/标程

    数据 http://files.cnblogs.com/htfy/data.zip 简要题解 桌球碰撞 纯模拟,注意一开始就在袋口和v=0的情况.v和坐标可以是小数.为保险起见最好用extended/ ...

  6. 队爷的讲学计划 CH Round #59 - OrzCC杯NOIP模拟赛day1

    题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的讲学计划 题解:刚开始理解题意理解了好半天,然后发 ...

  7. 队爷的Au Plan CH Round #59 - OrzCC杯NOIP模拟赛day1

    题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的Au%20Plan 题解:看了题之后觉得肯定是DP ...

  8. 队爷的新书 CH Round #59 - OrzCC杯NOIP模拟赛day1

    题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的新书 题解:看到这题就想到了 poetize 的封 ...

  9. CH Round #58 - OrzCC杯noip模拟赛day2

    A:颜色问题 题目:http://ch.ezoj.tk/contest/CH%20Round%20%2358%20-%20OrzCC杯noip模拟赛day2/颜色问题 题解:算一下每个仆人到它的目的地 ...

随机推荐

  1. php生成唯一订单号的方法

    第一种 $danhao = date('Ymd') . str_pad(mt_rand(1, 99999), 5, '0', STR_PAD_LEFT); 第二种 $danhao = date('Ym ...

  2. Classic BADI总结

    这里对sap Classic Badi 做一下总结,虽然已经是过时的技术了. Classic BADI的创建 Classic BADI的实施 Classic BADI的调用及运行原理 New BADI ...

  3. LOJ#503. 「LibreOJ β Round」ZQC 的课堂(容斥+FHQTreap)

    题面 传送门 题解 首先\(x\)和\(y\)两维互相独立,可以分开考虑,我们以\(x\)为例 我们把\(x\)做个前缀和,那么就是问有多少\(i\)满足\(s_is_{i-1}<0\),其中\ ...

  4. 1.2打印ASCII码

    描述 输入一个除空格以外的可见字符(保证在函数scanf中可使用格式说明符%c读入),输出其ASCII码. 输入一个除空格以外的可见字符.输出一个十进制整数,即该字符的ASCII码.样例输入 A 样例 ...

  5. 385 Mini Parser 迷你解析器

    Given a nested list of integers represented as a string, implement a parser to deserialize it.Each e ...

  6. HTML5 WEB Storage - localStorage存储位置在哪

    localStorage作为客户端浏览器持久化存储方案 这个是浏览器隔离的,每个浏览器都会把localStorage存储在自己的UserData中,如chrome一般就是 C:\Users\你的计算机 ...

  7. iPhone4到iPhone6的设计、制造工艺历程浅析

    这里就阐述一下我对这几个手机在设计和制造工艺上的一些看法. 这 5 个型号概括起来就三个外观: iPhone4 的三明治夹心设计. iPhone5 的三段式铝合金一体成型. iPhone6 的全金属一 ...

  8. BPI-M1P(全志A20)刷Android启动卡之后启动的过程

    http://blog.csdn.net/wb4916/article/details/78031511BPI-M1P(全志A20)刷Android启动卡之后启动的过程 BPI-M1P(全志A20)刷 ...

  9. Learning Face Age Progression: A Pyramid Architecture of GANs

    前言 作为IP模式识别的CNN初始模型是作为单纯判别式-模式识别存在的,并以此为基本模型扩展到各个方向.基本功能为图像判别模型,此后基于Loc+CNN的检测模型-分离式.end2end.以及MaskC ...

  10. 重启rsyncd

    systemctl  restart  rsyncd.service