洛谷 P4360 [CEOI2004]锯木厂选址

题面

---

看了半天才发现其实这不算一道DP

---

$O(n^2)$做法

大意很好理解，就不再赘述了
首先，我们可以肯定的是：锯木厂一定建在树旁边（要不然多跑一段冤枉路）
为了方便思考，我们倒序输入，将山脚下的树定为第一棵树
$f_i$表示将第2个锯木厂建在第$i$棵树旁边的总共最小花费
那么方程就出来了：
$f_i=\sum _{k=i}^n{a}_k(di{s}_k-di{s}_i)+\sum _{k=j}^{i-1}{a}_k(di{s}_k-di{s}_j)+\sum _{k=1}^{j-1}{a}_{k}di{s}_k$
再加上玄学前缀和优化
$s_i=\sum _{j=1}^i{a}_j\ast di{s}_j$
$su{m}_i=\sum _{j=1}^i{a}_j$
$f_i=s_n-s_{i-1}-di{s}_i(su{m}_n-su{m}_{i-1})+s_{i-1}-s_{j-1}-di{s}_j(su{m}_{i-1}-su{m}_{j-1})+s_{j-1}$
整理一下会发现方程还可以化简：
$f_i=s_n-di{s}_i(su{m}_n-su{m}_{i-1})-di{s}_j(su{m}_{i-1}-su{m}_{j-1})$

然后……
$我惊讶地AC了！！！$
$O(n^2)算法！！！$
$水得不像话！！！$

咳咳咳我们还是要来打一打正解的……

优化

老样子
$f_i=s_n-di{s}_i(su{m}_n-su{m}_{i-1})-di{s}_j(su{m}_{i-1}-su{m}_{j-1})$
$f_i=s_n-di{s}_{i}su{m}_n+di{s}_{i}su{m}_{i-1}-di{s}_{j}su{m}_{i-1}+di{s}_{j}su{m}_{j-1}$
$f_i-s_n+di{s}_{i}su{m}_n-di{s}_{i}su{m}_{i-1}=di{s}_j(su{m}_{j-1}-su{m}_{i-1})$
$di{s}_{j}su{m}_{j-1}=di{s}_{j}su{m}_{i-1}+f_i-s_n+di{s}_{i}su{m}_n-di{s}_{i}su{m}_{i-1}$
$y=di{s}_{j}su{m}_{j-1}$
$k=su{m}_{i-1}$
$x=di{s}_j$
$b=f_i-s_n+di{s}_{i}su{m}_n-di{s}_{i}su{m}_{i-1}$
另外$f$数组也可以去掉

代码：

#include<cstdio>
#include<algorithm>
using namespace std;
struct deque
{
/*/
	int list[110];//convinient for debug
/*/
	int *list;//avoid the stack-crashing
//*/

	int head,tail;//head~tail
	
/*/
	deque():head(1),tail(0){}
/*/
	deque():head(1),tail(0),list(new int[1110000]){}//avoid the stack-crashing
//*/
	int size(){return tail-head+1;}
	
	int front(){return list[head];}
	void push_front(int x){list[--head]=x;}
	void pop_front(){++head;}
	int front_2nd(){return list[head+1];}
	
	int back(){return list[tail];}
	void push_back(int x){list[++tail]=x;}
	void pop_back(){--tail;}
	int back_2nd(){return list[tail-1];}
};
#define y(p) (dis[p]*sum[p-1])
#define k(p) (sum[p-1])
#define x(p) (dis[p])
int dis[21000];
long long sum[21000],s[21000];
double slope(int a,int b)
{return (y(a)-y(b))/(double)(x(a)-x(b));}
int main()
{
	int n;
	scanf("%d",&n);
	for(int i=n;i>=1;i--)
		scanf("%lld%d",sum+i,dis+i);
	s[1]=sum[1]*dis[1];
	for(int i=2;i<=n;i++)
		dis[i]+=dis[i-1],s[i]=s[i-1]+sum[i]*dis[i],sum[i]+=sum[i-1];
	long long ans=0x7fffffffffffffffll;
	deque q;
	q.push_back(0);
	for(int i=1;i<=n;i++)
	{
		while(q.size()>1&&slope(q.front(),q.front_2nd())<k(i))q.pop_front();
		int j=q.front();
		ans=min(ans,s[n]-dis[i]*(sum[n]-sum[i-1])-dis[j]*(sum[i-1]-sum[j-1]));
		while(q.size()>1&&slope(q.back_2nd(),q.back())>slope(q.back(),i))q.pop_back();
		q.push_back(i);
	}
	printf("%lld",ans);
	return 0;
}