[Codeforces1155F][DP]Delivery Oligopoly

翻译

给你一张双联通图，问你最少留下多少条边使得他还是双联通
这里的双联通指两点能够通过至少两条边不相交路径互达

题解

貌似是一个经典问题
fuck考场掉分一点都不想写这场题解来着…
我们考虑成为答案的双联通图是可以怎么构造的
我们可以每次拿出一条链，然后使得答案的双连通图变成一个更小的双联通图
然后一个思路就出来了
预处理一个g[S][x][y]g[S][x][y]g[S][x][y]表示从xxx到yyy的链，经过了状态SSS中的点，能否走出这样的链
然后考虑f[S]f[S]f[S]表示状态SSS中的点作为双联通分量的最小边数
转移直接枚举子集转移
记录一下方案就好了…

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#include<set>
#define LL long long
#define mp(x,y) make_pair(x,y)
#define pll pair<long long,long long>
#define pii pair<int,int>
using namespace std;
inline int read()
{int f=1,x=0;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;
}
int stack[20];
inline void write(int x)
{if(x<0){putchar('-');x=-x;}if(!x){putchar('0');return;}int top=0;while(x)stack[++top]=x%10,x/=10;while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);putchar(' ');}
inline void pr2(int x){write(x);putchar('\n');}
const int MAXN=15;
const int MAXMASK=(1<<14);
struct edge{int x,y,next;}a[MAXN*MAXN*2];int len,last[MAXN];
void ins(int x,int y){len++;a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;}
int g[MAXMASK][MAXN][MAXN],f[MAXMASK];
int pre1[MAXMASK][3],pre2[MAXMASK][MAXN][MAXN][3];
int n,m,nxt[MAXN],bin[25];
void upd(int &x,int y){x=x<y?x:y;}
void solve2(int S,int x,int y)
{if(x==y)return ;int o=pre2[S][x][y][0];int nS=S,nx=x,ny=y;if(!o)pr1(x),pr2(pre2[S][x][y][1]),nS^=bin[x],nx=pre2[S][x][y][1];else pr1(y),pr2(pre2[S][x][y][1]),nS^=bin[y],ny=pre2[S][x][y][1];solve2(nS,nx,ny);
}
void solve1(int S)
{if(!pre1[S][0])return ;int e1=pre1[S][0],e2=pre1[S][1],msk=pre1[S][2];solve2(msk,e1,e2);int ne=S^msk;for(int i=1;i<=n;i++)if(bin[i]&(ne&nxt[e1])){pr1(i);break;}pr2(e1);for(int i=n;i>=1;i--)if(bin[i]&(ne&nxt[e2])){pr1(i);break;}pr2(e2);solve1(ne);
}
int ct[MAXMASK];
int main()
{bin[1]=1;for(int i=2;i<=20;i++)bin[i]=bin[i-1]<<1;for(int i=1;i<MAXMASK;i++)ct[i]=(ct[i>>1]+(i&1));n=read();m=read();for(int i=1;i<=m;i++){int x=read(),y=read();ins(x,y),ins(y,x);nxt[x]|=bin[y];nxt[y]|=bin[x];}memset(g,63,sizeof(g));int INF=g[0][0][0];for(int i=1;i<=n;i++)g[bin[i]][i][i]=0;for(int S=1;S<bin[n+1];S++)for(int x=1;x<=n;x++)if(S&bin[x])for(int y=1;y<=n;y++)if(S&bin[y]&&g[S][x][y]!=INF){for(int k=last[x];k;k=a[k].next)if(!(S&bin[a[k].y])){upd(g[S|bin[a[k].y]][a[k].y][y],g[S][x][y]+1);if(g[S|bin[a[k].y]][a[k].y][y]==g[S][x][y]+1)pre2[S|bin[a[k].y]][a[k].y][y][0]=0,pre2[S|bin[a[k].y]][a[k].y][y][1]=x;}for(int k=last[y];k;k=a[k].next)if(!(S&bin[a[k].y])){upd(g[S|bin[a[k].y]][x][a[k].y],g[S][x][y]+1);if(g[S|bin[a[k].y]][x][a[k].y]==g[S][x][y]+1)pre2[S|bin[a[k].y]][x][a[k].y][0]=1,pre2[S|bin[a[k].y]][x][a[k].y][1]=y;}}memset(f,63,sizeof(f));for(int i=1;i<=n;i++)f[bin[i]]=0;for(int S=1;S<bin[n+1];S++){for(int u=(S-1)&S;u;u=(u-1)&S)if(f[S^u]!=INF){int ps=S^u;for(int x=1;x<=n;x++)if((u&bin[x])&&(nxt[x]&ps))for(int y=x+1;y<=n;y++)if((u&bin[y])&&(nxt[y]&ps)&&g[u][x][y]!=INF){ upd(f[S],f[ps]+g[u][x][y]+2);if(f[S]==f[ps]+g[u][x][y]+2)pre1[S][0]=x,pre1[S][1]=y,pre1[S][2]=u;}if(ct[u]==1){for(int x=1;x<=n;x++)if(u&bin[x]&&ct[nxt[x]&ps]>=2){upd(f[S],f[ps]+2);if(f[S]==f[ps]+2)pre1[S][0]=x,pre1[S][1]=x,pre1[S][2]=u;}}}}pr2(f[bin[n+1]-1]);solve1(bin[n+1]-1);return 0;
}

[Codeforces1155F][DP]Delivery Oligopoly相关推荐

codeforces1155F Delivery Oligopoly
题面题意给出一张无向图,问至少保留几条边才能使此图是一个边双,并输出方案. 做法考虑怎么构造出一个边双,发现可以通过在一个小边双上加一条链,使它变为一个新的大边双,所以可以记f[u][v][i] ...
Codeforces 1155F Delivery Oligopoly dp（看题解）
看别人写的才学会的... 我们考虑刚开始的一个点, 然后我们枚举接上去的一条一条链, dp[mask]表示当前已经加进去点的状态是mask所需的最少边数. 反正就是很麻烦的一道题, 让我自己写我是写不 ...
Codeforces Round 63 (Rated for Div. 2) F. Delivery Oligopoly dp+图论状态转移
题目链接: https://codeforces.com/contest/1155/problem/F 题意: 现在给你一个 141414 个点的无向边双联通图,现在要你删掉一些边,使得留下来的边最少 ...
CF1155F Delivery Oligopoly
题意:给定简单无向图,求一个最小的边集使得这些点是边双,输出方案.n <= 14 解:考虑一个边双肯定是一条一条的链拼起来的.于是每次枚举一条链加上去就行了. 设fs表示点集s形成边双的最小边数 ...
【CodeForces】【状压DP】1155F-Delivery Oligopoly
CodeForces 1155F Delivery Oligopoly 题目大意给定一个已经是边双联通分量的图,要求删掉最多的边,使得最终得到的图是也是一个边双联通分量.输出保留的边. 分析考虑我 ...
java秒转换为年月日_SimpleDateFormat将月/日/年时分秒转换为年-月-日时:分:秒
String expirTime = "12 / 27 / 2018 12: 00: 00 AM"; SimpleDateFormat in = new SimpleDateFor ...
NCPC2018 D.Delivery Delays[二分答案+DP check]
Delivery Delays 题意 100010001000个点,500050005000条边的无向图,披萨店在111号店.100010001000份披萨订单,每个订单有下单时间,送达地点,披萨制作 ...
ZOJ3469 Food Delivery 区间DP
题意:有一家快餐店送外卖,现在同时有n个家庭打进电话订购,送货员得以V-1的速度一家一家的运送,但是每一个家庭都有一个不开心的值,每分钟都会增加一倍,值达到一定程度,该家庭将不会再订购外卖了,现在为了 ...
ZOJ ~ 3469 ~ Food Delivery （区间DP）
题意一个人要送外卖,现在,在X轴上有N个人,先输入n,x,v,表示n个人,骑手一开始在x位置,他的速度是v,然后N行输入X[i]和B[i]表示这个人在x[i]位置,每分钟的等待都会增加b[i]的愤怒 ...

[Codeforces1155F][DP]Delivery Oligopoly

[Codeforces1155F][DP]Delivery Oligopoly相关推荐

最新文章

热门文章