zoj-3726-Alice's Print Service【二分+rmq】

标签：

答案就是 q*p[i]（i为最大的s[i]<=q的位置，二分），s[i+1]*p[i+1],....s[n]*p[n] 之中最小的

要快速的查找 (i,n)区间的s*p查询手段有很多。这里我用的rmq。

/*
就是一个打印分段收费政策，印的越多，单张价格越低，
输入需要印刷的数量，求最小印刷费用

打印k页的资料，给出n中付费方案，
一次打印超过s1但不超过s2的每页收费p1，
超过s2不超过s3的收费p2.....数据保证
0=s1<s2<...<sn，p1>=p1>=p3>=...>=pn。
接下来m个查询，对于每个查询问最少花多少钱？
例如s1=0 s2=100  p1=20 p2=10 的时候，若要打印99页，
显然直接打印100页要更便宜一点..所以结果是1000..
*/
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <algorithm>
#include <math.h>
#define fmin(a,b) a>b?b:a 
using namespace std;
typedef long long int ll;
const int maxs=100009;
const int max_log=18;
ll s[maxs],p[maxs];
ll mn[maxs][max_log];
void init(ll n)
{
	for(ll i=1;i<=n;++i)
		mn[i][0]=s[i]*p[i];
	ll k=(ll)(log(n*1.0)/log(2.0));
	for(ll j=1;j<=k;++j)
		for(ll i=1;i<=n+1-(1<<j);++i)
			if(i+(1<<j)-1<=n)
				mn[i][j]=fmin(mn[i][j-1],mn[i+(1<<(j-1))][j-1]);		
}
ll get(ll l,ll r)
{
	ll k=(ll)(log((r-l+1)*1.0)/log(2.0));  
	return fmin(mn[l][k],mn[r-(1<<k)+1][k]);
}
void solve(ll q,ll n)
{
	ll rb=lower_bound(s+1,s+n+1,q)-s;
//	printf("%lld\n",rb);
	if(s[rb]>q||rb>n) --rb;
	ll res=q*p[rb];
	if(rb+1<=n)
		res=fmin(res,get(rb+1,n));
	printf("%lld\n",res);
}
int main()
{
	ll t; 
	ll n,m;
	scanf("%lld",&t);
	while(t--)
	{
		scanf("%lld%lld",&n,&m);
		init(n);
		for(ll i=1;i<=n;++i)
			scanf("%lld%lld",s+i,p+i);
		init(n);
		ll q;
		while(m--)
		{
			scanf("%lld",&q);
			solve(q,n);	
		}
	}
	return 0;
} 

zoj-3726-Alice's Print Service【二分+rmq】

标签：

原文地址：http://blog.csdn.net/a915800048/article/details/51366910