• 【SCOI2007】组队(单调性)


    题目链接

    大意

    给定(N)个人与三个常量(A,B,C),每个人有两个属性:(Hi),(Vi).
    现要让你选些人出来,定义(Hmin)为选出来的这些人中最小的(Hi)值,(Vmin)同理.
    对于选出来的这些人,每个人都要满足$$Acdot Hi+Bcdot Vile Acdot Hmin+Bcdot Vmin+C$$
    求最多能选出多少人。

    (满足(Nle 5000))

    思路

    对于给出的条件,简列一下就是:

    [egin{cases} ①~Acdot Hi+Bcdot Vile Acdot Hmin+Bcdot Vmin+C\ ②~Hige Hmin\ ③~Vige Vmin\ end{cases}]

    首先,易得一个(O(N^3))的算法:暴力枚举(Hmin,Vmin),再对每个人,看是否满足以上条件。

    考虑优化:
    首先观察数据范围,发现支持(O(N^2))的做法。
    那么先随便枚举一个(Hmin)出来(以下推导均已满足(Hminle Hi))。

    考虑在枚举(Vmin)时,如何统计:
    我们设(Si=Acdot Hi+Bcdot Vi),那么对于一个可以被选的点应满足

    [egin{cases} ①Vige Vmin\ ②Sile C+Acdot Hmin+Bcdot Vmin\ end{cases}]

    发现在(Vmin)值变大的时候,式子的右边都是单调递增的。
    所以我们按照(Vi)排序,再枚举一个(Vmin)出来。


    然后对于一个人,我们可以这样想,把它想成二维平面上的一个点((Vi,Si))
    那么就会有两个限制

    [egin{cases} ①Vige Vmin\ ②Sile C+Acdot Hmin+Bcdot Vmin\ end{cases}]

    设有两条直线:$$L1:X=Vmin$$$$L2:Y=C+Acdot Hmin+Bcdot Vmin$$

    那么在(Vmin)值变大的时候,(L1)右移,(L2)上移,如图:

    那么对应的,(L1)经过的点的(Vi)肯定都小于当前的(Vmin),故(L1)经过的都不合法。
    而当前合法的点肯定会被(L2)经过。(包含当前(L2)上点)

    所以,我们就可以在(Vmin)增大时,用什么Vis数组之类的动态维护答案,不断取Max值就行了。

    代码

    #include<cstdio>
    #include<vector>
    #include<algorithm>
    using namespace std;
    #define LL long long
    const int MAXN=5005;
    const int MAXV=10000;
    int N,Ans,Vis[MAXN];
    long long A,B,C;
    struct Node{int H,V,id;LL S;};
    Node s1[MAXN],s2[MAXN],s3[MAXN];
    bool cmpH(Node X,Node Y){return X.H<Y.H;}
    bool cmpV(Node X,Node Y){return X.V<Y.V;}
    bool cmpS(Node X,Node Y){return X.S<Y.S;}
    int main(){
    	//freopen("team.in","r",stdin);
    	//freopen("team.out","w",stdout);
    	scanf("%d%lld%lld%lld",&N,&A,&B,&C);
    	for(int i=1;i<=N;i++){
    		scanf("%d%d",&s1[i].H,&s1[i].V);
    		s1[i].S=s1[i].H*A+s1[i].V*B;s1[i].id=i;
    		s2[i]=s3[i]=s1[i];
    	}
    	sort(s2+1,s2+N+1,cmpV);
    	sort(s3+1,s3+N+1,cmpS);
    	for(int i=1;i<=N;i++){
    		int Hmin=s1[i].H;
    		int p1=1,p2=1,cnt=0;
    		for(int j=1;j<=N;j++)Vis[j]=0;
    		for(int j=1;j<=N;j++){
    			int Vmin=s2[j].V;
    			for(;p1<=N&&s3[p1].S<=C+Hmin*A+Vmin*B;p1++)
    				if(s3[p1].H>=Hmin)cnt+=(Vis[s3[p1].id]==0),Vis[s3[p1].id]=1;
    			for(;p2<=N&&s2[p2].V<Vmin;p2++)
    				if(s2[p2].H>=Hmin)cnt-=Vis[s2[p2].id],Vis[s2[p2].id]=1;
    			Ans=max(Ans,cnt);
    		}
    	}
    	printf("%d
    ",Ans);
    }
    /*
    H*A+V*B<=C+Hmin*A+Vmin*B
    H>=Hmin  V>=Vmin
    */
    
  • 相关阅读:
    强化学习_PolicyGradient(策略梯度)_代码解析
    leetcode_1053. Previous Permutation With One Swap
    leetcode_1052. Grumpy Bookstore Owner
    Tensorflow_入门学习_2_一个神经网络栗子
    通过批处理快速设置IP
    汇编、编译、反汇编、反编译的简单概念介绍
    Logistic回归
    基于概率论的分类方法:朴素贝叶斯
    决策树 预测隐形眼镜类型
    k-近邻算法2(kNN)手写识别系统
  • 原文地址:https://www.cnblogs.com/ftotl/p/11997363.html
Copyright © 2020-2023  润新知