标签：Land const int ll maxn dp USACO08MAR include Acquisition

嘟嘟嘟


只要会决策单调性，这题就是练手的


首先按矩形长排序，这样只用考虑宽了。
然后很容易搞出dp方程

\[dp[i] = min _ {j = 0} ^ {i - 1} (dp[j] + x[i] * max_{k = j + 1} ^ {i} y[k]) \]

找max可以用st表达到\(O(1)\)。
打表发现决策单调。然后就是正常的优化了。
二分的时候需要注意当前队列非空。要不然会像我一样，不开氧气AC,开了RE2个点。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
#define enter puts("") 
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define rg register
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 5e4 + 5;
inline ll read()
{
  ll ans = 0;
  char ch = getchar(), last = ' ';
  while(!isdigit(ch)) last = ch, ch = getchar();
  while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - '0', ch = getchar();
  if(last == '-') ans = -ans;
  return ans;
}
inline void write(ll x)
{
  if(x < 0) x = -x, putchar('-');
  if(x >= 10) write(x / 10);
  putchar(x % 10 + '0');
}

int n;
ll dp[maxn];
struct Node
{
  ll x, y;
  bool operator < (const Node& oth)const
  {
    return x < oth.x || (x == oth.x && y < oth.y);
  }
}t[maxn];
struct Node2
{
  int pos, L, R;
}q[maxn];
int l = 1, r = 0;

ll Max[20][maxn], b[maxn];
void init()
{
  for(int i = 1; i <= n; ++i) Max[0][i] = t[i].y;
  for(int j = 1; (1 << j) <= n; ++j)
    for(int i = 1; i + (1 << j) - 1 <= n; ++i)
      Max[j][i] = max(Max[j - 1][i], Max[j - 1][i + (1 << (j - 1))]);
  int x = 0;
  for(int i = 1; i <= n; ++i)
    {
      if((1 << (x + 1)) <= i) x++;
      b[i] = x;
    }
}
ll query(int L, int R)
{
  int k = b[R - L + 1];
  return max(Max[k][L], Max[k][R - (1 << k) + 1]);
}

ll w(int L, int R)
{
  return query(L, R) * t[R].x;
}

int solve(int x, Node2 a)
{
  int L = a.L, R = a.R;
  while(L <= R)
    {
      int mid = (L + R) >> 1;
      if(dp[x] + w(x + 1, mid) <= dp[a.pos] + w(a.pos + 1, mid))
	{
	  if(R != mid) R = mid;
	  else {L = mid; break;}
	}
      else
	{
	  if(L != mid + 1) L = mid + 1;
	  else break;
	}
    }
  return L;
}

int main()
{
  n = read();
  for(int i = 1; i <= n; ++i) t[i].x = read(), t[i].y = read();
  sort(t + 1, t + n + 1);
  init();
  q[++r] = (Node2){0, 1, n};
  for(int i = 1; i <= n; ++i)
    {
      while(q[l].R < i) l++;
      dp[i] = dp[q[l].pos] + w(q[l].pos + 1, i);
      q[l].L = i + 1;
      while(l <= r && dp[i] + w(i + 1, q[r].L) <= dp[q[r].pos] + w(q[r].pos + 1, q[r].L)) r--;
      int pos = i;
      if(l <= r) pos = solve(i, q[r]), q[r].R = pos - 1;
      if(pos <= n) q[++r] = (Node2){i, pos, n};
    }
  write(dp[n]), enter;
  return 0;
}

标签：Land,const,int,ll,maxn,dp,USACO08MAR,include,Acquisition
来源： https://blog.51cto.com/u_15234622/2831198

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。