#include<stdio.h>
#include<string.h>
#include<ctype.h>
#include<math.h>
#include<iostream>
#include<string>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre(){freopen("c://test//input.in","r",stdin);freopen("c://test//output.out","w",stdout);}
#define MS(x,y) memset(x,y,sizeof(x))
#define MC(x,y) memcpy(x,y,sizeof(x))
#define MP(x,y) make_pair(x,y)
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1,class T2>inline void gmax(T1 &a,T2 b){if(b>a)a=b;}
template <class T1,class T2>inline void gmin(T1 &a,T2 b){if(b<a)a=b;}
const int N=1e5+10,M=1e5+10,Z=1e9+7,ms63=1061109567;
int casenum,casei;
int n,m;
int x[N],y[N],cc[N];
int ww[M];
int first[M],id;
int w[N],c[N],nxt[N];
LL f[M];
bool e[M];
void ins(int x,int y,int z)
{
id++;
w[id]=y;
c[id]=z;
nxt[id]=first[x];
first[x]=id;
}
struct node
{
int x;LL v;
node(){}
node(int x_,LL v_){x=x_;v=v_;}
bool operator < (const node& b)const {return v>b.v;}
};
priority_queue<node>q;
void inq(int x,LL v)
{
if(v>=f[x])return;
f[x]=v;
q.push(node(x,v));
}
LL dijkstra()
{
for(int i=1;i<=n;i++)
{
ins(y[i],x[i],cc[i]);
if(ww[y[i]]==0)inq(y[i],0);
}
while(!q.empty())
{
int x=q.top().x;q.pop();
if(e[x])continue;e[x]=1;
for(int z=first[x];z;z=nxt[z])inq(w[z],f[x]+c[z]);
}
LL ans=0;
for(int i=1;i<=m;i++)
{
if(ww[i]==2)
{
if(f[i]==1e12)return -1;
ans+=f[i];
}
}
return ans;
}
const int L=1e6;
int Q[L],h,t;
void inQ(int x,LL v)
{
if(v>=f[x])return;
f[x]=v;
if(e[x])return;
e[x]=1;
Q[t++]=x;
}
LL spfa()
{
h=t=0;
for(int i=1;i<=n;i++)
{
ins(y[i],x[i],cc[i]);
if(ww[y[i]]==0)inQ(y[i],0);
}
while(h<t)
{
int x=Q[h++];e[x]=0;
for(int z=first[x];z;z=nxt[z])inQ(w[z],f[x]+c[z]);
}
LL ans=0;
for(int i=1;i<=m;i++)
{
if(ww[i]==2)
{
if(f[i]==1e12)return -1;
ans+=f[i];
}
}
return ans;
}
int main()
{
scanf("%d",&casenum);
for(casei=1;casei<=casenum;casei++)
{
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
first[i]=0;
f[i]=1e12;
e[i]=0;
}id=0;
for(int i=1;i<=n;i++)scanf("%d",&x[i]);
for(int i=1;i<=n;i++)scanf("%d",&y[i]);
for(int i=1;i<=n;i++)scanf("%d",&cc[i]);
for(int i=1;i<=m;i++)scanf("%d",&ww[i]);
//首先是dijkstra写法
//printf("Case #%d: %lld\n",casei,dijkstra());
//然后是SPFA写法
printf("Case #%d: %lld\n",casei,spfa());
}
return 0;
}
/*
【trick&&吐槽】
有1e5个点,边权也为1e5,是可能爆int的。
【题意】
有T(30)组数据。
对于每组数据有n(1e5)个村庄,m(1e5)个战场
对于村庄i,曹操可以选择支付c[i]*num元,(0<=c[i]<=1e5)
使得这个村庄派出num个士兵,在战场x[i]为曹操作战,在战场y[i]为袁绍作战。(1<=x[i],y[i]<=m)
对于每个战场,都有个战略价值(0,1,2)。
在战略价值为2的战场,曹操在该战场的士兵数必须严格比袁绍多
在战略价值为1的战场,曹操在该战场的士兵数必须不能比袁绍少(按照贪心原则,实际一定会相等)
在战略价值为0的战场,曹操在该战场的士兵数无所谓(按照贪心原则,实际一定会为0)
让你输出保障上述条件曹操至少需要花费的金钱。
如果无法保障这个条件,则输出-1。
【类型】
差分约束思想
网络流思想
最短路。
【分析】
首先有一个很显然的贪心——
如果一个战场的战略价值为2,那么在这个战场,恰有:曹兵-袁兵=1
如果一个战场的战略价值为1,那么在这个战场,恰有:曹兵=袁兵=0
这题我们可以从村庄向两个战场连边。
但是——如何更加简单地设置关系呢?
为何不考虑直接从曹战场向袁战场连边呢?
这题设计到分配,尽管数据巨大不能用网络流做,我们还是可以考虑先简化数据规模,引入思想。
我们对于每个村庄所提供的(x,y,z),从曹战场x向袁战场y连接一条流量无穷,费用为z的边。
然后,我们有一个关键性的问题就是,如何保证题目的要求呢?
首先我们简化问题,只去考虑合法性。
对于所有战略价值为2的战场,如果从这个样的点ST出发,沿着图,能够达到任意一个战略价值为0的点。
那说明,这个战略价值为2的战场,可以使得其上的曹兵比袁兵多。
为什么能说明这个呢?
很显然,我们对于一条边,使得ST的曹兵+1,下一个点的袁兵+1。
但是这样,对于下一个点是有影响的。
如果下个点重要度为0,显然这样已经可以了。
如果下个点重要度为1,我们还要继续把这1个人转移出去,一直转移到重要度为0的点上即可。
如果下个点重要度为2,我们依然一定要把这1个人转移出去。而只要转移出去了,下个点实际上其实并没有受到影响。
于是,我们枚举所有重要度为2的点,每个点沿着这个图转移一个人到重要度为0的战场上,这道题就可以AC了。
体现在网络流上,就是:
[超级源点->重要度为2的点,流量为1,费用为0]
[曹战场->袁战场,流量无限,费用为人的雇佣成本]
[重要度为0的点->超级汇点,流量无限,费用为0]
跑一个最小费用最大流,如果最大流=重要度为2的点数,那么这个最小费用就是答案了。
只是时间复杂度不允许我们跑最小费用最大流。
然而我们发现,基于这道题的特殊性,其实直接求:
枚举每个重要度为2的点,
对于每个这样的点,求它到重要度为0的点中距离最近的那个的距离。
然后这些距离求和即可。
但是这样是多源最短路。
然而我们发现,只要把边逆过来,初始化所有重要度为0的点为起点,都归为ZERO,求最短路。
然后累加所有重要度为2的点的到 ZERO的最短路即可。
【时间复杂度&&优化】
dijkstra O(nlogn)
SPFA O(kn)
*/
【HDU5545 2015 CCPC 南阳国赛F】【差分约束思想 费用流思想】The Battle of Guandu 官渡之战 重要战场人数多 士兵流向转化问题为最短路spfa+dijkstra双写
最新推荐文章于 2019-09-14 11:55:45 发布