。
每类志愿者可能是若干段,不满足那个...全幺模矩阵(全单位模矩阵)的条件,所以线性规划可能存在非整数解。 于是就可以用费用流水过去顺便拿个rank2 233.//20704kb 300ms#include#include #include #include #include //#define gc() getchar()#define MAXIN 300000#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)const int N=1007,M=1000005,INF=0x3f3f3f3f;int src,des,n,m,S[N],T[N],Enum,H[N],nxt[M],fr[M],to[M],cap[M],cost[M],pre[N],dis[N];char IN[MAXIN],*SS=IN,*TT=IN;inline int read(){ int now=0;register char c=gc(); for(;!isdigit(c);c=gc()); for(;isdigit(c);now=now*10+c-'0',c=gc()); return now;}inline void AddEdge(int u,int v,int w,int c){ fr[++Enum]=u, to[Enum]=v, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w, cost[Enum]=c; fr[++Enum]=v, to[Enum]=u, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0, cost[Enum]=-c;}bool SPFA(){ static std::queue q; static bool inq[N]; memset(dis,0x3f,sizeof dis); dis[src]=0, q.push(src); while(!q.empty()) { int x=q.front(); q.pop(); inq[x]=0; for(int v,i=H[x]; i; i=nxt[i]) if(cap[i] && dis[v=to[i]]>dis[x]+cost[i]) { dis[v]=dis[x]+cost[i], pre[v]=i; if(!inq[v]) inq[v]=1, q.push(v); } } return dis[des]