NOIP 2017 宝藏 包括部分分做法!

本文深入探讨了NOIP中一道名为“宝藏”的题目,详细分析了20分、40分、70分及100分的不同解题策略,包括最短路径算法的应用和优化技巧。

题目链接 P3959 宝藏

一个小小的开头

对于这道题本蒟讲讲自己的一些想法:分别是20分做法、40分做法、70分做法(暂时还没做出,会补上的)、和100分做法。
其实现在的我很是不明白当初(去年高二时)参加noip时这一道题为什么没有拿到45分?我真的搞不懂当初自己为啥为那么愚蠢(也许现在也差不多)n明明很小(当初没怎么接触过状压,所以就不晓得,这情有可原)而且明明很多边是没用的!但我还因为考场一紧张就不记得dij(我也不记得当时自己是不是想的prim)然后就觉得可以用floyd求任意点的最短路(可我当时不会dij也应该会SPFA啊,当时自己还特喜欢用SPFA)。然后因为45分以内的边长是相同的,所以我们就可以去乱搞一下,但我tema还是只得了25分!!!当时的自己是有多菜!(可能现在也好不到哪去)

20分做法

个人觉得这是给一些不知道最短路是啥的OIer的分数,求一下树上的路径,然后因为边长都一样所以最后答案就是每个点到跟节点深度*边长值就行了。不过还要以每一个节点为跟节点去算一次,球个最小值。很容易写,就没写代码了。

40分做法

根据20分做法带来的启发,我们只要对这个图以某一个点为起始点跑最短路(不管是dij还是floyd还是spfa都可以,只是dij将所有的都可以一次性求出来)然后对于没一个点我们就可像上述20分的做法去求答案了。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
int n,m,ans=0x3f3f3f3f,sum,vv,pos;
int a[15][15],dist[15],vis[15];
inline int read() {
    int x=0,w=1;char ch;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+(ch-48),ch=getchar();
    return x*w;
}
void dij(int s) {
    memset(vis,0,sizeof(vis));sum=0;
    memset(dist,0x3f,sizeof(dist));
    for(register int i=1;i<=n;++i) dist[i]=a[s][i];
    vis[s]=1;
    for(register int i=1;i<=n-1;++i) {
        int minn=0x3f3f3f3f;
        for(register int j=1;j<=n;++j) 
            if(!vis[j]&&dist[j]<minn) minn=dist[j],pos=j;
        vis[pos]=1;sum=sum+dist[pos]*vv;//边求就边计算答案了
        for(register int j=1;j<=n;++j) 
            if(!vis[j]&&dist[j]>a[pos][j]+dist[pos]) 
            	dist[j]=min(dist[j],dist[pos]+a[pos][j]);
    }
}
int main() {
    n=read();m=read();
    memset(a,0x3f,sizeof(a));
    for(register int i=1;i<=m;++i) {
        int u=read(),v=read(),w=read();vv=w;
        a[u][v]=a[v][u]=1;
    }
    for(register int i=1;i<=n;++i) {
        prim(i);
        ans=min(ans,sum);
    }
    cout<<ans<<endl;
    return 0;
}

70分做法

暂时还没作出,会尽快补上,请谅解。

100分做法

网上很多大佬都讲了,贴一个博客吧。https://www.luogu.org/blog/dazade8/solution-p3959

# P3959 [NOIP 2017 提高组] 宝藏 ## 目背景 NOIP2017 D2T2 ## 目描述 参与考古挖掘的小明得到了一份藏宝图,藏宝图上标出了 $n$ 个深埋在地下的宝藏屋, 也给出了这 $n$ 个宝藏屋之间可供开发的 $m$ 条道路和它们的长度。 小明决心亲自前往挖掘所有宝藏屋中的宝藏。但是,每个宝藏屋距离地面都很远,也就是说,从地面打通一条到某个宝藏屋的道路是很困难的,而开发宝藏屋之间的道路则相对容易很多。 小明的决心感动了考古挖掘的赞助商,赞助商决定免费赞助他打通一条从地面到某个宝藏屋的通道,通往哪个宝藏屋则由小明来决定。 在此基础上,小明还需要考虑如何开凿宝藏屋之间的道路。已经开凿出的道路可以任意通行不消耗代价。每开凿出一条新道路,小明就会与考古队一起挖掘出由该条道路所能到达的宝藏屋的宝藏。另外,小明不想开发无用道路,即两个已经被挖掘过的宝藏屋之间的道路无需再开发。 新开发一条道路的代价是 $L\times K$。其中 $L$ 代表这条道路的长度,$K$ 代表从赞助商帮你打通的宝藏屋到这条道路起点的宝藏屋所经过的宝藏屋的数量(包括赞助商帮你打通的宝藏屋和这条道路起点的宝藏屋) 。 请你编写程序为小明选定由赞助商打通的宝藏屋和之后开凿的道路,使得工程总代价最小,并输出这个最小值。 ## 输入格式 第一行两个用空格离的正整数 $n,m$,代表宝藏屋的个数和道路数。 接下来 $m$ 行,每行三个用空格离的正整数,别是由一条道路连接的两个宝藏屋的编号(编号为 $1-n$),和这条道路的长度 $v$。 ## 输出格式 一个正整数,表示最小的总代价。 ## 输入输出样例 #1 ### 输入 #1 ``` 4 5 1 2 1 1 3 3 1 4 1 2 3 4 3 4 1 ``` ### 输出 #1 ``` 4 ``` ## 输入输出样例 #2 ### 输入 #2 ``` 4 5 1 2 1 1 3 3 1 4 1 2 3 4 3 4 2 ``` ### 输出 #2 ``` 5 ``` ## 说明/提示 ![](https://cdn.luogu.com.cn/upload/pic/10868.png) 【样例解释 $1$】 小明选定让赞助商打通了 $1$ 号宝藏屋。小明开发了道路 $1 \to 2$,挖掘了 $2$ 号宝藏。开发了道路 $1 \to 4$,挖掘了 $4$ 号宝藏。还开发了道路 $4 \to 3$,挖掘了 $3$ 号宝藏。 工程总代价为 $1 \times 1 + 1 \times 1 + 1 \times 2 = 4 $。 【样例解释 $2$】 小明选定让赞助商打通了 $1$ 号宝藏屋。小明开发了道路 $1 \to 2$,挖掘了 $2$ 号宝藏。开发了道路 $1 \to 3$,挖掘了 $3$ 号宝藏。还开发了道路 $1 \to 4$,挖掘了 $4$ 号宝藏。 工程总代价为 $1 \times 1 + 3 \times 1 + 1 \times 1 = 5$。 【数据规模与约定】 对于 $ 20\%$ 的数据: 保证输入是一棵树,$1 \le n \le 8$,$v \le 5\times 10^3$ 且所有的 $v$ 都相等。 对于 $40\%$ 的数据: $1 \le n \le 8$,$0 \le m \le 10^3$,$v \le 5\times 10^3$ 且所有的 $v$ 都相等。 对于 $ 70\%$ 的数据: $1 \le n \le 8$,$0 \le m \le 10^3$,$v \le 5\times 10^3$。 对于 $ 100\%$ 的数据: $1 \le n \le 12$,$0 \le m \le 10^3$,$v \le 5\times 10^5$。 --- **根据目,修改我代码中错误部分** #include<bits/stdc++.h> #define int unsigned long long #define baka; cout<<p.v<<" "<<p.w<<" "<<p.time<<endl; using namespace std; int n, m ,ans = 0x3f3f3f3f3f3f3f3f3f3f3f3f3f3f3f3f,vis; struct Node{ int v, w, time, b; inline bool operator > (const Node &p) const{ if(w < 50*p.w) return w > p.w; return b > p.b; } }; vector<bool> vi; vector<int> head,nxt; vector<Node> to; priority_queue<Node, vector<Node>, greater<Node> > q; inline void add(int u,int v,int w){ Node p; p.v = v,p.w = w; nxt.push_back(head[u]), head[u] = to.size(), to.push_back(p); } void bfs(){ int temp = 0; Node p, pp; while(!q.empty()){ p = q.top(), q.pop(); if(vi[p.v]) continue; // baka; vi[p.v] = 1, vis++, temp += p.w, p.time ++; if(vis == n){ ans = min(temp, ans); priority_queue<Node, vector<Node>, greater<Node> > empty; swap(empty, q); return; } for(int i=head[p.v]; ~i; i = nxt[i]){ pp.v = to[i].v, pp.w = p.time * to[i].w; pp.time = p.time; pp.b = rand() % (n*3); q.push(pp); } } } signed main(){ ios::sync_with_stdio(0),cin.tie(0),cout.tie(0); srand(time(0)); cin>>n>>m; int u, v, w; vi.resize(n+1, 0), head.resize(n+1, -1); while(m--){ cin>>u>>v>>w; add(u,v,w); add(v,u,w); } Node p; for(int i=1;i<=n;i++){ for(int j=1;j<=1300;j++){ vis = 0, p.v = i,p.w = 0,p.time = 0; p.b = rand() % (n*3); q.push(p); bfs(); vi.clear(); vi.resize(n+1,0); } } cout<<ans; return 0; }
最新发布
11-26
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值