团圆

https://zybuluo.com/ysner/note/1248811

题面

有\(n\)个站台,每个站台距第一个站台的距离为\(s_i\)公里。一列火车，上有\(k\)个座位。\(q\)位乘客，乘坐火车区间为\((l_i,r_i)\)。
设代价为每公里火车上空置座位数之和。
最小化代价，并最大化乘车人数。
（输出代价占\(40\%\)、最大人数乘车方案占\(20\%\)、怎样安排都有车坐的人占\(40\%\)）

• \(30pts\ k=1\)
• \(100pts\ n\leq500,q\leq3500\)

解析

\(30pts\)算法

大力\(DP\)即可。
不过这个\(DP\)看起来好神仙。（通过挂链，把乘客编号挂在右端点上，以找到左端点进行转移）
记得一位\(bitset\)只占\(1/8\)字节。

#include<bits/stdc++.h>
#define ll long long
#define rg register
using namespace std;
const int MAXN=3510;
struct gao{   int Le,Ri;   }App[MAXN];
struct zhang{   int Maxs;   bitset<MAXN>Tra,Mst;   }Ans;
int n,k,q,A[MAXN],Maxv;
struct yi{
    zhang Dp[MAXN];   int First[MAXN],Next[MAXN];
	inline void Solve()
		{   for(rg int i=0;i<=n;i++) First[i]=-1;
			for(rg int i=1;i<=n;i++) A[i]+=A[i-1];
			for(rg int i=1;i<=q;i++)
				Next[i]=First[App[i].Ri],First[App[i].Ri]=i;
		    for(rg int i=1;i<=n;i++)
			{   Dp[i]=Dp[i-1];
				for(rg int j=First[i];j!=-1;j=Next[j])
				{   zhang Now=Dp[App[j].Le];	Now.Tra[j]=1,Now.Mst[j]=1;
					Now.Maxs+=A[App[j].Ri]-A[App[j].Le];
					if(Now.Maxs<Dp[i].Maxs)  continue ;
					if(Now.Maxs==Dp[i].Maxs&&Now.Tra.count()<Dp[i].Tra.count())  continue ;
					if(Now.Maxs>Dp[i].Maxs||Now.Tra.count()>Dp[i].Tra.count()) Dp[i]=Now;
					Dp[i].Mst&=Now.Mst;
				}
			}
			printf("%d\n",Maxv-Dp[n].Maxs);
			for(rg int i=1;i<=q;i++) putchar(Dp[n].Tra[i]?'Y':'N');
			printf("\n%d\n",(int)Dp[n].Mst.count());
			for(rg int i=1;i<=q;i++) if(Dp[n].Mst[i]) printf("%d ",i);
			printf("(%d %d)\n",Dp[2].Maxs,(int)Dp[2].Tra.count());
		}
}TP2;
int main()
{   
	n=Read();   k=Read();
	for(rg int i=2;i<=n;i++) A[i]=Read();
	for(rg int i=n;i>=2;i--) A[i]-=A[i-1],Maxv+=A[i]*k;
	q=Read();
	for(rg int i=1;i<=q;i++)
		App[i].Le=Read(),App[i].Ri=Read();
	if(k==1)  TP2.Solve();
}

\(40\%\)算法

首先这是一个网络流经典模型。
应在各个站台间连容量为\(k\)，费用为\(0\)的边，然后各个区间两端连容量为\(1\)，费用为\(dis[r]-dis[l]\)的边（看作两种不同的选择）。
跑最大费用最大流即可。

跑最大费用与最小费用的区别是，\(dis\)数组初值置为\(-inf\)，然后每次把\(dis\)数组往上更新。
乘车的乘客就是满流的、容量减为\(0\)的带费用边所代表的乘客。

\(60\%\)算法

注意到有一种套路，如果要同时最小化或最大化两个量\(a,b\)，则等价于最小化或最大化\(aM+b\)。
但是，\(M\)必须大到足以区分\(a,b\)。一般来说，\(M>max\{abs(a-b)\}\)。

在当前题目中，\(a\)对应代价，\(b\)对应乘车人数。
于是把带费用边的费用改为\((dis[r]-dis[l])*M+b\)即可。

\(100pts\)算法

针对每个乘客的区间，可以在残量网络上跑一跑最大费用最大流（即找一找有没有不走当前乘客区间的、收益相同的方案）。
若该区间内，收益仍然等于该区间带费用边的费用，就说明不一定要走一开始用的带费用边，且方案仍最优；即可不选该乘客。
反之必选。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#define re register
#define il inline
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int mod=1e9+7,N=8005,inf=1e9,M=1e6;
struct Edge{int to,nxt,w;ll c;}e[N*100];
struct data{int l,r;}a[N];
int n,k,h[N],cnt=1,q,S,T,pv[N],pe[N],gu[N],p[N],top,sta[N];
ll dis[N],ans;
bool vis[N];
queue<int>Q;
il void add(re int u,re int v,re int w,re ll c)
{
  e[++cnt].to=v;e[cnt].nxt=h[u];e[cnt].w=w;e[cnt].c=c;h[u]=cnt;
  e[++cnt].to=u;e[cnt].nxt=h[v];e[cnt].w=0;e[cnt].c=-c;h[v]=cnt;
}
il ll gi()
{
   re ll x=0,t=1;
   re char ch=getchar();https://blog.csdn.net/yypony/article/details/17260153
   while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
   if(ch=='-') t=-1,ch=getchar();
   while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
   return x*t;
}
il int SPFA(re int S)
{
  memset(dis,-63,sizeof(dis));dis[S]=0;vis[S]=1;Q.push(S);
  while(!Q.empty())
    {
      re int u=Q.front();Q.pop();
      for(re int i=h[u];i+1;i=e[i].nxt)
	{
	  re int v=e[i].to;
	  if(e[i].w&&dis[v]<dis[u]+e[i].c)
	    {
	      dis[v]=dis[u]+e[i].c;
	      pe[v]=i;pv[v]=u;
	      if(!vis[v]) vis[v]=1,Q.push(v);
	    }
	}
      vis[u]=0;
    }
  return dis[T]>dis[0];
}
int main()
{
  freopen("reunion.in","r",stdin);
  freopen("reunion.out","w",stdout);
  memset(h,-1,sizeof(h));
  n=gi();k=gi();S=n+1;T=S+1;
  fp(i,2,n) p[i]=gi();
  q=gi();
  fp(i,1,q) a[i].l=gi(),a[i].r=gi(),add(a[i].l,a[i].r,1,(p[a[i].r]-p[a[i].l])*M+1),gu[i]=cnt-1;
  add(S,1,k,0);fp(i,1,n-1) add(i,i+1,k,0);add(n,T,k,0);
  while(SPFA(S))
    {
      re ll sum=1e18;
      for(re int i=T;i!=S;i=pv[i]) sum=min(sum,e[pe[i]].w);
      for(re int i=T;i!=S;i=pv[i]) e[pe[i]].w-=sum,e[pe[i]^1].w+=sum;
      ans+=sum*dis[T];
    }
  printf("%lld\n",p[n]*k-ans/M);
  fp(i,1,q) if(e[gu[i]].w) putchar('N');else putchar('Y');puts("");
  fp(i,1,q)
    {
      if(e[gu[i]].w) continue;
      SPFA(a[i].l);
      if(dis[a[i].r]<(p[a[i].r]-p[a[i].l])*M+1) sta[++top]=i;
    }
  printf("%d\n",top);
  fp(i,1,top) printf("%d ",sta[i]);puts("");
  fclose(stdin);
  fclose(stdout);
  return 0;
}