1. 链接

    Desiderium

  2. 题意

      给定n条线段,从中选取若干条,共有2n种选法(因为每一条线段有两种方法:选或者不选).

      每一种选法都对应一个长度,也就是所选线段的并集长度.

      求这2n种选法长度之和.

  3. 解法:

    把这n条线段进行离散化,使得任意一条线段都可以由若干条元线段组成.

    什么叫元线段呢?把全部的x坐标进行排序,去重,就得到了很多元线段.

    统计每条元线段出现的次数,也就是说它被这n条线段里面的多少个线段覆盖.

    如合统计呢?假设元线段被m条线段覆盖,那么有n-m条线段不覆盖它.

    所以,该元线段被2n种选法中的2n-2n-m种选法覆盖.

    ans=累加(元线段的长度*元线段使用的次数).

  4. 写法:

    主要是如何统计元线段出现的次数.

    • 建立线段树,把没一条线段进行插入O(nlgn)
    • 树状数组O(nlgn)
    • 普通数组O(n),这种方法最好,不仅简单,而且快.
  5. 代码(树状数组版):

    #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;
    }
    
  6. 代码(普通数组版):

     #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的更多相关文章

  1. HDU 5481 Desiderium 动态规划

    Desiderium Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=548 ...

随机推荐

  1. Jmeter默认报告优化

    一.本文目的: 之前写了两篇文章搭建持续集成接口测试平台(Jenkins+Ant+Jmeter)和ANT批量执行Jmeter脚本,功能实现上都没有什么问题,但是最后生成的报告有一点小问题,虽然不影响使 ...

  2. x01.Game.Main: 从零开始

    一切从零开始,一切皆有可能. 浅墨,90后,<逐梦之旅>深入浅出,堪比大师. 1.安装 DXSDK_June10.exe 或更新版本. 2.运行 vs2012,新建 VC Win32 空项 ...

  3. 怎么使PUTTY一直连接

    如何才能保证PUTTY一直连接,即使我们好长时间不去敲命令,也让它一直连接着? PuTTY的设置方法是:在Connection里面有个Seconds between keepaliaves.这里就是每 ...

  4. OpenCV 之 边缘检测

    上一篇 <OpenCV 之 图像平滑> 中,提到的图像平滑,从信号处理的角度来看,实际上是一种“低通滤波器”. 本篇中,数字图像的边缘,因为通常都是像素值变化剧烈的区域 (“高频”),故可 ...

  5. Visual Studio 代码折叠快捷键(摘要)

    代码编辑器的展开和折叠代码确实很方便和实用.以下是展开代码和折叠代码所用到的快捷键,很常用: Ctrl + M + O: 折叠所有方法 Ctrl + M + M: 折叠或者展开当前方法 Ctrl +  ...

  6. poj2486Apple Tree[树形背包!!!]

    Apple Tree Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 9989   Accepted: 3324 Descri ...

  7. Linux下mysql新建账号及权限设置

    http://www.cnblogs.com/eczhou/archive/2012/07/12/2588187.html 1.权限赋予 说明:mysql部署在服务器A上,内网上主机B通过客户端工具连 ...

  8. 烈焰SWF解密

    SWF 解密 是用UE编辑器 改 SWF开头 的AA AA AA ,改成43 57 53 就解密了

  9. dp入门问题

    昨天晚上的rank彻底废了..一个星期没敲代码完全没手感.作为总结,贴一道昨天浪费了我两小时的dp.http://acm.dirring.com/problem.php?cid=1003&pi ...

  10. uva131 The Psychic Poker Player

    The Psychic Poker Player Time Limit: 3000MS     64bit IO Format: %lld & %llu Description In 5-ca ...