hdu5481 Desiderium
链接
题意
给定n条线段,从中选取若干条,共有2n种选法(因为每一条线段有两种方法:选或者不选).
每一种选法都对应一个长度,也就是所选线段的并集长度.
求这2n种选法长度之和.
解法:
把这n条线段进行离散化,使得任意一条线段都可以由若干条元线段组成.
什么叫元线段呢?把全部的x坐标进行排序,去重,就得到了很多元线段.
统计每条元线段出现的次数,也就是说它被这n条线段里面的多少个线段覆盖.
如合统计呢?假设元线段被m条线段覆盖,那么有n-m条线段不覆盖它.
所以,该元线段被2n种选法中的2n-2n-m种选法覆盖.
ans=累加(元线段的长度*元线段使用的次数).
写法:
主要是如何统计元线段出现的次数.
- 建立线段树,把没一条线段进行插入O(nlgn)
- 树状数组O(nlgn)
- 普通数组O(n),这种方法最好,不仅简单,而且快.
代码(树状数组版):
#include<iostream> #include<string.h> #include<stdio.h> #include<math.h> #include<stdlib.h> #include<algorithm> using namespace std; #define re(i,n) for(int i=0;i<n;i++) typedef long long ll; const int maxn = 1e5 + 7; const int mod = 1e9 + 7; int n; struct Re{ int l, r; }a[maxn]; int x[maxn << 1], xi; int tr[maxn << 1]; #define lson(x) x<<1,f,mid #define rson(x) x<<1|1,mid+1,r int lowbit(int x){ return x&-x; } void ins(int x, int v){ for (int i = x; i > 0; i -= lowbit(i)){ tr[i] += v; } } int query(int x){ int ans = 0; for (int i = x; i < xi; i += lowbit(i)){ ans += tr[i]; } return ans; } int data[maxn]; int f(int m){ return data[n] - data[n - m]; } void init(){ data[0] = 1; re(i, maxn - 2){ data[i + 1] = ((ll)data[i] << 1) % mod; } } int cnt[maxn << 1]; int main(){ //freopen("in.txt", "r", stdin); init(); int T; cin >> T; while (T--){ scanf("%d", &n); xi = -1; re(i, n){ scanf("%d%d", &a[i].l, &a[i].r); x[++xi] = a[i].l, x[++xi] = a[i].r; } ++xi; sort(x, x + xi); xi = unique(x, x + xi) - x; memset(tr, 0, sizeof(tr)); re(i, n){ int beg = lower_bound(x, x + xi, a[i].l) - x; int over = lower_bound(x, x + xi, a[i].r) - x; ins(beg, -1), ins(over, 1); } cnt[0] = 0; for (int i = 1; i < xi; i++)cnt[i] = query(i); //re(i, xi)printf("(%d,%d) ", i, cnt[i]); //puts(""); ll ans = 0; re(i, xi - 1){ ll len = (ll)x[i + 1] - x[i]; int ge = cnt[i + 1]; //cout << len << " " << ge << endl; ans = (len*f(ge) + ans) % mod; } if (ans < 0)ans += mod; printf("%lld\n", ans); } return 0; }代码(普通数组版):
#include<iostream> #include<algorithm> #include<stdio.h> using namespace std; typedef long long ll; #define re(i,n) for(int i=0;i<n;i++) ; ; struct Node{ int l, r; }a[maxn]; int n; ], xi; ]; int data[maxn]; void init(){ data[] = ; ; i < maxn-; i++){ data[i] = (data[i - ] << ) % mod; } } int main(){ //freopen("in.txt", "r", stdin); int T; cin >> T; init(); while (T--){ scanf("%d", &n); xi = -; re(i, n)scanf("%d%d", &a[i].l, &a[i].r), x[++xi] = a[i].l, x[++xi] = a[i].r; sort(x, x + xi + ); xi = unique(x, x + xi + ) - x; re(i, xi)cnt[i] = ; re(i, n){ int f = lower_bound(x, x + xi, a[i].l) - x, t = lower_bound(x, x + xi, a[i].r) - x; cnt[t]++, cnt[f]--; } ; i >= ; i--){ cnt[i] += cnt[i + ]; } ll ans = ; re(i, xi-){ ] - x[i]; ]; ans = (ans + (ll)len*(data[n] - data[n - m])) % mod; } )ans += mod; printf("%lld\n", ans); } ; }
hdu5481 Desiderium的更多相关文章
- HDU 5481 Desiderium 动态规划
Desiderium Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=548 ...
随机推荐
- Jmeter默认报告优化
一.本文目的: 之前写了两篇文章搭建持续集成接口测试平台(Jenkins+Ant+Jmeter)和ANT批量执行Jmeter脚本,功能实现上都没有什么问题,但是最后生成的报告有一点小问题,虽然不影响使 ...
- x01.Game.Main: 从零开始
一切从零开始,一切皆有可能. 浅墨,90后,<逐梦之旅>深入浅出,堪比大师. 1.安装 DXSDK_June10.exe 或更新版本. 2.运行 vs2012,新建 VC Win32 空项 ...
- 怎么使PUTTY一直连接
如何才能保证PUTTY一直连接,即使我们好长时间不去敲命令,也让它一直连接着? PuTTY的设置方法是:在Connection里面有个Seconds between keepaliaves.这里就是每 ...
- OpenCV 之 边缘检测
上一篇 <OpenCV 之 图像平滑> 中,提到的图像平滑,从信号处理的角度来看,实际上是一种“低通滤波器”. 本篇中,数字图像的边缘,因为通常都是像素值变化剧烈的区域 (“高频”),故可 ...
- Visual Studio 代码折叠快捷键(摘要)
代码编辑器的展开和折叠代码确实很方便和实用.以下是展开代码和折叠代码所用到的快捷键,很常用: Ctrl + M + O: 折叠所有方法 Ctrl + M + M: 折叠或者展开当前方法 Ctrl + ...
- poj2486Apple Tree[树形背包!!!]
Apple Tree Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9989 Accepted: 3324 Descri ...
- Linux下mysql新建账号及权限设置
http://www.cnblogs.com/eczhou/archive/2012/07/12/2588187.html 1.权限赋予 说明:mysql部署在服务器A上,内网上主机B通过客户端工具连 ...
- 烈焰SWF解密
SWF 解密 是用UE编辑器 改 SWF开头 的AA AA AA ,改成43 57 53 就解密了
- dp入门问题
昨天晚上的rank彻底废了..一个星期没敲代码完全没手感.作为总结,贴一道昨天浪费了我两小时的dp.http://acm.dirring.com/problem.php?cid=1003&pi ...
- uva131 The Psychic Poker Player
The Psychic Poker Player Time Limit: 3000MS 64bit IO Format: %lld & %llu Description In 5-ca ...