P4544 [USACO10NOV] Buying Feed G

最新推荐文章于 2024-10-01 23:05:47 发布

内测人员

最新推荐文章于 2024-10-01 23:05:47 发布

阅读量408

点赞数 4

文章标签：算法

本文链接：https://blog.csdn.net/m0_71227570/article/details/136635660

版权

# [USACO10NOV] Buying Feed G

## 题目描述

约翰开车来到镇上，他要带$K$吨饲料回家。运送饲料是需要花钱的，如果他的车上有$X$吨饲料，每公里就要花费$X^2$元，开车D公里就需要$D\times X^2$元。约翰可以从$N$家商店购买饲料，所有商店都在一个坐标轴上，第$i$家店的位置是$X_i$，饲料的售价为每吨$C_i$元，库存为$F_i$。

约翰从坐标$X=0$开始沿坐标轴正方向前进，他家在坐标$X=E$上。为了带$K$吨饲料回家，约翰最少的花费是多少呢？假设所有商店的库存之和不会少于$K$。

举个例子，假设有三家商店，情况如下所示：

|坐标|$X=1$|$X=3$|$X=4$|$E=5$|
|:-:|:-:|:-:|:-:|:-:|
|库存|$1$|$1$|$1$|
|售价|$1$|$2$|$2$|

如果$K=2$，约翰的最优选择是在离家较近的两家商店购买饲料，则花在路上的钱是$1+4=5$，花在商店的钱是$2+2=4$，共需要$9$元。

## 输入格式

第$1$行:三个整数$K,E,N$ 第$2..N+1$行:第$i+1$行的三个整数代表,$X_i,F_i,C_i$.

## 输出格式

一个整数,代表最小花费

## 样例 #1

### 样例输入 #1

```
2 5 3
3 1 2
4 1 2
1 1 1
```

### 样例输出 #1

```
9
```

## 提示

$1 \leq K \leq 10000 , 1 \leq E \leq 500 , 1 \leq N \leq 500$

$0 < Xi < E, 1 \leq Fi \leq 10000, 1 \leq C_i \leq 10^7$

#include<iostream>
#include<stack>
#include<list>
#include<set>
#include<vector>
#include<algorithm>
#include<math.h>
#include<numeric>
#include<map>
#include<cstring>
#include<queue>
#include<iomanip>
#include<cmath>
#include<queue>
#include <bitset>
#include<unordered_map>
	#ifndef local
	#define endl '\n'
#endif */
#define mkp make_pair
using namespace std;
using std::bitset;
typedef long long ll;
typedef long double ld;
const int inf=0x3f3f3f3f;
const ll MAXN=2e6+10;
const ll INF=1e18;
const ll N=2e5+100;
const ll mod=1e9+7;
const ll hash_p1=1610612741;
const ll hash_p2=805306457;
const ll hash_p3=402653189;
//-----------------------------------------------------------------------------------------------------------------*/
// ll head[MAXN],net[MAXN],to[MAXN],edge[MAXN]/*流量*/,cost[MAXN]//费用;
/* 
void add(ll u,ll v,ll c,ll s){
	to[++cnt]=v;net[cnt]=head[u];edge[cnt]=c;cost[cnt]=s;head[u]=cnt;
	to[++cnt]=u;net[cnt]=head[v];edge[cnt]=0;cost[cnt]=-s;head[v]=cnt;
}
struct elemt{
	int p,v;
};
-----------------------------------
求[1,MAXN]组合式和逆元 
ll mi(ll a,ll b){
	ll res=1;
	chile(b){
		if(b%2){
			res=res*a%mod;
		}	
		a=a*a%mod;
		b/=2;
	}
	return res;
}
ll fac[MAXN+10],inv[MAXN+10]
ll C(int m,int n){//组合式C(m,n); 
	if(!n){
		return 1;
	}
	return fac[m]*(inv[n]*inv[m-n]%mod)%mod;
}
fac[0]=1;inv[0]=1;
for(ll i=1;i<=MAXN;i++){
	fac[i]=(fac[i-1]*i)%mod;
	inv[i]=mi(fac[i],mod-2);
}
---------------------------------
 unordered_map<int,int>mp;
//优先队列默认小顶堆 ， greater<int> --小顶堆  less<int> --大顶堆  
priority_queue<elemt,vector<elemt>,comp>q;
struct comp{
	public:
		bool operator()(elemt v1,elemt v2){
			return v1.v<v2.v;
		}
};
	set<int>::iterator it=st.begin();
*/
//emplace_back()  等于push_back(),但效率更高,传输pair时emplace_back(i,j)==push_back({i,j}) 
// vector<vector<int>>edge; 二维虚拟储存坐标 
//-----------------------------------------------------------------------------------------------------------------*/
  //map<int,bool>mp[N]; 
  //emplace_back() 
struct elemt{
	ll x,f,c;
	bool operator<(const elemt v){
		return x<v.x;//按坐标升序
	}
};//把商店定义为结构体，方便排序
//(x+d)^2=x^2+2xd+d^2;
ll dp[510][10010];//当前在第i个商店，已经有j吨饲料的花费
elemt a[510];
int main(){
/*cout<<setiosflags(ios::fixed)<<setprecision(8)<<ans<<endl;//输出ans（float）格式控制为8位小数（不含整数部分）*/
/*cout<<setprecision(8)<<ans<<endl;//输出ans（float）格式控制为8位小数（含整数部分）*/
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);//同步流
ll k,e,n;
cin>>k>>e>>n;
for(int i=1;i<=n;i++){
	cin>>a[i].x>>a[i].f>>a[i].c;
}
sort(a+1,a+n+1);//按坐标升序排序
a[n+1].x=e;
memset(dp,inf,sizeof(dp));
deque<pair<ll,ll>>q;//开两维分别维护dp[i-1][j]-a[i].c*j和k,队列内按dp[i-1][j]-a[i].c*j升序
for(int j=0;j<=a[1].f;j++){//初始化
	dp[1][j]=a[1].c*j+(a[2].x-a[1].x)*j*j;
}
for(int i=2;i<=n;i++){
	while(q.size()){
		q.pop_back();//把之前的清空
	}
	for(int j=0;j<=k;j++){
		while(q.size()&&q.back().first>=dp[i-1][j]-a[i].c*j){
			q.pop_back();//装的货物没有当前多，费用系数还比当前高（不如当前优）的可以踢出去了
		}
		q.push_back({dp[i-1][j]-a[i].c*j,j});//当前点放入队列
		while(q.size()&&q.front().second+a[i].f<j){//该商店全买完也不够j的踢出去
			q.pop_front();
		}
		dp[i][j]=min(dp[i][j],q.front().first+(a[i+1].x-a[i].x)*j*j+a[i].c*j);//最前面的就是最优的，
		//cout<<"i="<<i<<" j="<<j<<" dp="<<dp[i][j]<<" fr="<<q.front().first<<endl;
		
	}
}
//ll ans=dp[n][k]+(e-a[n].x)*k*k;//回到家的花费
cout<<dp[n][k]<<endl;
	return 0;
}