HDOJ题目页面传送门

给定一个无向带权图\(G=(V,E),|V|=n,|E|=m\),求边权之和最大的有\(s\)个节点的链的边权之和,即求\(\max\limits_{\forall i\in[1,s],\forall j\in(i,s],a_i\ne a_j,\forall i\in[1,s),(a_i,a_{i+1},l_i)\in E}\left\{\sum\limits_{i=1}^{s-1}l_i\right\}\)。如果没有符合要求的链,输出\(\texttt{impossible}\)。

\(n\in\left[2,10^4\right],m\in\left[1,10^4\right],s\in[2,6]\)。本题多测,最多有\(35\)组数据,\(\max(n,m)>100\)的最多有\(5\)组。时限\(15\mathrm s\)。

先说\(1\)种不是正解的玄学方法

暴搜。。。时间复杂度\(\mathrm O(\mathrm C_n^s)\),但很显然非常跑不满,因为这是个稀疏图。而且还可以最优性剪枝优化,假设所有边中最大的边权为\(longest\),当前的最大答案为\(ans\),当前还剩\(rst\)个节点(条边)要搜,目前的边权之和为\(now\),那么如果\(now+rst\times longest\le ans\),便可立即停止搜索。这样就可以水过去了。。。

这不是主要内容,时间复杂度到底是多少就不研究了,代码也不贴了。


正解

考虑对一个比原问题弱一点的问题求解:给每个节点\(i\)染一个\([0,s)\)的颜色\(col_i\),求边权之和最大的有\(s\)个节点且这\(s\)个节点的颜色是\(0\sim s-1\)的一个排列的链的边权之和。这个问题很好求,状压DP就可以了。设\(dp_{i,j}\)表示边权之和最大的以\(i\)为\(1\)个端点,上面节点颜色互不相同且bitmask为\(j\)的链的边权之和,边界为\(dp_{i,\{j\}}=\begin{cases}0&i=j\\-\infty&i\ne j\end{cases}\),状态转移方程为\(dp_{i,j}=\begin{cases}\max\limits_{(i,k,l)\in E}\left\{dp_{k,j-\{col_i\}}+l\right\}&col_i\in j\\-\infty&col_i\notin j\end{cases}\),最终答案为\(\max\limits_{i=1}^{n}\{dp_{i,[1,n]\cap\mathbb Z}\}\)。DP顺序要注意,要按照\(j\)的递增顺序DP,不能按照\(i\)的顺序(我就在上面WA过)。(简单DP

那么既然这个弱问题这么好求,那我们就强制给每个节点染色。假设我们给每个节点随机染色,那我们来分析一下弱问题的答案就是原问题的正确答案的概率:如果原问题所求得的最长链,在弱问题中,上面所有节点的颜色正好组成\(0\sim s-1\)的一个排列,那么弱问题的答案就是原问题的答案,而在所有\(s^s\)种最长链的颜色序列中,有\(s!\)种是\(0\sim s-1\)的排列,所以概率为\(\dfrac{s!}{s^s}\)。当\(s=6\)时,概率最小为\(\dfrac{6!}{6^6}\approx1.5\%\)。那么我们只需要在\(35\)组数据中每组随机染色、求解足够多次,就有足够大的的概率在每组中都至少有一次弱问题的答案是原问题的答案,这样把每次弱问题的答案\(\max\)起来即可。我们希望这个概率至少为\(99\%\),设随机\(x\)次,那么可以列出不等式\(\left(1-\left(1-\dfrac{6!}{6^6}\right)^x\right)^{35}\ge99\%\)。不难解出\(x\ge524.397\cdots\)。但由于HDOJ比较卡常,随机\(525\)次会TLE,于是我们只能用正确率换时间,只随机\(300\)次,这样才能勉强卡进\(15\mathrm s\)的时限。此时正确率为\(\left(1-\left(1-\dfrac{6!}{6^6}\right)^{300}\right)^{35}\approx71.8\%\),侥幸过。

btw,随机的时候不能写rand()*rand()%s。第一,在某些系统下两个rand()乘起来会爆int;第二,把两个随机数乘起来,期望因数个数会增加很多,\(\bmod s\)的结果等于\(0\)的概率也就会增加很多,就不均匀了。可以这样写:abs(rand()*rand()*rand()+rand())%s

UPD 2020.4.12:你看看,ycx,你看看,上面那段删除线写的是smjbwy,狗屁不通。。看来我还是太年轻了。昨天经hsc神仙的提醒才知道,abs(rand()*rand()*rand()+rand())这样的写法会生成出来分布极度不均匀的随机数,虽然在本题里没事,但显然不应该这么写。不如直接用mt19937,先定义一个随机数生成器mt19937 rng;,那么\([l,r]\)之间的随机数就是uniform_int_distribution<>(l,r)(rng)

下面贴正解的代码:UPDed 2020.4.12

#include<bits/stdc++.h>
using namespace std;
#define int long long//爆int
#define pb push_back
#define mp make_pair
#define X first
#define Y second
#define uid uniform_int_distribution
const int inf=0x3f3f3f3f3f3f3f3f;
mt19937 rng(20060617/*信仰优化*/);
int ppc(int x){return __builtin_popcount(x);}
const int N=10000,S=6;
int n/*节点数*/,m/*边数*/,s/*要求的链的节点数*/;
vector<pair<int,int> > nei[N+1];//邻接表
int col[N+1];//颜色
int ans;//原问题答案
int dp[N+1][1<<S];
void random(){
for(int i=1;i<=n;i++)col[i]=uid<>(0,s-1)(rng);//随机染色
for(int i=1;i<=n;i++)for(int j=1;j<1<<s;j++)dp[i][j]=j==1<<col[i]?0:-inf;//处理边界,顺便初始化
for(int j=1;j<1<<s;j++)for(int i=1;i<=n;i++)if(j&1<<col[i]&&ppc(j)>1/*大小为1的bitmask是边界*/)//枚举状态
for(int k=0;k<nei[i].size();k++){//转移
int x=nei[i][k].X,len=nei[i][k].Y;
dp[i][j]=max(dp[i][j],dp[x][j^1<<col[i]]+len);
}
for(int i=1;i<=n;i++)ans=max(ans,dp[i][(1<<s)-1]);//把每次弱问题的答案max起来就有大概率是原问题的答案
}
void mian(){
cin>>n>>m>>s;
for(int i=1;i<=n;i++)nei[i].clear();//数据不清空,爆零两行泪
while(m--){
int x,y,z;
scanf("%lld%lld%lld",&x,&y,&z);
nei[x].pb(mp(y,z));nei[y].pb(mp(x,z));
}
ans=-inf;//初始化
int randomnum=300;//随机次数
while(randomnum--)random();//随机
if(ans>=0)cout<<ans<<"\n";
else puts("impossible");
}
signed main(){
int testnum;//数据组数
cin>>testnum;
while(testnum--)mian();
return 0;
}

HDOJ 6664 Andy and Maze的更多相关文章

  1. 【HDOJ6664】Andy and Maze(color coding)

    题意:给定一张n点m边的无向带权图,问从任意结点出发,不能走已经经过的点,共经过k个点的最长路径的值 n,m<=1e4,k<=6 思路:color coding算法 考虑每次给每个点随机编 ...

  2. 2019DX#8

    Solved Pro.ID Title Ratio(Accepted / Submitted)   1001 Acesrc and Cube Hypernet 7.32%(3/41)   1002 A ...

  3. hdoj 5094 Maze 【BFS + 状态压缩】 【好多坑】

    Maze Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 100000/100000 K (Java/Others) Total Sub ...

  4. HDU 4067 hdoj 4067 Random Maze 最小费用流

    给出n个点,m条边,入口s和出口t,对于每条边有两个值a,b,如果保留这条边需要花费:否则,移除这条边需要花费b. 题目要求用最小费用构造一个有向图满足以下条件: 1.只有一个入口和出口 2.所有路都 ...

  5. 【HDOJ】1484 Basic wall maze

    BFS. /* 1484 */ #include <iostream> #include <queue> #include <string> #include &l ...

  6. BFS+贪心 HDOJ 5335 Walk Out

    题目传送门 /* 题意:求从(1, 1)走到(n, m)的二进制路径值最小 BFS+贪心:按照标程的作法,首先BFS搜索所有相邻0的位置,直到1出现.接下去从最靠近终点的1开始, 每一次走一步,不走回 ...

  7. 最大流增广路(KM算法) HDOJ 1533 Going Home

    题目传送门 /* 最小费用流:KM算法是求最大流,只要w = -w就可以了,很经典的方法 */ #include <cstdio> #include <cmath> #incl ...

  8. BFS HDOJ 1728 逃离迷宫

    题目传送门 /* BFS:三维BFS,加上方向.用dp[x][y][d]记录当前需要的最少转向数 */ #include <cstdio> #include <algorithm&g ...

  9. 离散化+BFS HDOJ 4444 Walk

    题目传送门 /* 题意:问一个点到另一个点的最少转向次数. 坐标离散化+BFS:因为数据很大,先对坐标离散化后,三维(有方向的)BFS 关键理解坐标离散化,BFS部分可参考HDOJ_1728 */ # ...

随机推荐

  1. soeasy的键盘鼠标事件

    在web自动化中,我们可能会遇到需要通过键盘或者鼠标去操作某些元素,那么我们就需要用到键盘事件和鼠标事件了,今天对键盘和鼠标操作进行一个总结 鼠标事件 鼠标事件需要引入ActionChains类,查看 ...

  2. 谷歌学术出现We're sorry解决办法

    出现这个的原因应该是同ip段的或者就是这个ip曾经是个google的黑名单ip,因为恶意爬取谷歌学术了.解决办法就是申请Hurricane Electric Free IPv6 Tunnel Brok ...

  3. Elasticsearch 动态修改replica配置、增删replica

    1. 获取当前所有index配置 curl -XGET http://localhost:9200/_settings 2. 获取某些index的配置 curl -XGET http://localh ...

  4. A9G开发

    巨坑 不支持联通 发短信 上电,准备好后使用 结束标志 取消换行  以16进制发  0x1a 实际测试 中间可以有很多行内容 短信内容结束换行 最后必须以 16进制 发送 0x1a  带和不带换行都可 ...

  5. Linux中关于samba的几个问题

    一.用smbclient命令登录成功但看不了文件 原因:SELinux的阻挡 解决:1.关闭SELinux  :  setenforce 0   (临时生效,重启后失效) 或vi /etc/sysco ...

  6. Softmax与Sigmoid函数的联系

    译自:http://willwolf.io/2017/04/19/deriving-the-softmax-from-first-principles/ 本文的原始目标是探索softmax函数与sig ...

  7. 【JZOJ6239】【20190629】智慧树

    题目 一颗\(n\)个节点的树,每个点有一个权值\(a_i\) 询问树上连通块权值之和对 \(m\) 取模为$ x $ 的方案数 答案对\(950009857\) 取模,满足\(m | 9500098 ...

  8. vue中使用element-ui自定义主题后,vue-cli跑不起来了

    环境:vue-cli 2.x版本 自己在官网配置了主题并放到了项目中https://element.eleme.cn/#/zh-CN/theme 然后,我的脚手架在我的电脑中休息了几天,就跑不通了呢! ...

  9. klass-oop

    (1)Klass Klass 简单来说就是 Java 类在 HotSpot 中的 C++ 对等体,主要用于描述对象实例的具体类型.一般 JVM 在加载 class 文件时,会在方法区创建 Klass ...

  10. hanlp添加词典不起作用

    不起作用的原因很多,这里列举几个 这里的hanlp我虽然用的maven建立的但是要添加自定义词典,所以没有用maven引入的方式,而是下载了data+hanlp.jar文件          data ...