poj 2373 Diving the Path

dp方程为dp[i]=min{dp[j]}+1,2*a<=i-j<=2*b，表示一个喷头控制范围的右端点为i时，最少需要的喷头数（i肯定为偶数)

直接搞转移是L^2的，但j是在由i确定的一个范围内的，所以可以用单调队列优化

另外，这题比较麻烦的是如果i在某头牛的区间内，i处就不能放喷头，我们需要知道哪些地方不能放置，可以对区间先预处理一下，能合并的就合并，然后用一个指针p指向当前区间，如果i大于等于区间右端点，指针右移，如果i大于区间左端点，不能放置

另外需要注意的是队头元素可能是inf，加1后大于inf，最后判断无解时，应判断dp[l]是否是大于等于inf,在这wa了，调了别人的代码发现dp[i]=min(dp[i],pos[head]+1)可以AC，若写成dp[i]=pos[head]+1就WA，才发现这个问题

代码：

#include<iostream>
#include<memory.h>
#include<string>
#include<cstdio>
#include<algorithm>
#include<math.h>
#include<stack>
#include<queue>
#include<vector>
#include<map>
#include<ctime>
using namespace std;
const int MAX=1000005;
const int inf=1<<30;
struct node 
{
	int l,r;
}cow[1005];
int flag[MAX],dp[MAX],pos[MAX],id[MAX];
int cmp(node a,node b)
{
	if(a.l!=b.l)
		return a.l<b.l;
	else
		return a.r<b.r;
}
int main()
{
	int i,j,k,n,l,p,a,b,head,tail;
	while(scanf("%d%d",&n,&l)!=EOF)
	{
		scanf("%d%d",&a,&b);
		for(i=0;i<n;i++)
		{
			scanf("%d%d",&cow[i].l,&cow[i].r);
		}
		sort(cow,cow+n,cmp);
		j=0;
		for(i=0;i<n;i++)
		{
			if(i==0)
			{
				cow[j].l=cow[i].l; cow[j++].r=cow[i].r;
			}
			else
			{
				if((cow[j-1].l<=cow[i].l&&cow[j-1].r>cow[i].l))
				{
					cow[j-1].l=min(cow[j-1].l,cow[i].l); cow[j-1].r=max(cow[i].r,cow[j-1].r);
				}
				else
				{
					cow[j].l=cow[i].l; cow[j++].r=cow[i].r;
				}
			}
		}
		n=j;
		for(i=1;i<=l;i++)
			dp[i]=inf;
		dp[0]=0;
		head=tail=0;
		pos[0]=0;
		id[0]=0;
		p=0;
		for(i=2;i<=l;i+=2)
		{
			if(p<n&&cow[p].r<=i)
				p++;
			if(p<n&&cow[p].l<i)
				continue;
			while(head<=tail&&(i-id[head]>2*b))
			{
				head++;
			}
			if(head>tail||i-id[head]<2*a)
				dp[i]=inf;
			else
				dp[i]=pos[head]+1;//pos[head]有可能大于等于inf 
			while(head<=tail&&pos[tail]>dp[i])
					tail--;
			pos[++tail]=dp[i];
			id[tail]=i;
		}
		if(dp[l]>=inf)
			dp[l]=-1;
		printf("%d\n",dp[l]);
	}
	return 0;
}