题目链接

http://acm.hdu.edu.cn/showproblem.php?pid=5550

题意

一撞大楼有N层楼,然后每层楼都有一部分人喜欢打羽毛球,一部分人喜欢打乒乓球

但是每层楼只能选择建一个羽毛球馆或者建一个乒乓球馆 那么每个人到它喜欢的球馆的距离就是一个权值

求出怎么规划 使得所有人到它喜欢的球馆的距离之和最小

思路

其实当时在训练的时候 有在想

当时训练的时候的思路是

dp[i][j] i 表示第几层 j 表示状态 0 表示乒乓球馆 1 表示羽毛球馆

然后每次转移的是从前面那一层转移过来的

因为那时候是这样想的,假设前面的规划是最优解,也就是说 dp[i - 1][0] 或者 dp[i - 1][1] 保存的是最优状态

那么我再加一层 是不是可以根据 这个最优状态转移过来,我加的这一层也就只有两种选择,0 或者 1

那么我只要用一个数组来维护前面状态离它最近的乒乓球馆和羽毛球馆的位置, 然后更新一下就好了。

然后jsw告诉我这样是错的,,然后我现在也不知道这样为什么是错的,。。 它是说 应该是从1 -> i - 1 都更新一遍,才是对的。。。 然后他还是没告诉我,为什么我这样是错的。。。。

然后查了题解,果然是要这样更新,

dp[i][j] i 表示第几层,j 表示状态

dp[i][j] = min(dp[k][j ^ 1] + valuesum(k + 1, i, j)) k < i

valuesum 表示 区间[i,. j] 的人到达他们最近的球馆的距离总和

其实更新是这样的

是 k 层 和 i + 1 层的状态是一样的,, 然后 k + 1 -> i 层的状态是一样的

那么 k + 1 -> i 层的和其楼层状态不同的人 前一半就去 k 层 后一半就去 i + 1 层 这样是最优的

如果暴力更新会T 就要用到前缀和 研究了好久。。

一个是人数的前缀和,一个是距离的前缀和

人数的前缀和很简单,就是前i层的总人数

我们先从低维度进行思考

如果我知道一个人在x 层 ,那么我要求他去y层的距离是多少 那么很显然是不是 abs(x - y)

那么很多人 其实也是一样的,,我们只需要记录每层有多少人,然后求到Y的距离 就是 abs(所有人数所在的楼层和 - 人数 * Y)

sum 数组存放的是人数前缀和, value 数组存放的是每层楼有多少人的前缀和

AC代码

#pragma comment(linker, "/STACK:102400000,102400000")

#include <cstdio>
#include <cstring>
#include <ctype.h>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <deque>
#include <vector>
#include <queue>
#include <string>
#include <map>
#include <stack>
#include <set>
#include <list>
#include <numeric>
#include <sstream>
#include <iomanip>
#include <limits> #define CLR(a, b) memset(a, (b), sizeof(a));
#define pb push_back
#define bug puts("***bug***");
#define fi first
#define se second
#define L(on) ((on)<<1)
#define R(on) (L(on) | 1)
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define syn_close ios::sync_with_stdio(false); cin.tie(0);
#define sp system("pause");
//#define gets gets_s using namespace std; typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef long double ld;
typedef pair <int, int> pii;
typedef pair <ll, ll> pll;
typedef pair <string, int> psi;
typedef pair <string, string> pss;
typedef pair <double, int> pdi; const double PI = acos(-1.0);
const double EI = exp(1.0);
const double eps = 1e-8; const int INF = 0x3f3f3f3f;
const ll INFLL = 0x3f3f3f3f3f3f3f3fll;
const int maxn = (int)4e3 + 10;
const int MAXN = (int)1e2 + 10;
const int MOD = (int)1e9 + 7; int n;
pll arr[maxn]; // fi 0 se 1
ll sum[2][maxn], value[2][maxn], dp[2][maxn]; int readint()
{
int num; scanf("%d", &num);
return num;
} void input()
{
n = readint();
CLR(sum, 0); CLR(value, 0); CLR(dp, 0x3f);
for (int i = 1; i <= n; i++)
{
arr[i].fi = (ll)readint(), arr[i].se = (ll)readint();
sum[0][i] = sum[0][i - 1] + arr[i].fi;
sum[1][i] = sum[1][i - 1] + arr[i].se;
value[0][i] = value[0][i - 1] + arr[i].fi * i;
value[1][i] = value[1][i - 1] + arr[i].se * i;
}
} ll leftcal(int l, int r, int vis)
{
ll x = value[vis][r] - value[vis][l - 1];
ll y = sum[vis][r] - sum[vis][l - 1];
return abs(x - y * (l - 1));
} ll rightcal(int l, int r, int vis)
{
ll x = value[vis][r] - value[vis][l - 1];
ll y = sum[vis][r] - sum[vis][l - 1];
return abs(x - y * (r + 1));
} ll valuesum(int l, int r, int vis)
{
if (l == 1) return rightcal(l, r, vis);
if (r == n) return leftcal(l, r, vis);
int mid = (l + r) >> 1;
return leftcal(l, mid, vis) + rightcal(mid + 1, r, vis);
} void solve()
{
for (int i = 1; i < n; i++)
for (int j = 0; j < 2; j++)
dp[j][i] = valuesum(1, i, j);
for (int i = 1; i <= n; i++)
for (int j = 1; j < i; j++)
for (int k = 0; k < 2; k++)
dp[k][i] = min(dp[k][i], dp[k ^ 1][j] + valuesum(j + 1, i, k));
printf("%lld\n", min(dp[0][n], dp[1][n]));
} int main()
{
int t = readint();
for (int tt = 1; tt <= t; tt++)
{
printf("Case #%d: ", tt);
input(); solve();
}
}

HDU - 5550 Game Rooms 【DP+前缀和】的更多相关文章

  1. HDU 5550 - Game Rooms(DP + 前缀和预处理)

    链接: http://acm.hdu.edu.cn/showproblem.php?pid=5550 题意: 一个大楼有n(2≤n≤4000)层,每层可以建一个乒乓球房或者一个游泳房,且每种房间在大楼 ...

  2. HDU 2829 区间DP & 前缀和优化 & 四边形不等式优化

    HDU 2829 区间DP & 前缀和优化 & 四边形不等式优化 n个节点n-1条线性边,炸掉M条边也就是分为m+1个区间 问你各个区间的总策略值最少的炸法 就题目本身而言,中规中矩的 ...

  3. HDU 1011 树形背包(DP) Starship Troopers

    题目链接:  HDU 1011 树形背包(DP) Starship Troopers 题意:  地图中有一些房间, 每个房间有一定的bugs和得到brains的可能性值, 一个人带领m支军队从入口(房 ...

  4. HDU 1559 最大子矩阵 (DP)

    题目地址:pid=1559">HDU 1559 构造二维前缀和矩阵.即矩阵上的点a[i][j]表示左上方的点为(0,0),右下方的点为(i,j)的矩阵的和.然后枚举每一个矩阵的左上方的 ...

  5. hdu 2296 aC自动机+dp(得到价值最大的字符串)

    Ring Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submis ...

  6. HDU 4778 状压DP

    一看就是状压,由于是类似博弈的游戏.游戏里的两人都是绝对聪明,那么先手的选择是能够确定最终局面的. 实际上是枚举最终局面情况,0代表是被Bob拿走的,1为Alice拿走的,当时Alice拿走且满足变换 ...

  7. HDOJ(HDU).3466 Dividing coins ( DP 01背包 无后效性的理解)

    HDOJ(HDU).3466 Dividing coins ( DP 01背包 无后效性的理解) 题意分析 要先排序,在做01背包,否则不满足无后效性,为什么呢? 等我理解了再补上. 代码总览 #in ...

  8. HDOJ(HDU).2546 饭卡(DP 01背包)

    HDOJ(HDU).2546 饭卡(DP 01背包) 题意分析 首先要对钱数小于5的时候特别处理,直接输出0.若钱数大于5,所有菜按价格排序,背包容量为钱数-5,对除去价格最贵的所有菜做01背包.因为 ...

  9. HDOJ(HDU).2602 Bone Collector (DP 01背包)

    HDOJ(HDU).2602 Bone Collector (DP 01背包) 题意分析 01背包的裸题 #include <iostream> #include <cstdio&g ...

随机推荐

  1. ES6 动态计算属性名

    在ES5之前,如果属性名是个变量或者需要动态计算,则只能通过 对象.[变量名] 的方式去访问. <script type="text/javascript"> var ...

  2. 3、Linux内核模块学习

    一.内核模块的学习   内核的整体框架是非常的大,包含的组件也是非常多,如何将需要的组件包含在内核中呢?选择一,就是将所有的组件全部编译进内核,虽然需要的组件都可以使用,但是内核过分庞大,势必带来效率 ...

  3. WCF实现上传图片功能

    初次学习实现WCF winform程序的通信,主要功能是实现图片的传输. 下面是实现步骤: 第一步: 首先建立一个类库项目TransferPicLib,导入wcf需要的引用System.Service ...

  4. NUTCH2.3 hadoop2.7.1 hbase1.0.1.1 solr5.2.1部署(三)

     Precondition: hadoop 2.7.1 hbase 0.98.13 solr 5.2.1 / Apache Solr 4.8.1 http://archive.apache.org ...

  5. Sphinx之配置文件

    # # Sphinx configuration file sample # # WARNING! While this sample file mentions all available opti ...

  6. Angular $httpProvider

    timeout超时响应 .factory('timestampMarker', ["$rootScope", function () { var timestampMarker = ...

  7. bzoj4069【APIO2015】巴厘岛的雕塑

    4069: [Apio2015]巴厘岛的雕塑 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 192  Solved: 89 [Submit][Stat ...

  8. 【JMeter4.0学习(十一)】之JMeter对(Mysql、Oracle)数据库性能测试脚本开发

    一.MySQL数据库链接: 注:下面所产生的问题一律参考详见:<[JMeter4.0]之遇到的问题总结(持续更新)>(包括Mysql.Orcale) 准备:引包,包路径一定要放对位置,参考 ...

  9. php 生成下载连接

    public function showdownload(){ $file_url=$_GET['url']; $new_name='激活码'; if(!isset($file_url)||trim( ...

  10. 阿里云服务器---centos编译安装ffmpeg

    环境 系统环境:CentOS release 6.7 (Final) 需求 编译安装ffmpeg 获取依赖 安装依赖包 yum install -y autoconf automake cmake f ...