这道题。。。神题。

首先看到数据范围,一眼状压 dp 。然后? 没了。

理性分析,这里说断掉任意一条边图依然连通,即整个图构成一个边双(而不是点双)。

之前用 fire (机房里的随机算法总称)之所以过不掉大概是我想错性质了(边双的条件直接变成 任意点的度大于2),

再打一遍 fire 不知道能不能过(可能关于判边双会十分麻烦,且复杂度较高)。

在这里我们考虑一个边双加上一条链(链的两端都在边双的集合里面),然后这个图依旧是边双。

那么在这里,一个点既是边双又是链,只不过要特判点 这一特殊情况。

然后我们可以处理出每一条链的最短长度,再进行边双最优解转移。

 //by Judge
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int inf=1e9+;
const int M=(<<)+;
#ifndef Judge
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
#endif
inline void cmin(int& a,int b){if(a>b)a=b;}
char buf[<<],*p1=buf,*p2=buf;
inline int read(){ int x=,f=; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
for(;isdigit(c);c=getchar()) x=x*+c-''; return x*f;
} int n,m,pat,msk,head[],log[M],h1[][M],h2[][M],g[][][M],f[M];
struct edge{ int to,val,next; }e[];
inline void add(int x,int y,int z){
e[++pat]={y,z,head[x]},head[x]=pat;
e[++pat]={x,z,head[y]},head[y]=pat;
}
#define lowbit(x) (x&-x)
int main(){ int T=read();
while(T--){ n=read(),m=read(),
msk=(<<n)-,pat=;
memset(head,,sizeof(head));
for(int i=;i<n;++i) log[<<i]=i;
for(int i=,a,b,c;i<m;++i)
a=read()-,b=read()-,
c=read(),add(a,b,c);
memset(h1,0x0f,sizeof(h1)),
memset(h2,0x0f,sizeof(h2)),
memset(g,0x0f,sizeof(g));
memset(f,0x0f,sizeof(f));
#define v e[i].to
for(int i,S=,a;S<=msk;++S){
for(int ta=msk^S;ta;ta-=lowbit(ta)){
for(a=log[lowbit(ta)],i=head[a];i;i=e[i].next) if((S>>v)&){
if(e[i].val<h1[a][S]) h2[a][S]=h1[a][S],h1[a][S]=e[i].val;
else cmin(h2[a][S],e[i].val);
}
}
}
for(int i=;i<n;++i) g[i][i][<<i]=f[<<i]=;
for(int S=;S<=msk;++S) for(int ta=S,a;ta;ta-=lowbit(ta))
for(int tb=S,b;tb;tb-=lowbit(tb)){
a=log[lowbit(ta)],b=log[lowbit(tb)]; if(a==b) continue;
for(int i=head[b];i;i=e[i].next) if((S>>v)&)
cmin(g[a][b][S],g[a][v][S^(<<b)]+e[i].val);
}
for(int S=,a,b;S<=msk;++S) if(S^lowbit(S)){
for(int T=(S-)&S;T;T=(T-)&S) for(int ta=T;ta;ta-=lowbit(ta))
for(int tb=T;tb;tb-=lowbit(tb)){ a=log[lowbit(ta)],b=log[lowbit(tb)];
if(a==b) cmin(f[S],f[S^T]+g[a][b][T]+h1[a][S^T]+h2[a][S^T]);
else cmin(f[S],f[S^T]+g[a][b][T]+h1[a][S^T]+h1[b][S^T]);
}
}
if(f[msk]>=f[]) puts("impossible");
else printf("%d\n",f[msk]);
} return ;
}

【BZOJ3590】[Snoi2013]Quare 状压DP的更多相关文章

  1. BZOJ3590 SNOI2013Quare(状压dp)

    可能作为最优解的边双都可以这样生成:初始时边双内只有一个点,每次选取边双内部两点(可以相同)和一个当前不在边双内的点集,以该两点为起止点找一条链(当然如果两点相同就是个环)将点集串起来,加入边双.状压 ...

  2. BZOJ 1087: [SCOI2005]互不侵犯King [状压DP]

    1087: [SCOI2005]互不侵犯King Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3336  Solved: 1936[Submit][ ...

  3. nefu1109 游戏争霸赛(状压dp)

    题目链接:http://acm.nefu.edu.cn/JudgeOnline/problemShow.php?problem_id=1109 //我们校赛的一个题,状压dp,还在的人用1表示,被淘汰 ...

  4. poj3311 TSP经典状压dp(Traveling Saleman Problem)

    题目链接:http://poj.org/problem?id=3311 题意:一个人到一些地方送披萨,要求找到一条路径能够遍历每一个城市后返回出发点,并且路径距离最短.最后输出最短距离即可.注意:每一 ...

  5. [NOIP2016]愤怒的小鸟 D2 T3 状压DP

    [NOIP2016]愤怒的小鸟 D2 T3 Description Kiana最近沉迷于一款神奇的游戏无法自拔. 简单来说,这款游戏是在一个平面上进行的. 有一架弹弓位于(0,0)处,每次Kiana可 ...

  6. 【BZOJ2073】[POI2004]PRZ 状压DP

    [BZOJ2073][POI2004]PRZ Description 一只队伍在爬山时碰到了雪崩,他们在逃跑时遇到了一座桥,他们要尽快的过桥. 桥已经很旧了, 所以它不能承受太重的东西. 任何时候队伍 ...

  7. bzoj3380: [Usaco2004 Open]Cave Cows 1 洞穴里的牛之一(spfa+状压DP)

    数据最多14个有宝藏的地方,所以可以想到用状压dp 可以先预处理出每个i到j的路径中最小权值的最大值dis[i][j] 本来想用Floyd写,无奈太弱调不出来..后来改用spfa 然后进行dp,这基本 ...

  8. HDU 1074 Doing Homework (状压dp)

    题意:给你N(<=15)个作业,每个作业有最晚提交时间与需要做的时间,每次只能做一个作业,每个作业超出最晚提交时间一天扣一分 求出扣的最小分数,并输出做作业的顺序.如果有多个最小分数一样的话,则 ...

  9. 【BZOJ1688】[Usaco2005 Open]Disease Manangement 疾病管理 状压DP

    [BZOJ1688][Usaco2005 Open]Disease Manangement 疾病管理 Description Alas! A set of D (1 <= D <= 15) ...

随机推荐

  1. layui(四)——table组件常见用法总结

    table是 layui 最核心的组成之一.它用于对表格进行一些列功能和动态化数据操作,涵盖了日常业务所涉及的几乎全部需求.支持固定表头.固定行.固定列左/列右,支持拖拽改变列宽度,支持排序,支持多级 ...

  2. TPS和QPS的区别和理解【转】

    QPS:Queries Per Second意思是“每秒查询率”,是一台服务器每秒能够相应的查询次数,是对一个特定的查询服务器在规定时间内所处理流量多少的衡量标准. TPS:是Transactions ...

  3. Golang入门教程(六)关键字和数据类型

    在 Go 编程语言中,数据类型用于声明函数和变量. 数据类型的出现是为了把数据分成所需内存大小不同的数据,编程的时候需要用大数据的时候才需要申请大内存,就可以充分利用内存. 一.25个关键字 二.18 ...

  4. AJAX工作原理与缺点

    1.概念:什么是AJAXAJAX全称为“Asynchronous JavaScript and XML”(异步JavaScript和XML),是一种创建交互式网页应用的网页开发技术.2.为什么要使用他 ...

  5. python中的顺序表

    Python中的list和tuple两种类型采用了顺序表的实现技术,tuple是不可变类型,即不变的顺序表,因此不支持改变其内部状态的任何操作,而其他方面,则与list的性质类似. list的基本实现 ...

  6. Android MediaPlayer播放raw资源封装类

    import android.content.Context; import android.media.MediaPlayer; import xxxx.R; public class MediaU ...

  7. Tooltip导致的无法访问已释放对象

    最近C#项目中遇到了一个无法访问已释放对象问题,经过反复测试,最终发现问题出在控件Tootip上,因为tootip内部有一个定时器,如果在窗口销毁时,鼠标移动到控件上恰好产生了一个tooltip,就会 ...

  8. IPv6地址测试宏

  9. 堆栈Stack介绍

    堆栈(Stack)代表了一个后进先出的对象集合.当您需要对各项进行后进先出的访问时,则使用堆栈.当您在列表中添加一项,称为推入元素,当您从列表中移除一项时,称为弹出元素. Stack 类的方法和属性 ...

  10. datatable转换为list<model> 映射

    using System; using System.Collections.Generic; using System.Data; using System.Linq; using System.R ...