题解 [NOIP2017 提高组]宝藏
这是蓝书上状压的例题啊,怎么会出现在模拟赛里
不过就算原题我也没把握写对
核心思路:
先令\(dp[s]\)为当前状态为\(s\)时的总花费最小值,\(cnt[s][i]\)为这个方案中由根节点(赞助商打通的节点)到\(i\)节点最少经过的房间数
此题转移较为麻烦,转移时需要知道当前每个点经过了几个房间
考虑对每个状态开个数组记录到此状态每个点经过了几个房间
注意这个\(cnt\)数组是\(dp\)值的附属物,所以\(dp\)值更新的时候\(cnt[]\)需要跟着覆写过去
然后就是一个状压求最小值的板子了
怎么可能这么简单,这里的转移覆写的后效性怎么说?
然而事实是,这里每次转移取最小值更新根本没有后效性
所以考场上直接猜个结论开始写可以有100pts拿 (弥天大雾
好了这个后效性的有无是此题难点
因为如果真的有后效性的话我们并不知道如何消除它,所以我们试着去证明这里不必处理后效性,杀敌于无形之中
首先考虑我们所说这里的后效性是什么:
放张图:
令赞助商打通的是节点1
那么对于这个节点3,它会被1-2-3这条路径更新,\(cnt=3\)
而1-3这条路径虽然权值更大,但其\(cnt=2\)
对于后续的更新无法判断谁更优,
那这里直接由1-2-3把3更新掉不就有后效性了吗?
然而我们更新3的实质是\(dp[0111] = min(dp[0011], dp[0101])\)
我们所想象的后效性是在节点1,2,3均打通后向4转移会有后效性
谁说我们一定要先打通1,2,3才能向4转移了?
考虑节点4的更新,\(dp[1111] = min(dp[0111], dp[1101])\)(其它状态不合法,已舍去)
由方程知这里可以也考虑了1-3-4-2的方案,推广到所有情况,会考虑少一条边的所有可能,并没有后效性
同样由方程知,我们在状压时实际考虑了所有点向这个点所有可能的连边方案,
而我们所想的「后效性」,实质上是「在已固定了几条连边后出现的后效性」
则我们直接列出方程
\]
而这里的cnt数组是dp值的附属物,在转移的同时要把\(cnt[s]\)覆写过去
然后按方程转移即可
哦对了还有这个题里全是重边,建议开邻接矩阵存边
Code:
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define ld long double
#define usd unsigned
#define ull unsigned long long
//#define int long long
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf, 1, 1<<21, stdin)), p1==p2?EOF:*p1++)
char buf[1<<21], *p1=buf, *p2=buf;
inline int read() {
int ans=0, f=1; char c=getchar();
while (!isdigit(c)) {if (c=='-') f=-f; c=getchar();}
while (isdigit(c)) {ans=(ans<<3)+(ans<<1)+(c^48); c=getchar();}
return ans*f;
}
int n, m;
ll dp2[13][1<<14], val[15];
int cnt[13][1<<14][13];
int head[20], size;
int mp[15][15];
inline void add(int s, int t, int w) {mp[s][t]=min(mp[s][t], w); mp[t][s]=min(mp[t][s], w);}
signed main()
{
#ifdef DEBUG
freopen("1.in", "r", stdin);
#endif
memset(dp2, 127, sizeof(dp2));
memset(mp, 127, sizeof(mp));
n=read(); m=read();
int lim=1<<n;
for (int i=1,u,v,w; i<=m; ++i) {u=read(); v=read(); w=read(); add(u, v, w); add(v, u, w);}
for (int begin=0; begin<n; ++begin) {
dp2[begin][1<<begin] = 0;
cnt[begin][1<<begin][begin] = 1;
ll *dp=dp2[begin];
for (int s=1; s<lim; ++s) {
for (int j=0; j<n; ++j) {
if (!(s&(1<<j))) continue;
for (int i=1,v; i<=n; ++i) {
if (j+1==i || mp[j+1][i]>500010) continue;
v = i-1;
if (s&(1<<v)) continue;
if (dp[s]+1ll*mp[j+1][i]*(cnt[begin][s][j]) < dp[s|(1<<v)]) {
dp[s|(1<<v)] = min(dp[s|(1<<v)], dp[s]+1ll*mp[j+1][i]*(cnt[begin][s][j]));
for (int k=0; k<n; ++k) cnt[begin][s|(1<<v)][k] = cnt[begin][s][k];
cnt[begin][s|(1<<v)][v] = cnt[begin][s][j]+1;
}
}
}
}
}
ll ans=INF;
for (int i=0; i<n; ++i) ans = min(ans, dp2[i][lim-1]);
printf("%lld\n", ans);
return 0;
}
题解 [NOIP2017 提高组]宝藏的更多相关文章
- NOIP2017[提高组] 宝藏 题解
解析 我们观察范围可以发现n非常的小,(一般来说不是搜索就是状压dp)所以说对于这题我们可以用记忆化搜索或者dp,我们发现起点不同那么最终答案也就不同,也就是说答案是跟起点有关的,于是我们便可以想到去 ...
- [NOIp2017提高组]宝藏
#include<cstdio> #include<cctype> #include<algorithm> inline int getint() { regist ...
- [NOIP2017 提高组] 宝藏
考虑到这种对于某种操作顺序有一个权值. 且这个权值有一个\(O(n)\)或者更好的复杂度求出. 求最值. 那可以用模拟退火. #include<iostream> #include< ...
- 【题解】NOIP2017 提高组 简要题解
[题解]NOIP2017 提高组 简要题解 小凯的疑惑(数论) 不讲 时间复杂度 大力模拟 奶酪 并查集模板题 宝藏 最优解一定存在一种构造方法是按照深度一步步生成所有的联通性. 枚举一个根,随后设\ ...
- [NOIp2017提高组]列队
[NOIp2017提高组]列队 题目大意 一个\(n\times m(n,m\le3\times10^5)\)的方阵,每个格子里的人都有一个编号.初始时第\(i\)行第\(j\)列的编号为\((i-1 ...
- JZOJ 5196. 【NOIP2017提高组模拟7.3】B
5196. [NOIP2017提高组模拟7.3]B Time Limits: 1000 ms Memory Limits: 262144 KB Detailed Limits Goto Pro ...
- JZOJ 5197. 【NOIP2017提高组模拟7.3】C
5197. [NOIP2017提高组模拟7.3]C Time Limits: 1000 ms Memory Limits: 262144 KB Detailed Limits Goto Pro ...
- JZOJ 5195. 【NOIP2017提高组模拟7.3】A
5195. [NOIP2017提高组模拟7.3]A Time Limits: 1000 ms Memory Limits: 262144 KB Detailed Limits Goto Pro ...
- JZOJ 5184. 【NOIP2017提高组模拟6.29】Gift
5184. [NOIP2017提高组模拟6.29]Gift (Standard IO) Time Limits: 1000 ms Memory Limits: 262144 KB Detailed ...
随机推荐
- Elasticsearch-02-入门:集群、节点、分片、索引及常用API
2. 基础入门 2.1 重要概念 2.1.1 集群和节点 1)cluster Elasticsearch集群是由一个或多个节点组成,通过其集群名称来进行唯一标识.节点在搜索到集群之后,通过判断自身的 ...
- 【笔记】Python编程 从入门到实践 第二版(基础部分)
1 字符串相关函数 .title() # 将字符串每个单词的首字母大写 .upper() #不改变字符串变量的值 .lower() #不改变字符串变量的值 f"{var} ,字符串" ...
- Django基础08篇 filter&tag
1.Django自带的过滤器filter views.py中代码 def template_tags(request): import datetime content = '三胖content三胖c ...
- MySQL服务不见 - 解决方法
因为开发需要,今天安装了PHPStudy服务.导致以前的MySQL服务在服务表里面不见了.通过查阅网址的资料解决了,那么赶快记录下来 1. 确认当前的系统是管理员身份 2. 切换到MySQL数据库的安 ...
- 正则表达式的模式匹配----V客学院技术分享
正则表达式是由一个字符序列形成的搜索模式. 你在文本中搜索数据时,你可以用搜索模式来描述你要查询的内容. 正则表达式可以是一个简单的字符,或一个更复杂的模式. 正则表达式可用于所有文本搜索和文本替换的 ...
- SpringCloud学习之【Eureka实现服务注册与发现】
这段时间开始整理之前的SpringCloud实践笔记,这里感谢翟永超大佬的文章SpringCloud从入门到精通的指导. 项目结构 服务注册中心 注意: 1.SpringCloud与SpringBoo ...
- SSM中如何上传图片
1.文件配置 2.jsp页面 文件的name值不能跟数据库列名一致 3.控制层收集数据转发到逻辑层 4.逻辑层处理把用户信息存到数据库 5.注册成功后跳到jsp页面进行展示
- CreateWindow() -- 创建普通的窗口
(1)函数原型 1 HWND CreateWindow( 2 LPCTSTR lpClassName, //pointer to register class name 3 LPCTSTR lpWin ...
- 带标签的for循环
for循环可以加标签,使用break或者continue时,若存在多层嵌套循环可指定标签的for循环 public class ForLabel { public static void main(S ...
- LinkedHashMap 的实现原理
LinkedHashMap 概述 HashMap 是无序的,HashMap 在 put 的时候是根据 key 的 hashcode 进行 hash 然后放入对应的地方.所以在按照一定顺序 put 进 ...