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 ...
随机推荐
- 企业邮箱在Android(安卓)系统手机上POP3/IMAP协议的设置方法
此处以三星(系统版本4.4.2)为例,介绍下使用安卓系统自带的客户端如何设置pop/imap协议方式方法 以下我们将使用test@zhuyuming.so 为测试案例,请您操作时更换成您自己的邮箱账号 ...
- nginx日志切割脚本
#!/bin/bash ip=`ifconfig eth0 | grep "inet addr" | cut -f 2 -d ":" | cut -f 1 -d ...
- 自罚一杯-PHP基础(一)
hi 懒癌是最可怕的癌了吧,在即将迎来本命年时,给自己定下每天必干的几件成习惯的事情,望不再鸽!! php是我一直在学习的东西,但总以为php简单,难点在需要掌握大量的相关知识,其实并不完全对,php ...
- 今天有事-MySQL
hi 今天有事,一会儿要去耍,能学多少是多少吧 1.MySQL -----子查询与连接(二)----- ----子查询 子查询,是指出现在其他SQL语句内的SELECT子句 注意:子查询指嵌套在查询内 ...
- 高级c++头文件bits/stdc++.h
用这种方法声明头文件只需两行代码 #include<bits/stdc++.h> using namespace std; 这个头文件包含以下等等C++中包含的所有头文件: #includ ...
- python读取excel并制表输出
源码如下: #!/usr/bin/python #coding=UTF-8 import xlrd import sys from texttable import Texttable def she ...
- 【Unity3d】Ray射线初探-射线的原理及用法
http://www.xiaobao1993.com/231.html 射线是一个无穷的线,开始于origin并沿着direction方向. 当射线碰到物体后.它就会停止发射. 在屏幕中拉一个CUBE ...
- 搞netty
开始搞netty了 查了下资料 在使用NIO 的时候,最好不要配置 SO_LINGER,假设设置了该參数,在 close的时候如缓冲区有数据待写出,会抛出 IOException. // 在netty ...
- MyBatis配置文件解析
MyBatis配置文件解析(概要) 1.configuration:根元素 1.1 properties:定义配置外在化 1.2 settings:一些全局性的配置 1.3 typeAliases:为 ...
- SQLite 增、删、改、查
1.iOS中实现SQLite的增.删.改.查 http://www.jianshu.com/p/0b9b78e704a4. 2.用数据库实现收藏功能 http://www.jianshu.com/p ...