zoukankan      html  css  js  c++  java
  • 防晒

    题目

    题目

    题解

    又是一个简(e)单(xin)题,思路不难想,但是证明是真的难啊,不过yxc的视频讲解也是真的好啊。

    我的思路是对于(r)从小到大排序,然后对于防晒霜也从小到大排序,然后对于目前的防晒霜看看有没有区间包括它的,包括就选。(优先选(r)小的)

    当然,题解大量的思路都是(l)递减,然后看看这个区间有没有防晒霜,有优先选(SPF)值大的。

    为什么是对的呢?(下面都按照大多数人的思想走)

    一个比较形象的思路就是:由于(l)递减,(SPF)越大的防晒霜在后面选择会越来越少,所以优先选大的,而且对于(x,y(SPF[x]<SPF[y]))防晒霜而言,后面如果覆盖了(y)号防晒霜那也肯定覆盖了(x)号防晒霜,所以选(SPF)最大的防晒霜是最优的。

    在这里插入图片描述

    那么严谨的证明是什么呢?

    我们把问题化成二分图匹配的问题,然后证明我们这样子做没有增广路径就行了(不懂的看视频或者上网看),而对于(cover[i]),我们就认为有(cover[i])个这样的点就行了。

    (c_i)表示第(i)头牛,(s_i)表示第(i)个防晒霜,那么如果存在增广路径的话,我们找到其中的一条最短的增广路径(注:(a_i)仅仅表示第(i)个点的下标,并无实际意义):

    (s_{a_{0}}→c_{a_{1}}→s_{a_{2}}→...→c_{a_k})
    那么可以说明(c_{a_{0}})(s_{a_k})是没有匹配过的点。

    考虑一下(SPF[a_0])(SPF[a_2])的大小关系,(c_{a_1})包括这两个防晒霜,根据算法大的优先,所以(SPF[a_0]≤SPF[a_2])

    接下来我们再来研究一下(a_1)(a_3)的大小关系。

    有没有可能(a_1<a_3)

    也就是说我们的(s_{a_2})跑去找了后面的牛。

    在这里插入图片描述
    由于的大优先,这样的话下面那个蓝色的防晒霜也就是(SPF[a_4])肯定也是大于(SPF[a_0]),那这样子的话为什么不直接让(s_{a_0})去占(s_{a_4})的位置把(s_{a_4})弹出啊,这样还能少两个点:(s_{a_2})(c_{a_1}),这样就违反了我们最小增广路径的设定了。

    所以(a_1>a_3)

    而后面的也可以像这样子无限推下去,知道(a_k=1),等会!(1)号牛没有匹配?可是(1)号牛优先级最高啊,矛盾,所以不成立。

    还有人问有没有可能是(c_{a_{0}}→s_{a_{1}}→c_{a_{2}}→...→s_{a_k})这样子的最短增广路径呢,但是由于前面的优先级比后面的优先级高,所以(a_0)不可能小于(a_2),所以(a_0>a_2)。抢完之后呢,(c_{a_2})就十分的迷茫了,它要找那个防晒霜呢,假设找到了(x)号防晒霜,(SPF[x]<SPF[a_2]),那么肯定是被是被([a_2+1,a_0-1])的牛给匹配了(不然(a_0)就会匹配到(x)),那么这样,我们我们可以直接让(c_{a_0})直接去抢(c_{a_4})的啊,这样还能少(a_1)(a_2)呢,又违反了最短,而如果(SPF[x]≥SPF[a_2]),确实可以成立,但是又跟上面第一个证明一样,这两条结论一成立,无限推下去知道找到(1)号牛,矛盾。

    所以是对的。

    巧妙的把贪心化成二分匹配进行证明,妙啊。

    当然我的思路其实是和它差不多的,所以也可以认为是对的。

    代码

    时间复杂度:(O(n^2+m))

    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #define  N  3100
    using  namespace  std;
    struct  node
    {
    	int  x,y;
    }a[N],b[N];int  n,m;
    bool  v[N];
    inline  bool  cmp1(node  x,node  y){return  x.y==y.y?(x.x<y.x):(x.y<y.y);}
    inline  bool  cmp2(node  x,node  y){return  x.x<y.x;}
    int  main()
    {
    	scanf("%d%d",&n,&m);
    	for(int  i=1;i<=n;i++)scanf("%d%d",&a[i].x,&a[i].y);
    	for(int  i=1;i<=m;i++)scanf("%d%d",&b[i].x,&b[i].y);
    	sort(a+1,a+n+1,cmp1);
    	sort(b+1,b+m+1,cmp2);
    	int  ans=0;
    	for(int  i=1;i<=m;i++)
    	{
    		for(int  k=1;k<=b[i].y;k++)
    		{
    			bool  bk=0;
    			for(int  j=1;j<=n;j++)
    			{
    				if(!v[j]  &&  a[j].x<=b[i].x  &&  a[j].y>=b[i].x)
    				{
    					v[j]=1;
    					ans++;
    					bk=1;
    					break;
    				}
    			}
    			if(!bk)break;
    		}
    	}
    	printf("%d
    ",ans);
    	return  0;
    }
    

    最后

    这里说一下,其实找到区间内的最小的数字或者找数字可以匹配的(r)最小的区间编号是可以直接用权值线段树的,懒得打了,是可以优化到(O(nlogn))

    当然也可以用https://www.acwing.com/solution/content/785/所用的二分法来找,也是(O(nlogn))的。

  • 相关阅读:
    Nginx禁止IP,只允许指定域名访问
    预防vsphere勒索病毒,适用于 ESXi 6.x 中的 OpenSLP 安全漏洞 (CVE-2019-5544) 的权宜措施 (76372)
    解决 vCenter root 密码过期无法登陆 User password expired
    停用Veritas Smart Meter-概念与操作
    centos7.0 没有netstat 和 ifconfig命令问题
    MAVEN剪除传递依赖
    Spring+SpringMVC+MyBatis+Maven框架整合
    CentOS7配置MongoDB
    利用mybatis-generator自动生成代码
    idea提交本地项目到git
  • 原文地址:https://www.cnblogs.com/zhangjianjunab/p/13418881.html
Copyright © 2011-2022 走看看