// 此博文为迁移而来,写于2015年7月14日,不代表本人现在的观点与看法。原始地址:http://blog.sina.com.cn/s/blog_6022c4720102w6ft.html
1、题目
2、TAG
NOIP提高组;二分答案;前缀和。
3、分析
首先题面一定要看懂,看清(上次我就是没有看懂)。说白了,就是使W取一个最合适的值,使每一个所给区间内满足条件的矿石的数量乘上价值之和,最后计算总和使与S最近。
我们根据分值分布来分析算法:
1、30分算法:O(n^3),首先O(n)进行W值选取的枚举,然后O(n^2)判断+计算;
2、50分算法:O(n^2 log n)。根据每次选取W值我们易发现,所得的总值是一个单调函数,随W增大而减小,且一定存在零点。故为了求得最小值,我们需要得到零点位于那两个整数之间。这样用二分显然比枚举快。O(log n)的二分加上O(n^2)的判断+计算;
3、100分算法:200000的数据显然只能O(n log n)或者O(n)了,但是这道题O(n)很明显不现实。我们的任务在于,如何消去O(n^2)这么大的判断与计算?这里,可以采取前缀和优化——由于每次计算检验值之和的时候,W值固定,故可以用O(n)的时间扫描一遍,用a[i]表示前i个矿石满足条件的个数;b[i]表示前i个满足条件的矿石价值之和。故每次询问一个区间的时候,可以O(1)直接相减!这样,复杂度就降到了O(n log n)。(这其实是一个很简单的优化但是我想了很久,看来姿势不够)
4、(这个可以无视)70分算法:存在这个分数段就必然存在这种算法对吧。。。这一定是给那些用了前缀和却没有用二分的。。。虽然我觉得很不可思议。
还有一个很重要的数据范围——权值会爆INT。
4、代码
-----------------------------------------------------------------------------------------------------------------
#include<cstdio>
#include<cstring>
#define MAXN 200005
#define INF 1ll<<60
typedef long long ll;
ll max(ll a,ll b) { return (a>b)?a:b; }
ll min(ll a,ll b) { return (a<b)?a:b; }
ll abs(ll a) { return (a<0)?-a:a; }
ll n,m,w[MAXN],v[MAXN],l[MAXN],r[MAXN],maxw,nl,nr;
ll s,ans=INF,s1[MAXN],s2[MAXN];
void init()
{
scanf("%lld %lld %lld",&n,&m,&s);
for (int i=1;i<=n;i++) { scanf("%lld %lld",&w[i],&v[i]); maxw=max(maxw,w[i]); }
for (int i=1;i<=m;i++) scanf("%lld %lld",&l[i],&r[i]);
}
ll calc(int now)
{
memset(s1,0,sizeof(s1)); memset(s2,0,sizeof(s2));
ll ret=0;
for (int i=1;i<=n;i++) if (w[i]>=now) { s1[i]=s1[i-1]+1; s2[i]=s2[i-1]+v[i]; } else { s1[i]=s1[i-1]; s2[i]=s2[i-1]; }
for (int i=1;i<=m;i++) ret+=(s1[r[i]]-s1[l[i]-1])*(s2[r[i]]-s2[l[i]-1]);
return ret;
}
int main()
{
init();
nl=1,nr=maxw;
while (nl<=nr)
{
int mid=(nl+nr)/2;
ll ret=calc(mid);
if (ret-s《0) { nr=mid-1; ans=min(ans,s-ret); }
else { nl=mid+1; ans=min(ans,ret-s); }
}
printf("%lld",ans);
return 0;
}
---------------------------------------------------------------------------------------------------