7045. 2021.04.05【2021省赛模拟】数学考试(sleep)题解

Jose何塞

于 2021-04-06 10:38:31 发布

阅读量149

点赞数

分类专栏：省选组网络流

本文链接：https://blog.csdn.net/wangyuda123456789/article/details/115454557

版权

省选组同时被 2 个专栏收录

13 篇文章 0 订阅

订阅专栏

网络流

3 篇文章 0 订阅

订阅专栏

网络流的好题目
三次函数的形式将我迷惑了，以为是一道数学题
其实这道题的突破口是它奇怪的限制
应该自然联想到网络流！
这里是最小割
不考虑限制，就是对于每个函数每个可能的取值 $x(l\leq x\leq r)$ ，从 $S$ 到 $T$ 连一条链，断掉一条边代表你选了 $f (x)$ 作为这个函数的值,容量为 $f (x)$
那这样岂不是变成了“最大割”问题？
容量取个相反数不就是最小割？
当然，你会连 $r - l + 2$ 条边，值却只有 $r - l + 1$ 个，所以 $S$ 的出边或 $T$ 的入边有一条边容量是 $i n f$ ，我将 $S$ 的出边容量定为 $i n f$ ，那么每个点的出边被断代表选这个点
考虑如何加上限制，其实有个最小割这个思路，限制也不难想。对于一组 $(u, v, d)$ , $u_x向v_{x-d}$ 连一条容量为 $i n f$ 的边, $i n f$ 不可能被割掉，所以选择了 $u_x$ 就必须选择 $v_{x-d}$ 或以后的边割（选 $v_{x-d}$ 前的边无法改变连通性）
如果 $v_{x-d}$ 没值了怎么办？
这就代表 $u_x$ 根本不能选，因为没有对应的 $v_{x-d}$ 满足限制
向 $T$ 连一条 $i n f$ 边即可
然后跑个最小割就能过
实现细节粗节：
你会发现，我们的容量有可能是负数！
加上一个大值 $C$ 使之变成正整数。
答案即为最小割 $- n * C$ （我们知道 $1$ 个函数 $1$ 个取值，也就是一共只会被割 $n$ 条边）
怎么判无解？
当最小割 $\geq inf$ 时,说明你割掉了一条 $i n f$ 边，这是无解的情况
当然 $i n f$ 值要设的足够大

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 700005
#define ll long long
#define C 200000000
#define D 205
#define mem(a) memset(a,0,sizeof(a))
#define S 400001
#define T 400002
#define inf 100000000000
using namespace std;
ll cur[N],gap[N],dis[N],f[N<<1][4],q[N];
ll i,j,k,m,n,o,p,l,s,t,times,k1,k2,k3,k4,last,x,y,d,ans,r;
ll a[105][205],tot;
struct node{
	ll l,r;
}b[10005];
void insert(ll x,ll y,ll z)
{
	if (z!=inf) z=-z+C; 
	f[++t][1]=y,f[t][2]=q[x],f[t][3]=z,q[x]=t;
	f[++t][1]=x,f[t][2]=q[y],f[t][3]=0,q[y]=t;
}
ll dfs(ll t,ll flow)
{
	if (t==T) return flow;
	ll have=0;
	for (ll k=cur[t];k;k=f[k][2])
	{
		if (dis[f[k][1]]+1==dis[t]&&f[k][3])
		{
			cur[t]=k;
			ll now=dfs(f[k][1],min(flow-have,f[k][3]));
			f[k][3]-=now,f[k^1][3]+=now,have+=now;
			if (flow==have) return flow;
		}
	}
	cur[t]=q[t];
	if (!(--gap[dis[t]])) dis[S]=T; 
	++gap[++dis[t]];
	return have;
}
ll cal(ll x) {return k1*x*x*x+k2*x*x+k3*x+k4;}
int main()
{
	freopen("sleep.in","r",stdin);
	freopen("sleep.out","w",stdout);
	scanf("%lld",&times);
	while (times--)
	{
		mem(cur),mem(gap),mem(dis),mem(f),mem(q),t=tot=1;
		scanf("%lld%lld",&n,&m);
		for (i=1;i<=n;i++)
		{
			scanf("%lld%lld%lld%lld%lld%lld",&k1,&k2,&k3,&k4,&l,&r);
			insert(S,++tot,inf);a[i][l+D]=last=tot;
			for (j=l+1,k=2;j<=r;j++,k++)
			{
				insert(last,++tot,cal(j-1));
				a[i][j+D]=last=tot;
			}
			insert(tot,T,cal(r));
			b[i]=(node){l,r};
		}p=0;
		for (i=1;i<=m;i++)
		{
			scanf("%lld%lld%lld",&x,&y,&d);
			o=1e9;
			for (j=b[x].l;j<=b[x].r;j++)
			{
				if (j-d<=b[y].r)
				{
					o=1;
				}
				if (j-d>=b[y].l&&j-d<=b[y].r)
				{
					insert(a[x][j+D],a[y][j-d+D],inf);
				}
				if (j-d>=b[y].r+1) insert(a[x][j+D],T,inf);
			}
			if (o==1e9) 
			{
				p=1;//break;
			}
		}
/*		if (p)
		{
			puts("mei ji ge");
			continue;
		}*/
		ans=0;
		while (dis[S]<T) 
			ans+=dfs(S,inf);
		if (ans>=inf)
		{
			puts("mei ji ge");
			continue;
		}
		printf("%lld\n",n*C-ans);
	}
	return 0;
}