Minimum-cost Flow(最小费用最大流)

题目：https://ac.nowcoder.com/acm/contest/5666/H
思路：
有1e5次询问，如果对于每次询问都用一次最小费用最大流肯定会超时
1.以容量为1建图求最小费用最大流，所求出的每条增广路的流量都是1(因为所有边的容量都是1)。因为所有边的容量都是一样的而且所有边的花费不变，所以所有询问的增广路应该是一样的。
2.对于每个询问容量为u/v,流量为1：cost(u/v,1)=cost(u,v)/v
（容量和流量同时改变一个倍数，所选择的的增广路不变）
3. 每条边的容量变为u,所求的流量变为v，每条增广路的流量变为u。
对于每个询问，从前往后选择增广路，直到所有增广路流量之和v。
v=u*a+b
a表示选择前a条增广路,每条流量为a。
b表示a+1条增广路的流量为b。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#define ll long long
using namespace std;
const int M=210;
const int N=100;
int h[N],pre[N],tot;
int edge[M],cost[M],ne[M],ver[M];
int v[N],dis[N],incf[N];
int s,t,maxflow,ans;
vector<int> st,sum;
ll gcd(ll a,ll b)
{
	return b?gcd(b,a%b):a;
}
void add(int a,int b,int c,int z)
{
	ver[tot]=b;edge[tot]=c;cost[tot]=z;ne[tot]=h[a];h[a]=tot++;
	ver[tot]=a;edge[tot]=0;cost[tot]=-z;ne[tot]=h[b];h[b]=tot++;
}
bool spfa()
{
	memset(v,0,sizeof v);
	memset(dis,0x3f,sizeof dis);
	queue<int> q;
	q.push(s);dis[s]=0;v[s]=1;
	incf[s]=0x3f3f3f3f;
	while(q.size())
	{
		int x=q.front();q.pop();v[x]=0;
		for(int i=h[x];i!=-1;i=ne[i])
		{
			if(!edge[i]) continue;
			int y=ver[i];
			if(dis[y]>dis[x]+cost[i])
			{
				
				dis[y]=dis[x]+cost[i];
				incf[y]=min(incf[x],edge[i]);
				pre[y]=i;
				if(v[y]==0)
				{
					q.push(y);v[y]=1;
				}
			}
		}
	}
	if(dis[t]==0x3f3f3f3f) return 0;
	return 1;
}
void update()
{
	int x=t;
	while(x!=s)
	{
		int i=pre[x];
		edge[i]-=incf[t];
		edge[i^1]+=incf[t];
		x=ver[i^1];
	}
	st.push_back(dis[t]);
	sum.push_back(sum.back()+dis[t]);
	maxflow+=incf[t];
	ans+=dis[t]*incf[t];
}
int main()
{
	int n,m,q;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		tot=0;maxflow=0;ans=0;
		memset(h,-1,sizeof h);
		int a,b,c;
		for(int i=1;i<=m;i++)
		{
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,1,c);
		}
		s=1;t=n;
	    st.clear();
	    sum.clear();
		st.push_back(0); 
		sum.push_back(0);
		while(spfa())
		{
			update();
		}
		scanf("%d",&q);
		for(int i=1;i<=q;i++)
		{
			ll a,b;
			scanf("%lld%lld",&a,&b);
			if(maxflow*a<b)
			{
				printf("NaN\n");continue;	
			}
			ll x=b%a;
			ll y=b/a;
			ll g=gcd(b,sum[y]*a+st[y+1]*x);
			printf("%lld/%lld\n",(sum[y]*a+st[y+1]*x)/g,b/g);
		}
	}
	
}