bzoj1597：dp+斜率优化

标签：

日。。。。 一开始没有把他倒过来然后一直乱一直乱，没办法了就重新写，并且把他倒过来，然后就和做过的差不多了。。。

----------------------------------------------------------------------------------------

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define rep(i,n) for(int i=1;i<=n;i++)
#define down(i,n) for(int i=n;i;i--)
#define ll long long
#define clr(x,c) memset(x,c,sizeof(x))
int read(){
 int x=0,f=1;char c=getchar();
 while(!isdigit(c)){
  if(c==‘-‘) f=-1;
  c=getchar();
 }
 while(isdigit(c)){
  x=x*10+c-‘0‘;
  c=getchar();
 }
 return x*f;
}
struct data{
 int x,y;
 bool operator<(const data&rhs)const {
   return x<rhs.x||x==rhs.x&&y<rhs.y;}
};
data e[50005];
ll f[50005];
int q[50005],t[50005];
double pd(int i,int j){
 return (double)(f[i]-f[j])/(e[t[j+1]].y-e[t[i+1]].y);
}
int main(){
 int n=read();
 rep(i,n)  {
  e[i].x=read(),e[i].y=read();
 }
 sort(e+1,e+n+1);
 int cur=0,nmax=0;
 down(i,n){
  if(e[i].y>nmax){
   nmax=e[i].y;
   t[++cur]=i;
  }
 }
 for(int l=1,r=cur;l<r;l++,r--) swap(t[l],t[r]);//就是这里，看了师兄的才知道，错的时候总以为是什么新的技巧。。。然并卵
 /*rep(i,cur){
  printf("%d %d\n",e[t[i]].x,e[t[i]].y);
 }*/
 n=cur;
 int ql=0,qr=1;q[0]=0;
 rep(i,n){
  while(ql+1<qr&&pd(q[ql+1],q[ql])<e[t[i]].x) ql++;
  int tmp=q[ql];
  f[i]=f[tmp]+(ll)(e[t[tmp+1]].y)*e[t[i]].x;
  while(ql+1<qr&&pd(q[qr-1],q[qr-2])>=pd(i,q[qr-1])) qr--;
  q[qr++]=i;
 }
 printf("%lld\n",f[n]);
 return 0;
}

----------------------------------------------------------------------------------------

Description

Input

Output

Sample Input

Sample Output

HINT

Source

bzoj1597：dp+斜率优化

标签：

原文地址：http://www.cnblogs.com/fighting-to-the-end/p/5096993.html