APIO2007-2015题解大集合（2010年篇）
好消息：这一篇没有上古警告！
特别行动队
经典斜率优化DP。
关于斜率优化dp的式子有两种推法：第一种是按照直线Y=kx+b的形式，把与i有关的看成b，与i,j都相关的看成kx，其中与i相关的是k，与j相关的是x，然后画个图来推；第二种是考虑两个决策点j,k,对它们的决策做差，然后把i有关的移到一边，会形成(f[j]-f[k])/(j-k)>p*i（与i相关的一坨）的形式。其中第一种能够应用的情况比较广泛，因为第二种推法到最后可能会出现移项时不等式符号不好判断的问题。
这个题以第二种推法为例，第一种推法留给读者思考。
设f[i]为所求最大战斗力，则
f[i]=f[j]+a*(s[i]-s[j])^2+b(s[i]-s[j])+c
于是对于j优于k，有
f[j]+a*s[i]^2+a*s[j]^2-2as[i]s[j]+b*s[i]-b*s[j]+c>f[k]+a*s[i]^2+a*s[k]^2-2as[i]s[k]+b*s[i]-b*s[k]+c
f[j]-f[k]>a*(s[k]^2-s[j]^2)+2as[i](s[j]-s[k])+b(s[j]-s[k])
(f[j]+a*s[j]^2-b*s[j])-(f[k]+a*s[k]^2-b*s[j])>2as[i]*(s[j]-s[k])
然后把s[j]-s[k]除过去，就是所求的式子。照着这个式子写斜率优化即可。另外注意斜率用double很慢，有些毒瘤出题人可能会卡几个点必须用分子分母分开
除法
转
乘法
才能过……
#include
#include
#include
#include
#include
#include
#include
#define LL long long
using namespace std;
LL n,a,b,c,head,tail,q[1000005],f[1000005],BEGIN[1000005],s[1000005];
double Slope(LL j,LL k)
{return double(f[k]-f[j]+a*(s[k]*s[k]-s[j]*s[j])-b*(s[k]-s[j]))/double(2*a*(s[k]-s[j]));
}
void Solve()
{f[0]=0,head=1,tail=1;//2as单增
for(LL i=1;i<=n;i++)
{while(headSlope(q[tail],i))tail--;
q[++tail]=i;
}
}
void Input()
{scanf("%lld%lld%lld%lld",&n,&a,&b,&c);
s[0]=0;
for(LL i=1;i<=n;i++)
{scanf("%lld",&BEGIN[i]);
s[i]=s[i-1]+BEGIN[i];
}
}
int main()
{Input();
Solve();
printf("%lld\n",f[n]);
return 0;
}
/*
4
-1 10 -20
2 2 3 4
*/ 巡逻
此题其实是一道NOIP难度的题目……这个题有三种做法：
1、树上的插头DP。这个题目可以转化为求树上不相交的最长链，于是可以用插头乱搞……然而我并不会做……
2、考虑k=1的情况，则只需要找树的最长链即可。考虑k=2，则现在图是一个环套树，仍然需要在图上找出最能节省路程的一条边。这个时候把树的最长链提出来，由于成环的那条路只能走一次，另一次是不能走那条边的，所以还是可以当成一条链来处理。于是有两种情况：
(1)在最长链上某个节点的子树内部进行拼接：直接对最长链上每个点子树进行最长链DP然后取最长+次长即可。
(2)中间要经过一段最长链，易知此时路程应该为节省路程-中间走过的路程长度。写出状态转移方程不难发现这是一个可以用单调栈解决的简单DP，具体方程留给读者自己思考。
3、先求一次最长链，把最长链上长度设为-1，再求一次最长链即可。（代码最后良心地给了几组数据，要的拿走就好……）
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
int cnt=0,ans,wer,n,k,a,b,h[200005]={0},vis[200005],dist[200005],from[200005],stack[200005],D[200005],Maxx[200005],Sec[200005];
int qu[500005];
struct node{int next,to;}edge[400005];
queueq;
void Addedge(int x,int y)
{cnt++,edge[cnt].to=y,edge[cnt].next=h[x],h[x]=cnt;
}
void Input()
{scanf("%d%d",&n,&k);
for(int i=1;iMaxx[x])Sec[x]=Maxx[x],Maxx[x]=Maxx[y]+1;
else if(Maxx[y]+1>Sec[x])Sec[x]=Maxx[y]+1;
}
wer=min(wer,ans-Maxx[x]-Sec[x]+1);
}
void Dp()//链上的DP
{int head=1,tail=1;qu[1]=0;
for(int i=1;i<=stack[0];i++)
{wer=min(wer,ans-(qu[head]+Maxx[stack[i]]-D[stack[i]]-1));
while(head<=tail&&qu[tail]maxx)maxx=dist[i],rec=i;
S=rec;
memset(vis,0,sizeof(vis)),memset(from,0,sizeof(from));
dist[S]=0,from[S]=0;
Dfs(S);
maxx=0;
for(int i=1;i<=n;i++)
if(dist[i]>maxx)maxx=dist[i],rec=i;
ans=2*(n-1)-maxx+1;
wer=999999999;
if(k==1){printf("%d\n",ans);return;}
memset(vis,0,sizeof(vis));
while(!q.empty())q.pop();
int now=rec;
while(now)stack[++stack[0]]=now,vis[now]=1,dist[now]=0,now=from[now];
for(int i=1,j=stack[0];i 信号覆盖
#include
#include
#include
#include
#include
#include
#include
#define LL long long
using namespace std;
LL n,ans1,ans2,top;
struct node{LL x,y;double ang;}point[1505],temp[1505];
void Input()
{scanf("%lld",&n);
if(n==3){printf("3.000000\n");exit(0);}
for(LL i=1;i<=n;i++)
{scanf("%lld%lld",&point[i].x,&point[i].y);
}
}
bool Comp(node x,node y)
{return x.ang=0)
{R=R%top+1,num++;
if(R==i)break;
}
ret-=num*(num-1)/2,num--;
}
return ret;
}
void Solve()
{ans1=0;
for(LL i=1;i<=n;i++)
ans1+=Inner(i);
ans2=n*(n-1)*(n-2)*(n-3)/24-ans1;
double wer=(double(2*ans2+ans1))/double(n*(n-1)*(n-2)/6.0);
printf("%.6lf\n",wer+3);
}
int main()
{Input();
Solve();
return 0;
}