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 ...
随机推荐
- 项目搭建系列之一:使用Maven搭建SpringMVC项目
约定电脑都安装了eclipse,且已配置好Maven以及eclipse插件. 1.Eclipse 2.maven 3.Eclipse 需要安装maven插件.url:maven - http://do ...
- 用php生成静态html页面(通用2种方法)
因为每次用户点击动态链接的时候都会对服务器发送数据查询的要求 对于一个访问量可能达百万千万级别的网站来说 这无疑是服务器一个大大的负担 所以把动态数据转换成静态html页面就成了节省人力物力的首选 因 ...
- node命令
- 免费股票数据API接口
免费股票数据API接口提供沪深.香港.美国股市信息. 1.沪深股市 2.香港股市 3.美国股市 4.香港股市列表 5.美国股市列表 6.深圳股市列表 7.沪股列表 API文档:https://www. ...
- css中单位em和rem
一.介绍 1.em w3cschool中给出css中尺寸单位如下: 单位 描述 % 百分比 in 英寸 cm 厘米 mm 毫米 em 1em 等于当前的字体尺寸. 2em 等于当前字体尺寸的两倍. 例 ...
- USB hacker Collection
This blog contains some ideas and tricks about USB hacking.
- [每日菜单]lunch menu for Wednesday, February 24 2016
HIGHLIGHTS ** Boiled Broccoli原味西兰花 Steamed Turbot with Chili剁椒蒸多宝鱼 * Rye Bread黑麦面包 CHINESE CU ...
- Unity 5.3 安装完没有Android(安卓)或IOS Module(模块)?
Unity5.3 模块独立 Unity5.3把各个模块分开来了,主程序安装包更轻巧,在官网下载的话,能下载到 Unity安装程序,Unity编辑器等一些资源package,其它的模块可以通过Unity ...
- jmeter 与 java http
jmeter 如果对java代码进行测试 1.eclips中创建一个项目,且写一个待测试的简单java代码 2.将jmeter路径下 x:\xx\lxx\Dowxxxxxx\apache-jmeter ...
- poj1416 Shredding Company
Shredding Company Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 5379 Accepted: 3023 ...