zoukankan      html  css  js  c++  java
  • 【洛谷P5785】[SDOI2012]任务安排

    【题目传送门】

    题解

    首先设计出 (n^2) 的 DP,(dp(i,j)) 表示前 (i) 个分成 (j) 段。
    后来发现空间不够,于是第二维可以压掉,反正每一次都是从前一维状态转移来。
    由于 (T) 的范围,新加入的斜率不一定单调递增,所以单调队列每次不能弹出队首元素,变成了单调栈。
    仍然是维护斜率单调递增的下凸包,但是转移的时候二分找到最优的转移点。

    代码

    #include<bits/stdc++.h>
    using namespace std;
    #define ll long long
    const int INF = 0x3f3f3f3f,N = 3e5+10;
    inline ll read()
    {
    	ll ret=0;char ch=' ',c=getchar();
    	while(!(c>='0'&&c<='9')) ch=c,c=getchar();
    	while(c>='0'&&c<='9') ret=(ret<<1)+(ret<<3)+c-'0',c=getchar();
    	return ch=='-'?-ret:ret;
    }
    ll dp[N],tmp[N];
    ll s,n;
    ll sumT[N],sumC[N];
    int q[N],l,r;
    inline double X(int i){return sumC[i];}
    inline double Y(int i){return dp[i]-sumC[i]*s;}
    inline int find(double t)
    {
    	int L=l,R=r;
    	if(l==r) return q[l];//这样特判比较简洁
    	while(L<R)
    	{
    		int mid=(L+R+1)>>1;
    		if((Y(q[mid])-Y(q[mid-1]))<=(X(q[mid])-X(q[mid-1]))*t) L=mid;
    		else R=mid-1;
    	}
    	return q[L];//注意二分的是q数组下标
    }
    int main()
    {
    	n=read(),s=read();
    	for(int i=1;i<=n;i++) 
    	{
    		sumT[i]=read(),sumT[i]+=sumT[i-1];
    		sumC[i]=read(),sumC[i]+=sumC[i-1];
    	}
    	memset(dp,0x3f,sizeof(dp));//和暴力一样的初始化
    	dp[0]=0;
    	l=1,r=1;
    	for(int i=1;i<=n;i++)
    	{
    		int j=find(sumT[i]);//二分找到转移点
    		dp[i]=Y(j)+sumT[i]*(sumC[i]-sumC[j])+sumC[n]*s;
    		while(l<r&&(X(i)-X(q[r-1]))*(Y(q[r])-Y(q[r-1]))>=(X(q[r])-X(q[r-1]))*(Y(i)-Y(q[r-1]))) r--;
    		q[++r]=i; 
    	}
    	printf("%lld
    ",dp[n]);
    	return 0;
    }
    
  • 相关阅读:
    ESAPI = Enterprise Security API
    WingIDE中调试GAE(google app engine)
    C# 发送Http请求 WebClient类
    [转]使用Google App Engine Helper for Django
    装上Window7
    最近遇到个关于接口的奇怪的问题
    JNDI概述(转载)
    Google App Engine 中通过自定义Django的filter解决时区问题
    C# string与byte[]互转
    Python天天美味(33) 五分钟理解元类(Metaclasses)[转]
  • 原文地址:https://www.cnblogs.com/conprour/p/15480825.html
Copyright © 2011-2022 走看看