参考资料:icpc暑期课
题意:
N 头牛在长为 L 的线段上吃草(???),每头牛活动范围(S,E),不同奶牛活动范围可以有重叠。现在需要安装喷水头,喷水头可调节,调节范围 [ A,B ],要求:线段上每个整点都恰好位于一个喷水头喷洒范围,每头奶牛的活动范围都要在一个喷水头喷洒范围内,任何喷水头喷水范围不可超过线段两端。现在求解最少安装几个喷头。
问题分析:
1. 从起点往终点安装喷水头,令 f(x) 表示,所安装的喷水头恰好覆盖 [0,x] 最少需要的喷水头
2. 显然x需要满足
① x是偶数,且 x 的位置没有奶牛
②
③ 当 x>2B时,存在,且Y满足上述条件,使得f(x)=f(Y)+1
解题思路:
递推计算 f(x):
- x是奇数,x<2A,x的地方可能有奶牛,f(x) 没得算,不妨定义为 ∞ 。
- ,且 x 位于任何奶牛活动范围之外,f(x) = 1。
- X > 2B:f(x) = 1 + min{ f(Y) } (,且 Y 位于任何奶牛活动范围之外)
关键:快速找到 [x-2B,x-2A] 中使得 f(Y) 最小的元素。(用优先队列去维护)
细节:
- 若坐标属于 [x-2B,x-2A] 的二元组(i,F(i))都保存priority_queue中,并根据 F(i) 排序,则队头元素可以保证是最小的
- 求x点的 F(x) 时,必须确保队列中包含所有属于 [x-2B,x-2A] 的点。而且不允许出现坐标大于 x-2A 的点(这样的点,如果出现在队头,对后续有用,对当前无用,算法无法继续),可以出现小于 x-2B 的点,直接抛弃即可
- 求出 x 点的 F 值后,将(x-2A+2,F(x-2A+2))放入队列,为求 F(x+2) 做准备
- 队列里只需要有偶数点
代码:
#include <cstdio>
#include <queue>
#include <cstring>
#include <iostream>
using namespace std;
const int inf=1<<31;
const int maxl=1001000;
const int maxn=1010;
int f[maxl];//f[L]就是答案
int cow[maxl];//cow[i]表示i点有奶牛
int n,l,a,b;
struct node
{
int f,x;
bool operator <(const node & t)const
{
return f>t.f;
}
node(int xx=0,int ff=0):x(xx),f(ff){}
};//f的值越小,越优先
priority_queue<node> fx;
int main()
{
cin>>n>>l;
cin>>a>>b;
a<<=1;
b<<=1;
memset(cow,0,sizeof(cow));
for(int i=0;i<n;i++)
{
int s,e;
cin>>s>>e;
++cow[s+1];//s+1进入奶牛区
--cow[e];//e退出奶牛区
}
int incow=0;//当前点位于多少个奶牛活动范围内
for(int i=0;i<=l;i++)//每个点是否有奶牛
{
f[i]=inf;
incow+=cow[i];
cow[i]=incow>0;
}
for(int i=a;i<=b;i+=2)
{
if(!cow[i])
{
f[i]=1;
if(i<=b+2-a)
fx.push(node(i,1));
}
}
for(int i=b+2;i<=l;i+=2)
{
if(!cow[i])
{
node tmp;
while(!fx.empty())
{
tmp=fx.top();
if(tmp.x<i-b)
{
fx.pop();
}
else
break;
}
if(!fx.empty())
f[i]=tmp.f+1;
}
if(f[i-a+2]!=inf)
{
fx.push(node(i-a+2,f[i-a+2]));
}
}
if(f[l]==inf)
{
cout<<-1<<endl;
}
else
cout<<f[l]<<endl;
return 0;
}
手写优先队列:
如果一个队列满足一下条件:
- 开始为空
- 每在队尾加入一个元素a之前,都从现有队尾往前删除元素,一直删到碰到小于a的元素为止,然后加入a
那么队列就是递增的,当然队头的元素,一定是队列中最小的。