題目鏈接

見上題。
每類志愿者可能是若干段，不滿足那個...全幺模矩陣(全單位模矩陣)的條件，所以線性規劃可能存在非整數解。
于是就可以用費用流水過去順便拿個rank2 233.

//20704kb 300ms #include <queue> #include <cstdio> #include <cctype> #include <cstring> #include <algorithm> //#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<int> 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]<INF; } int MCMF() {int res=0, mn=INF;for(int i=des; i!=src; i=fr[pre[i]])mn=std::min(mn,cap[pre[i]]);for(int i=des,v=pre[i]; i!=src; i=fr[v],v=pre[i])res+=mn*cost[v], cap[v]-=mn, cap[v^1]+=mn;return res; }int main() {n=read(), m=read(), Enum=1, src=1, des=n+2;for(int i=1; i<=n; ++i) AddEdge(i,i+1,INF-read(),0);for(int i=1,k,c; i<=m; ++i){k=read();for(int j=1; j<=k; ++j) S[j]=read(),T[j]=read();c=read();for(int j=1; j<=k; ++j) AddEdge(S[j],T[j]+1,INF,c);}AddEdge(n+1,des,INF,0);int res=0;while(SPFA()) res+=MCMF();printf("%d\n",res);return 0; }

轉載于:https://www.cnblogs.com/SovietPower/p/9364197.html

總結

以上是生活随笔為你收集整理的BZOJ.3265.志愿者招募加强版(费用流SPFA)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。