【洛谷】分梨子

传送门：点击打开链接

题意：

 

每个梨子都具有两个属性值，Ai和Bi。设在选出的梨子中，两个属性的最小值分别是A0和B0。对于所有被选出的梨子i，都满足C1*(Ai-A0)+C2*(Bi-B0)≤C3（其中，C1、C2和C3都是已知的常数）.

求最多可以挑选出多少个梨子。

数据范围：

 

对于20%的数据，N \le 5N≤5 ，1 \le B_{Min} \le B_{Max} \le 101≤BMin​≤BMax​≤10 。

对于40%的数据，N \le 10N≤10 ，1 \le B_{Min} \le B_{Max} \le 10^61≤BMin​≤BMax​≤106 。

对于100%的数据，N \le 12N≤12 ，0 \le a_i \le 5*10^50≤ai​≤5∗105 ，1 \le B_{Min} \le B_{Max} \le 10^{12}1≤BMin​≤BMax​≤1012 。

题解：

以后还会更新的

一道标签是dp但非常妙的题。先尝试了一下强势枚举，得了20分。代码如下（比较好懂，虽然我也不知道哪出了问题）

#include<cstdio>
#include<cstring>
#include<queue>
#include<vector> 
#include<algorithm>
using namespace std;
int n,a,b,c,ans=1,temp,cnt;
struct q{
    int x,y;
    bool operator <(const q& u)const{
       return x<u.x;
    } 
}t[2005];

struct p{
   int ch,ty;
   bool operator <(const p& u)const{
     return ch<u.ch;
   }
}pr[2005],op;

void update(int k)
{
    cnt=0;
    for(int i=k+1;i<=n;i++){//满足x[k]为x0 
        int qw=a*(t[i].x-t[k].x);
        int pw=b*(t[i].y-t[k].y);
        if(qw>c) break;
        if(pw<=c-qw && pw>=qw-c){
            pr[++cnt].ch=t[i].y-t[k].y;
            pr[cnt].ty=c-qw-pw;
        }
    }//以y[k]为y0 
    sort(pr+1,pr+cnt+1);
    for(int i=1;i<=cnt;i++) if(pr[i].ch>=0){temp=i;break;}
    //找到第一个y[i]>y[k]的数，然后倒序搜索。 
    cnt=cnt-temp+2; 
    ans=max(ans,cnt); //初始化最大值 
    for(int i=temp-1;i>=1;i--){
        cnt++;//取当前y为y0 
        int qw=b*(-pr[i].ch);
        int p=cnt,oppp=cnt; 
        while(pr[cnt].ty<qw){
         cnt--;if(oppp-p<i) p--;
        }
        if(cnt==p) ans=max(ans,p);
     //若满足cnt-i<p表示获得下端的y的大于舍去上端y的个数，取当前y[i]为y0  即有更新答案的可能
    }
}

int main(){
    scanf("%d",&n);
    scanf("%d%d%d",&a,&b,&c);
    for(int i=1;i<=n;i++){
    	scanf("%d%d",&t[i].x,&t[i].y);
    }
    sort(t+1,t+n+1);
    for(int i=1;i<n;i++){
        update(i);
    }
    printf("%d\n",ans);
    return 0;
}

正解如下：

看的是洛谷的题解。传送门：点击打开链接

其实就是优化的暴力，但就是对了，特别妙。

有一点不是很明白：满足x0[i]，y0[j]的枚举时，按理说如果x0本身就不是ij最小值，或y0本身不是ij最小值，那么跳过这种情况，应该时不会对最优解产生影响的，但是加上去（注释的两句）就只能得到一半的分。

 

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,a,b,c,ans,sum[2005];
int x[2005],y[2005];
struct P{
	int v,tag;
	bool operator <(const P& u)const{
		return v<u.v;
	}
	P(int v=0,int tag=0):v(v),tag(tag){}
}t[2005],r[2005];
int main(){
	scanf("%d",&n);
	scanf("%d%d%d",&a,&b,&c);
	for(int i=1;i<=n;i++){
		scanf("%d%d",&x[i],&y[i]);
		t[i]=P(y[i],i);
		r[i]=P(a*x[i]+b*y[i]-c,i);
	}
	sort(t+1,t+n+1);sort(r+1,r+1+n);
	for(int temp=0,i=1;i<=n;i++,temp=0){
		memset(sum,0,sizeof(sum));
		for(int k=0,j=1;j<=n;j++){
			//if(x[i]>x[t[j].tag] || t[j].v>y[i]) continue;
			//???
			for(;k<=n && r[k].v<=a*x[i]+b*t[j].v;k++){
				if(x[r[k].tag]>=x[i] && y[r[k].tag]>=t[j].v){
					temp++;
					sum[y[r[k].tag]]++;
				}
			}
			temp-=sum[t[j-1].v];
			sum[t[j-1].v]=0;
			ans=max(ans,temp);
		}
	}
	printf("%d\n",ans);
	return 0;
}