[JZOJ3348] 【NOI2013模拟】秘密任务
题目
题目大意
给你一个无向图,你要割掉一些边使得\(1\)到\(n\)的所有最短路径被阻截。
割掉一个边\((u,v)\)的代价为\(a_u\)或\(a_v\)(记为两种不同的方案)。
问最小代价及其唯一性。
思考历程
首先要将最短路图给建出来。
然后我就莫名其妙地想到了支配树,还在这个方向上思考了一阵子……
想不到怎么做……
然后我又想到之前某道题,将最大反链长度转化为最小链覆盖,然后用二分图匹配来实现。
我模仿着打了个KM算法,然后发现果然错了……
看来是不能直接生搬硬套的……
然后蓦然发现:这道题不就是个最小割吗!
于是就开始狂打……
最终的问题就是:如何判断唯一性……
匆匆忙忙地打了个错误的方法,交了上去。
正解
这题当然是最小割啦……
由于每条边选\(a_u\)和\(a_v\)算作两种方案,所以每条边要拆成两条来搞。
建图当然显然了。
至于判断是否有唯一性,可以把最小割后的残量网络的\(S\)集合和\(T\)集合找出来。
枚举\(S\)集合和\(T\)集合之间被割掉的边,求它们的权值和。
跟最小割进行对比,如果相等,则具有唯一性。
代码
using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <vector>
#define N 410
#define M 4010
int n,m;
int a[N];
struct EDGE{
int to,len;
EDGE *las;
} e[M*2];
int ne;
EDGE *last[N*2];
inline void link(int u,int v,int len){
e[ne]={v,len,last[u]};
last[u]=e+ne++;
}
long long dis[N+M];
struct EDGE2{
int to,c;
EDGE2 *las;
} e2[M*4];
int ne2;
EDGE2 *last2[N+M];
#define rev(ei) (e2+(int(ei-e2)^1))
inline void link2(int u,int v,int c){
e2[ne2]={v,c,last2[u]};
last2[u]=e2+ne2++;
e2[ne2]={u,0,last2[v]};
last2[v]=e2+ne2++;
}
int S,T;
inline void S_P(){
static int q[100001];
static bool inq[N];
int head=0,tail=1;
q[1]=1;
memset(inq,0,sizeof inq);
inq[1]=1;
memset(dis,127,sizeof dis);
dis[1]=0;
do{
int x=q[++head];
for (EDGE *ei=last[x];ei;ei=ei->las)
if (dis[x]+ei->len<dis[ei->to]){
dis[ei->to]=dis[x]+ei->len;
if (!inq[ei->to]){
inq[ei->to]=1;
q[++tail]=ei->to;
}
}
inq[x]=0;
}
while (head!=tail);
memset(last2,0,sizeof last2);
ne2=0;
S=1,T=n;
for (int i=1;i<=T;++i){
for (EDGE *ei=last[i];ei;ei=ei->las)
if (dis[i]+ei->len==dis[ei->to]){
++n;
link2(i,n,a[i]);
link2(n,ei->to,a[ei->to]);
}
}
}
bool BZ;
int gap[N+M];
EDGE2 *cur[N+M];
int dfs(int x,int s){
if (x==T)
return s;
int have=s;
for (EDGE2 *ei=cur[x];ei;ei=ei->las){
cur[x]=ei;
if (ei->c && dis[ei->to]+1==dis[x]){
int t=dfs(ei->to,min(ei->c,have));
ei->c-=t,rev(ei)->c+=t,have-=t;
if (!have)
return s;
}
}
cur[x]=last2[x];
if (!--gap[dis[x]])
BZ=0;
dis[x]++;
gap[dis[x]]++;
return s-have;
}
inline long long flow(){
long long res=0;
memset(gap,0,sizeof gap);
memset(dis,0,sizeof dis),
gap[0]=n;
BZ=1;
while (BZ)
res+=dfs(S,INT_MAX);
return res;
}
inline bool pd(int mincut){
static int q[N+M];
static int vis[N+M];
int head=0,tail=1;
q[1]=S;
memset(vis,0,sizeof vis);
vis[S]=1;
do{
int x=q[++head];
for (EDGE2 *ei=last2[x];ei;ei=ei->las)
if (ei->c && !vis[ei->to]){
vis[ei->to]=1;
q[++tail]=ei->to;
}
}
while (head!=tail);
head=0,tail=1;
q[1]=T;
vis[T]=2;
do{
int x=q[++head];
for (EDGE2 *ei=last2[x];ei;ei=ei->las)
if (rev(ei)->c && !vis[ei->to]){
vis[ei->to]=2;
q[++tail]=ei->to;
}
}
while (head!=tail);
long long sum=0;
for (int i=1;i<=tail;++i)
for (EDGE2 *ei=last2[q[i]];ei;ei=ei->las)
if (rev(ei)->c==0 && vis[ei->to]==1)
sum+=ei->c;
return sum==mincut;
}
int main(){
freopen("in.txt","r",stdin);
// freopen("out.txt","w",stdout);
int T;
scanf("%d",&T);
while (T--){
scanf("%d%d",&n,&m);
for (int i=1;i<n;++i)
scanf("%d",&a[i]);
a[n]=INT_MAX;
memset(last,0,sizeof last);
ne=0;
for (int i=1;i<=m;++i){
int u,v,len;
scanf("%d%d%d",&u,&v,&len);
link(u,v,len),link(v,u,len);
}
S_P();
long long mincut=flow();
if (pd(mincut))
printf("Yes %d\n",mincut);
else
printf("No %d\n",mincut);
}
return 0;
}
总结
这题的最大收获就是唯一性的判定了……
[JZOJ3348] 【NOI2013模拟】秘密任务的更多相关文章
- 【NOI2013模拟】坑带的树(仙人球的同构+圆方树乱搞+计数+HASH)
[NOI2013模拟]坑带的树 题意: 求\(n\)个点,\(m\)条边的同构仙人球个数. \(n\le 1000\) 这是一道怎么看怎么不可做的题. 这种题,肯定是圆方树啦~ 好,那么首先转为广义圆 ...
- [NOI2013模拟] BZOJ4705 棋盘游戏 解题报告(组合计数)
莫名打不开这道题的链接,请读者自行搜索 Description 有一个N*M的棋盘,初始每个格子都是白色的.行操作是指选定某一行,将这行所有格子的颜色取反(黑白互换).列操作是指选定某一列,将这列所有 ...
- [JZOJ3362] 【NOI2013模拟】数数
题目 题目大意 求区间\([A,B]\)有多少个数是"完美的". 一个数是"完美的",当且仅当这个数的各位能分成两个集合,使得两个集合中数字的和相等. \(B\ ...
- [JZOJ3347] 【NOI2013模拟】树的难题
题目 题目大意 给你一棵树,每个节点有三种黑.白.灰三种颜色. 你要割掉一些边(每条边被割需要付出一定的代价),使得森林的每棵树满足: 没有黑点或至多一个白点. 思考历程 这题一看就知道是一个树形DP ...
- [JZOJ3339]【NOI2013模拟】wyl8899和法法塔的游戏
题目 题目大意 给你一个数列,每次给出\(r,a,b\),你要找到\(l\in [a,b]\)使得\([l,r-1]\)的异或和最小, 并且要修改\(r\)位置的数. 思考历程 当我看到这题的时候,已 ...
- [JZOJ3337] 【NOI2013模拟】wyl8899的TLE
题目 题目大意 给你两个字符串\(A\)和\(B\),可以修改\(A\)中的一个字符,求修改后最长的\(A\)的前缀,使它是\(B\)的子串. 思考历程 看到这道题之后,第一眼想到的就是后缀自动机! ...
- 忘记秘密利用python模拟登录暴力破解秘密
忘记秘密利用python模拟登录暴力破解秘密: #encoding=utf-8 import itertools import string import requests def gen_pwd_f ...
- JZOJ 4272. 【NOIP2015模拟10.28B组】序章-弗兰德的秘密
272. [NOIP2015模拟10.28B组]序章-弗兰德的秘密 (File IO): input:frand.in output:frand.out Time Limits: 1000 ms M ...
- [jzoj 5781]【NOIP提高A组模拟2018.8.8】秘密通道 (最短路)
传送门 Description 有一副nm的地图,有nm块地,每块是下列四种中的一种: 墙:用#表示,墙有4个面,分别是前面,后面,左面,右面. 起点:用C表示,为主角的起点,是一片空地. 终点:用F ...
随机推荐
- SVM-SVR
高频率的接触到了SVM模型,而且还有使用SVM模型做回归的情况,即SVR.另外考虑到自己从第一次知道这个模型到现在也差不多两年时间了,从最开始的腾云驾雾到现在有了一点直观的认识,花费了不少时间.因此在 ...
- 矢量切片应用中geoserver与geowebcache分布式部署方案
在进行GIS项目开发中,常使用Geoserver作为开源的地图服务器,Geoserver是一个JavaEE项目,常通过Tomcat进行部署.而GeoWebCache是一个采用Java实现用于缓存WMS ...
- H5新增input表单、表单属性
新增表单 email,Email类型 url , Url类型 date,日期类型 time,时间类型 month,月类型 week,周类型 number,数字类型 tel,电话类型 search,搜索 ...
- 2018-2-13-win10-UWP-动画
title author date CreateTime categories win10 UWP 动画 lindexi 2018-2-13 17:23:3 +0800 2018-2-13 17:23 ...
- Spring 在Web中的应用
Spring 在Web中的应用 在web项目开发中,不会直接实例化ApplicationContext对象,如果想用到ApplicationContext,一般的步骤: 配置一个监听器ContextL ...
- jdk源码阅读-ConcurrentLinkedQueue(一)
说明 concurrentLinkedQueue为无界非阻塞队列,是线程安全的 内部结构为链表的形式, 内部使用cas保存线程安全.采用cas保证原子性 什么是CAS CAS 操作包含三个操作数 —— ...
- Dubbo中有哪些角色?
registry 注册中心. 是用于发布和订阅服务的一个平台.用于替代SOA结构体系框架中的ESB服务总线的. 发布 开发服务端代码完毕后, 将服务信息发布出去. 实现一个服务的公开. 订阅 客户端程 ...
- CSIC_716_20191128【多态、绑定与非绑定方法、isinstance与issubclass 】
多态 what: 同一个事物有不同的形态. 多态的目的:在不知道对象具体类型的情况下,统一对象调用方法的规范.(通俗讲,即规定实现同样功能的方法,一定要起同样的名字). 多态的表现形式之一就是继承,先 ...
- vue computed和methods 计算属性和侦听器
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...
- IDEA 注解开发流程
IDEA 注解开发流程 分为以下四步 分别是 1 导入依赖 2 配置文件 applicationContext.xml 3 在需要创建对象的类上添加注解 4 测试 详细步骤和代码如下 1 导入依 ...