Loj#3026-「ROIR 2018 Day1」管道监控【Trie,费用流】

正题

题目链接:https://loj.ac/p/3026


题目大意

给出nnn个点的一棵外向树,然后mmm个字符串和费用表示你每次可以花费这个费用覆盖路径字符串和给出字符串相等的路径,求覆盖所有边的最小花费(可以重复覆盖)

输出方案

1≤n≤500,1≤m≤105,∑∣S∣≤1061\leq n\leq 500,1\leq m\leq10^5,\sum |S|\leq 10^61n500,1m105,S106


解题思路

注意到nnn很小,可以考虑枚举计算所有路径的最小花费,先构一个TrieTrieTrie即可处理这部分。

然后考虑怎么覆盖的问题,和[NOI2008] 志愿者招募类似的,我们可以通过边调走流来实现覆盖问题。

先原点向每一个叶子连流量为1的边,设sizxsiz_xsizx表示xxx的子树中有多少个叶子那么xxxfaxfa_xfax连接流量为sizx−1siz_x-1sizx1的边,然后根向汇点连流量为siz1siz_1siz1的边。

这样就保证了每条边都至少有一个流量被调走,但是需要考虑单链上的重复覆盖所以还需要每个点向儿子连无向流量的边。


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define ll long long
using namespace std;
const ll N=510,M=1e6+10;
struct node{ll to,next,w;
}a[N];
ll n,m,t,tot,ans,MF,ls[N],siz[N],fa[N];
ll cnt,ch[M][26],cost[M],fail[M],mark[M];
queue<int> q;char s[M];
void addl(ll x,ll y,ll w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return;
}
void Insert(char *s,ll val,ll num){ll x=1,m=strlen(s);for(ll i=0;i<m;i++){ll c=s[i]-'a';if(!ch[x][c])ch[x][c]=++cnt;x=ch[x][c];}if(cost[x]>val)cost[x]=val,mark[x]=num;return;
}
void GetFail(){for(ll i=0;i<26;i++)q.push(ch[1][i]);while(!q.empty()){ll x=q.front();for(ll i=0;i<26;i++){if(!ch[x][i])ch[x][i]=ch[fail[x]][i];else{fail[ch[x][i]]=ch[fail[x]][i];q.push(ch[x][i]);}}}return;
}
namespace NF{struct node{ll to,next,w,c,typ;}a[N*N];ll tot=1,s,t,ls[N],f[N],mf[N],pre[N];bool v[N];void addl(ll x,ll y,ll w,ll c,ll mk){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;a[tot].c=c;//a[tot].typ=mk;a[++tot].to=x;a[tot].next=ls[y];ls[y]=tot;a[tot].w=0;a[tot].c=-c;a[tot].typ=mk;return;}bool SPFA(){memset(f,0x3f,sizeof(f));q.push(s);v[s]=1;f[s]=0;mf[s]=f[0];mf[t]=0;while(!q.empty()){ll x=q.front();v[x]=0;q.pop();for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(a[i].w&&f[x]+a[i].c<f[y]){f[y]=f[x]+a[i].c;pre[y]=i;mf[y]=min(mf[x],a[i].w);if(!v[y])v[y]=1,q.push(y);}}}return mf[t];}void Updata(){ll x=t;ans+=mf[t]*f[t];MF+=mf[t];while(x!=s){a[pre[x]].w-=mf[t];a[pre[x]^1].w+=mf[t];x=a[pre[x]^1].to;}return;}void GetAns(){while(SPFA())Updata();return;}void Print(){ll sum=0;for(ll i=2;i<=tot;i++)if(a[i].typ)sum+=a[i].w;printf("%lld\n",sum);for(ll i=2;i<=tot;i++)if(a[i].typ){while(a[i].w)printf("%lld %lld %lld\n",a[i^1].to,a[i].to,a[i].typ),a[i].w--;}}
}
void calc(ll x,ll p,ll s){if(!p)return;if(cost[p]<=1e9)NF::addl(x,s,n,cost[p],mark[p]);for(ll i=ls[x];i;i=a[i].next)calc(a[i].to,ch[p][a[i].w],s);return;
}
void solve(ll x){siz[x]+=(!ls[x]);if(!ls[x])NF::addl(NF::s,x,1,0,0);calc(x,1,x);for(ll i=ls[x];i;i=a[i].next){NF::addl(x,a[i].to,n,0,0);solve(a[i].to),siz[x]+=siz[a[i].to];}if(fa[x]&&siz[x]>1)NF::addl(x,fa[x],siz[x]-1,0,0);return;
}
signed main()
{memset(cost,0x3f,sizeof(cost));scanf("%lld%lld%lld",&n,&m,&t);for(ll i=2;i<=n;i++){ll x;char op[3];scanf("%lld%s",&fa[i],op);addl(fa[i],i,op[0]-'a');}cnt=1;for(ll i=1;i<=m;i++){ll c;scanf("%lld%s",&c,s);Insert(s,c,i);}NF::s=n+1;NF::t=n+2;solve(1);NF::addl(1,NF::t,siz[1],0,0);NF::GetAns();if(MF!=siz[1])return puts("-1")&0;printf("%lld\n",ans);if(t)NF::Print();return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318630.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

不止代码:区间圆数(ybtoj-数位DP)

文章目录题目描述解析理解万岁&#xff01;代码题目描述 解析 一寸山河一寸血 理解万岁&#xff01; 首先&#xff0c;这题统计[l,r]的个数&#xff0c;可以用[1,r]-[1,l-1]来实现 接下来就是如何统计出[1,n]的个数了 首先&#xff0c;用dp[pos][s0]来表示一个二进制pos位有s…

P3224 [HNOI2012]永无乡(并查集+权值线段树合并/平衡树)

[HNOI2012]永无乡 Code1 权值线段树天然支持merge&#xff0c;线段树上二分求第k小 #include<bits/stdc.h>using namespace std; using lllong long; template <class Tint> T rd() {T res0;T fg1;char chgetchar();while(!isdigit(ch)) {if(ch-) fg-1;chgetcha…

牛客题霸 [没有重复项数字的所有排列] C++题解/答案

牛客题霸 [没有重复项数字的所有排列] C题解/答案 题目描述 给出一组数字&#xff0c;返回该组数字的所有排列 例如&#xff1a; [1,2,3]的所有排列如下 [1,2,3],[1,3,2],[2,1,3],[2,3,1],[3,1,2], [3,2,1]. &#xff08;以数字在数组中的位置靠前为优先级&#xff0c;按字典…

.NET Core 2.2 新增部分功能使用尝鲜

前言美国当地时间12月4日&#xff0c;微软2019开发者大会中发布了一系列的重磅消息&#xff0c;包含了软硬件和开源社区的各种好消息是铺天盖地&#xff0c;作为一名普通的开发者&#xff0c;我第一时间下载了 .NET Core 2.2 的源码&#xff0c;针对发布说明逐条浏览&#xff0…

不止代码:机器分配(动态规划)

题目描述 解析 头疼 什么破题 就是一个dp寻找最优性决策的常规题 但是要输出过程&#xff0c;可以使用递归输出 但是&#xff01; 这题数据的意思是&#xff1a;存在a[i]a[i1]的情况&#xff0c;且在不影响利润的情况下&#xff0c;机器要尽可能的用完 这河里吗 可能是机器多了…

P7599-[APIO2021]雨林跳跃【二分,倍增,ST表】

正题 题目链接:https://www.luogu.com.cn/problem/P7599 题目大意 nnn棵树&#xff0c;在某棵树上时可以选择向左右两边第一棵比它高的树跳&#xff0c;现在qqq次询问从[A,B][A,B][A,B]中某个点出发跳到[C,D][C,D][C,D]中某个点的最少次数。 1≤n≤21051\leq n\leq 2\times 10…

2021中国大学生程序设计竞赛(CCPC)- 网络选拔赛 GCD on Sequence(线段树)

https://blog.csdn.net/qq_45863710/article/details/120121607 对于每个左端点 l 维护尚未确定 v 值的最大的右端点 mx[l] #include<bits/stdc.h>using namespace std; using lllong long; const int N100010; vector<int> vec[N]; int n,a[N],p[N]; ll ans[N]; …

牛客题霸 [ 验证IP地址] C++题解/答案

牛客题霸 [ 验证IP地址] C题解/答案 题目描述 编写一个函数来验证输入的字符串是否是有效的 IPv4 或 IPv6 地址 IPv4 地址由十进制数和点来表示&#xff0c;每个地址包含4个十进制数&#xff0c;其范围为 0 - 255&#xff0c; 用(".")分割。比如&#xff0c;172.1…

不止代码:洛谷P1064 金明的预算方案+P2014选课(依赖背包)

文章目录题目描述总结解析解法1解法2代码解法3代码题目描述 金明的预算方案 选课 金明今天很开心&#xff0c;家里购置的新房就要领钥匙了&#xff0c;新房里有一间金明自己专用的很宽敞的房间。更让他高兴的是&#xff0c;妈妈昨天对他说&#xff1a;“你的房间需要购买哪些物…

Docker最全教程——从理论到实战(六)

本篇教程持续编写了3个星期左右并且一直在完善、补充具体的细节和实践&#xff0c;预计全部完成需要1到2个月的时间。由于编写的过程中极其费时&#xff0c;并且还需要配合做一些实践&#xff08;有些实践存在一些坑&#xff0c;而且极其费时费事&#xff09;。因此目前产出的速…

P4332-[SHOI2014]三叉神经树【LCT】

正题 题目链接:https://www.luogu.com.cn/problem/P4332 题目大意 给出nnn个点的一棵有根三叉树&#xff0c;保证每个点的儿子个数为333或者000&#xff0c;每个叶子有一个权值000或111&#xff0c;每个非叶子节点的权值是它儿子中权值较多的那个&#xff0c;每次修改一个叶子…

牛客题霸 [ 环形链表的约瑟夫问题] C++题解/答案

牛客题霸 [ 环形链表的约瑟夫问题] C题解/答案 题目描述 据说著名犹太历史学家 Josephus 有过以下故事&#xff1a;在罗马人占领乔塔帕特后&#xff0c;39 个犹太人与 Josephus 及他的朋友躲到一个洞中&#xff0c;39 个犹太人决定宁愿死也不要被敌人抓到&#xff0c;于是决定…

2016 Multi-University Training Contest 1 1004 GCD(ST表+二分)

GCD 注意观察gcd⁡(al,al1,...,ar)\gcd(a_{l},a_{l1},...,a_{r})gcd(al​,al1​,...,ar​)&#xff0c;当lll固定不动的时候&#xff0c;rl...nrl...nrl...n时&#xff0c;我们可以容易的发现,随着rrr的増大&#xff0c;gcd⁡(al,al1,...,ar)\gcd(a_{l},a_{l1},...,a_{r})gcd(…

警卫站岗(树上dp)

题目描述 五一来临&#xff0c;某地下超市为了便于疏通和指挥密集的人员和车辆&#xff0c;以免造成超市内的混乱和拥挤&#xff0c;准备临时从外单位调用部分保安来维持交通秩序。 已知整个地下超市的所有通道呈一棵树的形状&#xff1b;某些通道之间可以互相望见。总经理要求…

.NET Core实战项目之CMS 第八章 设计篇-内容管理极简设计全过程

写在前面上一篇文章.NET Core实战项目之CMS 第七章 设计篇-用户权限极简设计全过程中我带着大家进行了权限部分的极简设计&#xff0c;也仅仅是一个基本的权限设计。不过你完全可以基于这套权限系统设计你的更复杂的权限系统&#xff0c;当然更复杂的权限系统要根据你的业务来进…

P7405-[JOI 2021 Final]雪玉【二分】

正题 题目链接:https://www.luogu.com.cn/problem/P7405 题目大意 nnn个点在坐标轴上&#xff0c;qqq次每次所有点向一个方向移动若干步&#xff0c;每个点的权值是它第一次覆盖的区间长度&#xff08;也就是一个区间只能贡献到第一次经过它的点&#xff09;。 求所有点的最终…

牛客题霸 [ 集合的所有子集] C++题解/答案

牛客题霸 [ 集合的所有子集] C题解/答案 题目描述 现在有一个没有重复元素的整数集合S&#xff0c;求S的所有子集 注意&#xff1a; 你给出的子集中的元素必须按升序排列 给出的解集中不能出现重复的元素 题解&#xff1a; 先sort一下&#xff0c;然后通过回溯实现选与不选…

2021“MINIEYE杯”中国大学生算法设计超级联赛(2)I love max and multiply(转化)

I love max and multiply Code 代码抄的std #include<bits/stdc.h> using namespace std; using lllong long; template <class Tint> T rd() {T res0;T fg1;char chgetchar();while(!isdigit(ch)) {if(ch-) fg-1;chgetchar();}while( isdigit(ch)) res(res<&…

人工智能第二课:认知服务和机器人框架探秘

这是《人工智能系列笔记》的第二篇&#xff0c;我利用周六下午完成课程学习。这一方面是因为内容属于入门级&#xff0c;并且之前我已经对认知服务和机器人框架比较熟悉。如有兴趣&#xff0c;请关注该系列 https://aka.ms/learningAI 但是学习这门课程还是很有收获&#xff0c…

洛谷P4322 最佳团伙(树上dp)

题目描述 洛谷传送门 题目描述 JSOI 信息学代表队一共有 N 名候选人&#xff0c;这些候选人从 1 到 N 编号。方便起见&#xff0c;JYY 的编号是 0 号。每个候选人都由一位编号比他小的候选人Ri推荐。如果 Ri0&#xff0c;则说明这个候选人是 JYY 自己看上的。 为了保证团队的…