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. x01.os.7: 傻子一样的等

    傻子一样的等 昨日出差,办完事后,下午 2:30,准备进长途汽车站买票回家,被一人拦住,说可以带我进去,车马上就要开了,买票来不及.我以为是汽车司机,就跟了进去.进去后,他打了个电话,说 3:30 车 ...

  2. Hadoop Kernel tunning

    /etc/security/limits.conf @ochadoop soft nofile 102642 @ochadoop hard nofile 102642 @ochadoop soft n ...

  3. 报表引擎API开发入门—简单程序数据集

    小编最近接的项目是有关报表开发的,很想把这部分知识分享出来.希望大家能够支持我!不多说,马上进入我们今天的话题. API基本知识 小编最近项目所做的是关于一个报表软件—FineReport报表开发的一 ...

  4. AWS国际版的Route 53和CloudFront

    注册AWS国际版账号后,却发现Route 53和CloudFront功能是无法使用的.于是提交了一个Service Request,得到的答复是这两个功能需要验证后才能激活. 在控制台中点击进入Rou ...

  5. What makes an excellent front-end developer?(for my English speech)

    What makes an excellent front-end developer? Let me please start this talking by saying that what is ...

  6. Golang tips

    1.go test 测试单个函数 go test -v -test.run Test* 2.

  7. Object-C中self和super的区别

    self与super的区别 原文CSDN evilotus 有所整理 **** 在ObjC中的类实现中经常看到这两个关键字"self"和"super",以以前o ...

  8. 嵌入式Linux驱动学习之路(二)u-boot体验

    u-boot工程简介 现在的u-boot支持PowerPC.ARM.X86.MIPS体系结构的上百种开发板,已经称为功能最多.灵活性最强,并且开发最积极的开源Bootloader.目前由DENX的Wo ...

  9. Github 安全类Repo收集整理

    作者:天谕链接:https://zhuanlan.zhihu.com/p/21380662来源:知乎著作权归作者所有.商业转载请联系作者获得授权,非商业转载请注明出处.刚好这两天对之前github上关 ...

  10. uva131 The Psychic Poker Player

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