【BZOJ4774】修路

Description

村子间的小路年久失修,为了保障村子之间的往来,法珞决定带领大家修路。对于边带权的无向图 G = (V, E),请选择一些边,使得1 <= i <= d, i号节点和 n - i + 1 号节点可以通过选中的边连通,最小化选中的所有边的权值和。

Input

第一行两个整数 n, m,表示图的点数和边数。接下来的 m行,每行三个整数 ui, vi, wi,表示有一条 ui 与 vi 之间,权值为 wi 的无向边。
1 <= d <= 4
2d <= n <= 10^4
0 <= m <= 10^4
1 <= ui, vi <= n
1 <= wi <= 1000

Output

一行一个整数,表示答案,如果无解输出-1

Sample Input

10 20 1
6 5 1
6 9 4
9 4 2
9 4 10
6 1 2
2 3 6
7 6 10
5 7 1
9 7 2
5 9 10
1 6 8
4 7 4
5 7 1
2 6 9
10 10 6
8 7 2
10 9 10
1 2 4
10 1 8
9 9 7

Sample Output

8

题解:设f[S][i]表示已经连通的关键点状态为S,当前位于点i的最小权值和。转移就是斯坦纳树。

再设g[S]表示已经连通的关键点状态为S的最小权值和。只有当S保证所有的关键点对的连通状态相同时,才可以从f[S][..]更新到g[S],最后对g数组跑枚举子集的DP即可。

P.S.我至今才会枚举子集的正确姿势~~~

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <utility>
#define mp(A,B) make_pair(A,B)
using namespace std;
int n,m,d,cnt,now,tot;
int f[1<<8][10010],to[20010],next[20010],val[20010],head[10010],vis[10010],p[20],Log[1<<8],ref[1<<8],g[1<<8];
priority_queue<pair<int,int> > q;
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-') f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+(gc^'0'),gc=getchar();
return ret*f;
}
inline void add(int a,int b,int c)
{
to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
}
inline void dij(int S,int x)
{
q.push(mp(-f[S][x],x));
int i,u;
now++;
while(!q.empty())
{
u=q.top().second,q.pop();
if(vis[u]==now) continue;
vis[u]=now;
for(i=head[u];i!=-1;i=next[i]) if(f[S][to[i]]>f[S][u]+val[i])
f[S][to[i]]=f[S][u]+val[i],q.push(mp(-f[S][to[i]],to[i]));
}
}
int main()
{
int i,a,b,c,x,y;
n=rd(),m=rd(),d=rd();
memset(head,-1,sizeof(head));
for(i=1;i<=m;i++) a=rd(),b=rd(),c=rd(),add(a,b,c),add(b,a,c);
for(i=0;i<2*d;i++) Log[1<<i]=i;
memset(f,0x3f,sizeof(f)),memset(g,0x3f,sizeof(g));
for(i=1;i<=d;i++)
{
f[1<<(i-1)][i]=0,dij(1<<(i-1),i);
f[1<<(d+i-1)][n-i+1]=0,dij(1<<(d+i-1),n-i+1);
}
for(x=1;x<(1<<(2*d));x++)
{
for(tot=0,y=x;y;y-=y&-y) p[tot++]=y&-y;
for(y=1;y<(1<<tot);y++)
{
ref[y]=ref[y^(y&-y)]|p[Log[y&-y]];
for(i=1;i<=n;i++) if(f[x][i]>f[ref[y]][i]+f[x^ref[y]][i])
f[x][i]=f[ref[y]][i]+f[x^ref[y]][i];
}
for(i=1;i<=n;i++) dij(x,i);
}
for(x=1;x<(1<<d);x++) for(i=1;i<=n;i++) g[x]=min(g[x],f[(x<<d)|x][i]);
for(x=1;x<(1<<d);x++)
{
for(tot=0,y=x;y;y-=y&-y) p[tot++]=y&-y;
for(y=1;y<(1<<tot);y++)
{
ref[y]=ref[y^(y&-y)]|p[Log[y&-y]];
g[x]=min(g[x],g[ref[y]]+g[x^ref[y]]);
}
}
if(g[(1<<d)-1]==0x3f3f3f3f) printf("-1");
else printf("%d",g[(1<<d)-1]);
return 0;
}
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <utility>
#define mp(A,B) make_pair(A,B)
using namespace std;
int n,m,d,cnt,now;
int to[6010],next[6010],val[6010],head[1010],f[1<<10][1010],g[1<<10],p[15],w[15],vis[1<<10];
priority_queue<pair<int,int> > q;
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-') f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+(gc^'0'),gc=getchar();
return ret*f;
}
inline void add(int a,int b,int c)
{
to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
}
inline void dij(int S,int x)
{
q.push(mp(-f[S][x],x));
int i,u;
now++;
while(!q.empty())
{
u=q.top().second,q.pop();
if(vis[u]==now) continue;
vis[u]=now;
for(i=head[u];i!=-1;i=next[i]) if(f[S][to[i]]>f[S][u]+val[i])
f[S][to[i]]=f[S][u]+val[i],q.push(mp(-f[S][to[i]],to[i]));
}
}
int main()
{
n=rd(),m=rd(),d=rd();
int i,j,a,b,c,x,y;
memset(head,-1,sizeof(head));
for(i=1;i<=m;i++) a=rd(),b=rd(),c=rd(),add(a,b,c),add(b,a,c);
memset(f,0x3f,sizeof(f)),memset(g,0x3f,sizeof(g));
for(i=0;i<d;i++) w[i]=rd(),p[i]=rd(),f[1<<i][p[i]]=0,dij(1<<i,p[i]);
for(x=1;x<(1<<d);x++)
{
for(y=(x-1)&x;y;y=(y-1)&x) for(i=1;i<=n;i++) f[x][i]=min(f[x][i],f[y][i]+f[x^y][i]);
for(i=1;i<=n;i++) dij(x,i);
}
for(x=1;x<(1<<d);x++) for(i=1;i<=n;i++) g[x]=min(g[x],f[x][i]);
for(x=1;x<(1<<d);x++)
{
for(i=0;i<d;i++) for(j=0;j<i;j++) if(w[i]==w[j]&&((x>>i)&1)!=((x>>j)&1)) g[x]=0x3f3f3f3f;
}
for(x=1;x<(1<<d);x++) for(y=(x-1)&x;y;y=(y-1)&x) g[x]=min(g[x],g[y]+g[x^y]);
printf("%d",g[(1<<d)-1]);
return 0;
}

【BZOJ4774/4006】修路/[JLOI2015]管道连接 斯坦纳树的更多相关文章

  1. BZOJ 4006 Luogu P3264 [JLOI2015]管道连接 (斯坦纳树、状压DP)

    题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4006 (luogu)https://www.luogu.org/probl ...

  2. BZOJ4006 JLOI2015 管道连接(斯坦纳树生成森林)

    4006: [JLOI2015]管道连接 Time Limit: 30 Sec Memory Limit: 128 MB Description 小铭铭最近进入了某情报部门,该部门正在被如何建立安全的 ...

  3. BZOJ4006: [JLOI2015]管道连接(斯坦纳树,状压DP)

    Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 1171  Solved: 639[Submit][Status][Discuss] Descripti ...

  4. 【bzoj4006】[JLOI2015]管道连接 斯坦纳树+状压dp

    题目描述 给出一张 $n$ 个点 $m$ 条边的无向图和 $p$ 个特殊点,每个特殊点有一个颜色.要求选出若干条边,使得颜色相同的特殊点在同一个连通块内.输出最小边权和. 输入 第一行包含三个整数 n ...

  5. 洛谷P3264 [JLOI2015]管道连接 (斯坦纳树)

    题目链接 题目大意:有一张无向图,每条边有一定的花费,给出一些点集,让你从中选出一些边,用最小的花费将每个点集内的点相互连通,可以使用点集之外的点(如果需要的话). 算是斯坦纳树的入门题吧. 什么是斯 ...

  6. bzoj 4006 [JLOI2015]管道连接——斯坦纳树

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4006 除了模板,就是记录 ans[ s ] 表示 s 合法的最小代价.合法即保证 s 里同一 ...

  7. 【BZOJ4774】修路(动态规划,斯坦纳树)

    [BZOJ4774]修路(动态规划,斯坦纳树) 题面 BZOJ 题解 先讲怎么求解最小斯坦纳树. 先明白什么是斯坦纳树. 斯坦纳树可以认为是最小生成树的一般情况.最小生成树是把所有给定点都要加入到联通 ...

  8. bzoj 4006 管道连接 —— 斯坦纳树+状压DP

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4006 用斯坦纳树求出所有关键点的各种连通情况的代价,把这个作为状压(压的是集合选择情况)的初 ...

  9. 【LuoguP3264】[JLOI2015] 管道连接(斯坦那树)

    题目链接 题目描述 小铭铭最近进入了某情报部门,该部门正在被如何建立安全的通道连接困扰.该部门有 n 个情报站,用 1 到 n 的整数编号.给出 m 对情报站 ui;vi 和费用 wi,表示情报站 u ...

随机推荐

  1. python三大神器之virtualenv pip, virtualenv, fabric通称为pythoner的三大神器。

    python三大神器之virtualenv   pip, virtualenv, fabric通称为pythoner的三大神器. virtualenv virtualenv------用来建立一个虚拟 ...

  2. 错误将UIViewController当做UITableViewController来用

  3. nodejs 模块以及加载机制,主要讨论找不到模块的问题

    最主要的一个思想,加载模块无非就是找到模块在哪,只要清楚了模块的位置以及模块加载的逻辑那么找不到模块的问题就迎刃而解了.本文只是综合了自己所学的知识点进行总结,难免出现理解错误的地方,请见谅. nod ...

  4. C/C++-中的sort排序用法

    转载于:http://www.cnblogs.com/luorende/p/6121906.htmlSTL中就自带了排序函数sortsort 对给定区间所有元素进行排序 要使用此函数只需用#inclu ...

  5. CSRF简单介绍及利用方法

    x00 简要介绍 CSRF(Cross-site request forgery)跨站请求伪造,由于目标站无token/referer限制,导致攻击者可以用户的身份完成操作达到各种目的.根据HTTP请 ...

  6. 关于在Android中访问和使用到上下文变量

    在监听器内部实现类中要引用上下文变量this的时候 一.采用类名.this的方法 FActivity.this 二.采用全局变量当做中间变量 1.先定义一个全局变量 private Context m ...

  7. 关于在Android中Activity页面跳转的方法

    一.无返回结果的页面跳转 1.创建两个类FActivity.java和SActivity.java 2.创建两个layout目录下的factivity.xml和sactivity.xml 3.在And ...

  8. (转)从海康7816的ps流里获取数据h264数据

    海康7816使用ps流来封装h.264数据,这里使用的解码器无法识别ps流,因此需要将h264数据从ps流里提取出来 对于ps流的规定可以参考13818-1文档 这里从7816里获取到一些数据取样 0 ...

  9. Android开发之经常使用的时间格式

    /**   * 获取如今时间   *    * @return 返回时间类型 yyyy-MM-dd HH:mm:ss   */ public static Date getNowDate() {   ...

  10. C语言中FILE是结构体,文件类型的指针

    c语言文件类型指针 我们在定义文件类型指针变量后,称作该指针指向该文件,但本质上,它不是指向一个存储文件信息的结构型变量么?那么我们在用各个函数对所谓的“文件指针”进行操作时,本质上是不是函数通过获取 ...