Disillusioning #1 水题+原题赛（被虐瞎）

https://vijos.org/tests/542c04dc17f3ca2064fe7718

好一场 水题 比赛啊

t1直接上暴力费用流10分QAQ，虽然一开始我觉得可以不用的，直接dfs可以得出最大流，但是写撮了就放弃了。

t2直接上暴力又是10分QAQ，虽然本来我就不会。。

t3直接上暴力还写撮了。。。。。。。。sad story

orz 小岛 orz zyf

t1：P1885  Phorni(Disillusioning)

听zyf讲完后无限仰慕。。。。最大流的确可以dfs出来orzQAQ。然后费用流就麻烦了。。

因为是树型的，显然流量是全部流向叶子的，而且叶子到根有且只有一条路径！sigh。。。这样就能保证我选了一条费用最小的路径增广一定是最优！！贪心。。。

sad。

那么我用树剖+线段树维护路径最小，至于官方题解的lct实在是仰慕orz。虽然以前听说过lct优化的网络流，但是今天第一次遇到orz。有时间去看看。。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
using namespace std;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
#define printarr1(a, b) for1(_, 1, b) cout << a[_] << '\t'; cout << endl
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; }
inline const int max(const int &a, const int &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }
#define rdm(i, x) for(int i=ihead[x]; i; i=e[i].next)
#define lc x<<1
#define rc x<<1|1
#define lson l, m, lc
#define rson m+1, r, rc
#define MID (l+r)>>1

const int N=100005, oo=~0u>>1;
int n, ihead[N], cnt, d[N], id[N], fa[N], top[N], son[N], size[N], tot, num, cs[N], upd[N];
struct ED { int to, next, cap, cost; }e[N];
struct TR { int mn, tag; }t[N<<2];
struct ID { int c, id; }a[N];
void add(int u, int v, int c, int w) { e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].to=v; e[cnt].cap=c; e[cnt].cost=w; }
void dfs(int x, int c, int t) {
	d[x]=0; int y; cs[x]=t;
	rdm(i, x) {
		y=e[i].to;
		dfs(y, e[i].cost+c, e[i].cap);
		d[x]+=min(d[y], e[i].cap);
	}
	if(!ihead[x]) d[x]=oo, a[++num].c=c, a[num].id=x;
}
void dfs1(int x) {
	size[x]=1;
	rdm(i, x) {
		int y=e[i].to; fa[y]=x;
		dfs1(y);
		if(size[y]>size[son[x]]) son[x]=y;
		size[x]+=size[y];
	}
}
void dfs2(int x, int t) {
	id[x]=++tot; top[x]=t; upd[tot]=cs[x]; //dbg(x); dbg(tot);
	if(son[x]) dfs2(son[x], t);
	rdm(i, x) if(e[i].to!=son[x]) dfs2(e[i].to, e[i].to);
}
void pushup(int x) { t[x].mn=min(t[lc].mn, t[rc].mn); }
void pushdown(int x) {
	if(t[x].tag) {
		int y=t[x].tag;
		t[x].tag=0;
		t[lc].mn+=y; t[rc].mn+=y;
		t[lc].tag+=y; t[rc].tag+=y;
	}
}
void build(int l, int r, int x) {
	t[x].mn=oo;
	if(l==r) { t[x].mn=upd[l]; return; }
	int m=MID;
	build(lson); build(rson);
	pushup(x);
}
void update(int l, int r, int x, int L, int R, int k) {
	pushdown(x);
	if(L<=l && r<=R) {
		t[x].tag+=k;
		t[x].mn+=k;
		return;
	}
	int m=MID;
	if(L<=m) update(lson, L, R, k); if(m<R) update(rson, L, R, k);
	pushup(x);
}
int query(int l, int r, int x, int L, int R) {
	pushdown(x);
	if(L<=l && r<=R) return t[x].mn;
	int m=MID, ret=oo;
	if(L<=m) ret=query(lson, L, R); if(m<R) ret=min(ret, query(rson, L, R));
	return ret;
}
int getmin(int x) {
	int ret=oo;// dbg(x);
	while(top[x]!=1) {
		ret=min(ret, query(1, n, 1, id[top[x]], id[x]));
		x=fa[top[x]];
	}
	ret=min(ret, query(1, n, 1, 1, id[x])); //dbg(ret);
	return ret;
}
void change(int x, int k) {
	while(top[x]!=1) {
		update(1, n, 1, id[top[x]], id[x], k);
		x=fa[top[x]];
	}
	update(1, n, 1, 1, id[x], k);
}
bool cmp(const ID &a, const ID &b) { return a.c<b.c; }

int main() {
	read(n);
	rep(i, n-1) {
		int u=getint(), v=getint(), c=getint(), w=getint();
		add(u, v, c, w);
	}
	dfs(1, 0, oo);
	dfs1(1); dfs2(1, 1); build(1, n, 1);
	sort(a+1, a+1+num, cmp);
	int flow=d[1], ans=0;
	//for1(i, 1, cnt)
	for1(i, 1, num) {
		int f=min(flow, getmin(a[i].id)); // dbg(f); dbg(a[i].id);
		ans+=f*a[i].c;
		change(a[i].id, -f);
		flow-=f; if(flow==0) break;
	}
	printf("%d\n%d\n", d[1], ans);
	return 0;
}

无限仰膜zyf神犇，千古神犇zyf！