题解 P2738 【[USACO4.1]篱笆回路Fence Loops】
这题是我期中测试的一题水题,然而英文题目太长了不想读...后面考完被同学提醒后20分钟切了(心塞)
切完看了波题解,发现貌似我的方法跟大家都不一样呢...
常规做法: \(Floyd\)
这个有三页的题解了,本蒟蒻在此不多讲
时间复杂度 \(O(n^3)\)
我的解法: \(queue+dijkstrra\)
可以发现,如果一个点绕了一圈后跑回远点,那么他所跑的路就是一圈周长
用这个思路,我们暴力枚举每个点跑 \(dijkstra\) .每次绕了一圈回到原点的时候(肯定比原点要大),我们更新答案,不更新原点(本来也不需要更新)
问题来了, \(dijkstra\) 不会从自己一条边出发,再从对吗跑回自己呢?
我们可以观察到一个点:左边端点所连的边跟右边端点所连的边没有任何关系.基于此,不管我们在点 \(i\) 的哪一个端点,我们只要确认去的点 \(j\) 跟点 \(i\) 所连接的端点的位置.例如我们在他的左端点,我们只将右端点放入队列.反之依然.(记得不要在端点之间建边)
int to = 1;//假设他在右端点(1为右端点,0为左端点)
for (int k : adj[v][0]) if (k==qf) goto abcd;假设在左端点找到,那么就之间下一步(因为如果他在左端点,他要去的地方是右端点,to表示的不是他现在的地方,而是要去的地方)
to = 0;//如果搜遍左端点都没找到,那么证明他现在在右端点,要去左端点
abcd:;
if (dist[v][to]>dist[qf][qs]+len[v]){//裸的dijkstra
dist[v][to] = dist[qf][qs]+len[v];
if (!inq[v][to]){inq[v][to] = true;q.push(make_pair(v,to));}
}
最后的答案更新方法:(在 \(dijkstra\) 里面更新)
if (v==pos && to==curr) {
ans = min(ans,dist[qf][qs]+len[pos]);
}
完整代码:
#include <iostream>
#include <algorithm>
#include <queue>
#include <vector>
using namespace std;
#define pp pair<int,int>
#define f first
#define s second
int n;
int len[105],ans = 1e9;
int dist[105][2];
bool inq[105][2];
vector<int> adj[105][2];
void dfs(int pos,int curr){
queue<pp> q;
memset(dist,0x3f3f,sizeof(dist));
dist[pos][curr] = 0;
q.push(make_pair(pos,curr));//起点(记得将起点状态放进去
while(!q.empty()){
int qf = q.front().f,qs = q.front().s;q.pop();
inq[qf][qs] = false;
for (int v : adj[qf][qs]){
int to = 1;
for (int k : adj[v][0]) if (k==qf) goto abcd;
to = 0;
abcd:;//上面讲的转移方式
if (dist[v][to]>dist[qf][qs]+len[v]){
dist[v][to] = dist[qf][qs]+len[v];
if (!inq[v][to]){inq[v][to] = true;q.push(make_pair(v,to));}
}
if (v==pos && to==curr) {
ans = min(ans,dist[qf][qs]+len[pos]);
}//如果他现在要去原点,那么更新答案
}
}
}
int main(){
cin >> n;
for (int i=0;i<n;i++){
int a,b,c,d; cin >> a >> b >> c >> d;
len[a] = b;
for (int j=0;j<c;j++) {
int t; cin >> t;
adj[a][0].push_back(t);
}
for (int j=0;j<d;j++){
int t; cin >> t;
adj[a][1].push_back(t);
}//分开左右端点建边
}
for (int i=1;i<=n;i++) {dfs(i,0);dfs(i,1);}
cout << ans;
}
复杂度仍然是 \(O(n^3)\)
为什么呢?其实就是一个小地方:在寻找左右端点的时候最坏情况会做n次.这个转移用数组可以优化成 \(O(1)\) .加两行代码后可以变成 \(O(n^2)\) 然而过了就懒得改了
接下来将那位Java大佬的思路了(已得到授权):
大佬的地址
其实区别也不大,就是将 \(queue\) 改成了 \(priority\) \(queue\) ,再加上一些小细节的区别
存图方法: 用 \(hashset\) 来代替数组,保证能在 \(O(logn)\) 的速度找到这个数
更新答案方式:他将每个点的距离改为 \(inf\) (包括原点),在原点出发时直接更新而不是取最小值.跑完之后只需要取原点的距离就是最终答案
求是否在端点上:
if(j==i) {//去左端点
if(!seg[i].hs2.contains(e.v))continue;//这个点不在右端点就不做
dist[i] = Math.min(dist[i], seg[i].l+e.w);}//看看自己能不能更新
if(seg[j].l+e.w < dist[j]) {
dist[j] = seg[j].l + e.w;
pq.add(new Status(e.v,j,dist[j]));
}//pq更新下一个数
}
完整代码:
import java.io.*;
import java.util.*;
public class Main {
private static StreamTokenizer st;
private static int nextInt()throws IOException{
st.nextToken();
return (int)st.nval;
}
public static void main(String args[]) throws IOException{
st = new StreamTokenizer(new BufferedReader(new InputStreamReader(System.in)));
int N = nextInt();
Segment[] seg = new Segment[N];
for(int i = 0; i < N; ++i) {
int s = nextInt()-1, l = nextInt(), n1 = nextInt(), n2 = nextInt();
seg[s] = new Segment(l);
for(int j = 0; j < n1; ++j)seg[s].hs1.add(nextInt()-1);
for(int j = 0; j < n2; ++j)seg[s].hs2.add(nextInt()-1);
}//记图
//bfs
PriorityQueue<Status> pq = new PriorityQueue<>();
int ans = Integer.MAX_VALUE;
for(int i = 0; i < N; ++i) {
//go from hs1 and back to hs2
int dist[] = new int[N];
Arrays.fill(dist, Integer.MAX_VALUE);
for(int j : seg[i].hs1) {
pq.add(new Status(i,j,seg[j].l));
dist[j] = seg[j].l;
}
while(!pq.isEmpty()) {
Status e = pq.poll();
if(e.w != dist[e.v])continue;
if(seg[e.v].hs1.contains(e.u)) {
//往右端点走的情况
for(int j : seg[e.v].hs2) {
if(j==i) {
if(!seg[i].hs2.contains(e.v))continue;
dist[i] = Math.min(dist[i], seg[i].l+e.w);
}
if(seg[j].l+e.w < dist[j]) {
dist[j] = seg[j].l + e.w;
pq.add(new Status(e.v,j,dist[j]));
}
}
}else {
//往左端点走的情况
for(int j : seg[e.v].hs1) {
if(j==i) {
if(!seg[i].hs2.contains(e.v))continue;
dist[i] = Math.min(dist[i], seg[i].l+e.w);
}
if(seg[j].l+e.w < dist[j]) {
dist[j] = seg[j].l + e.w;
pq.add(new Status(e.v,j,dist[j]));
}
}
}
}
ans = Math.min(ans, dist[i]);
}
PrintWriter pw = new PrintWriter(new BufferedWriter(new OutputStreamWriter(System.out)));
pw.println(ans);
pw.close();
}
private static class Segment{
int l;
HashSet<Integer> hs1, hs2;
public Segment(int l) {
this.l = l;
hs1 = new HashSet<>();
hs2 = new HashSet<>();
}
}
private static class Status implements Comparable<Status>{
int u, v, w;
public Status(int u, int v, int w) {
this.u = u;
this.v = v;
this.w = w;
}
@Override
public int compareTo(Status o) {
return w-o.w;
}
}
}
复杂度 \(O(n^2log^2n)\) 第一个 $log $在 \(pq\) ,第二个 \(log\) 在 \(set\).然而由于 \(priority\) \(queue\) \(dijkstra\) 的性质,在正常情况下跑不满这个时间,预估时间 \(O(n^2)\).
然而这位大佬用的java,于是在洛谷时间被虐的体无完肤
题解 P2738 【[USACO4.1]篱笆回路Fence Loops】的更多相关文章
- 洛谷P2738 [USACO4.1]篱笆回路Fence Loops
P2738 [USACO4.1]篱笆回路Fence Loops 11通过 21提交 题目提供者该用户不存在 标签USACO 难度提高+/省选- 提交 讨论 题解 最新讨论 暂时没有讨论 题目描述 ...
- 洛谷 P2738 [USACO4.1]篱笆回路Fence Loops
P2738 [USACO4.1]篱笆回路Fence Loops 题目描述 农夫布朗的牧场上的篱笆已经失去控制了.它们分成了1~200英尺长的线段.只有在线段的端点处才能连接两个线段,有时给定的一个端点 ...
- [USACO4.1]篱笆回路Fence Loops
题目:USACO Training 4.1(在官网上提交需加文件输入输出).洛谷P2738. 题目大意:给你一张图里的边集,让你求出这张图的最小环. 解题思路:求最小环很简单,用Floyd即可.最重要 ...
- USACO 4.1 Fence Loops(Floyd求最小环)
Fence Loops The fences that surround Farmer Brown's collection of pastures have gotten out of contro ...
- USACO 4.1 Fence Loops
Fence Loops The fences that surround Farmer Brown's collection of pastures have gotten out of contro ...
- 编程算法 - 篱笆修理(Fence Repair) 代码(C)
篱笆修理(Fence Repair) 代码(C) 本文地址: http://blog.csdn.net/caroline_wendy 题目: 把一块木板切成N块, 每次切两块, 分割的开销是木板长度, ...
- USACO4.1 Fence Loops【最小环&边->点转化】
数据不是很大,如果要转换为正常的那种建图方式的话,可以给点进行标号,用一个二维数组存这两条边相交的那个点的标号,方便处理.一定要注意不要同一个点使用不同的编号也不要不同的点使用同一个编号(这不是废话嘛 ...
- USACO Section 4
前言 好久没更新这个系列了,最近闲的无聊写一下.有两题搜索懒得写了. P2737 [USACO4.1]麦香牛块Beef McNuggets https://www.luogu.com.cn/probl ...
- 浅谈Floyd的三种用法 By cellur925
Floyd大家可能第一时间想到的是他求多源最短路的n³算法.其实它还有另外两种算法的嘛qwq.写一发总结好了qwq. 一.多源最短路 放段代码跑,注意枚举顺序,用邻接矩阵存图.本质是一种动规. 复杂度 ...
随机推荐
- ES6 之 Object 的方法总结
阅读:Object 1.obj的"." 或 "[]"方法 读取对象的属性或方法 对象属性的读取:ES6中被Proxy的get(target, propKey, ...
- python 第一节 脚本 import from reload exec
环境Ubuntu 14.04, 不写交互式命令行了,直接脚本开始. # first Python script import sys print(sys.platform) print(2**4) x ...
- 联系我们地图坐标展示js
<script type="text/javascript" src="http://api.map.baidu.com/api?v=2.0&ak=6d88 ...
- (HN)AHOI2018 转盘
题意: 有 \(n\) 个格子围成一圈,每个格子里有一个物品,每个物品的出现时间为 \(T_i\) .开始时选择一个格子为起点,每个单位时间可以向前走一格或不动,若当前格的物品已出现则将其标记.有 \ ...
- Inception Score
转载 https://www.jiqizhixin.com/articles/2019-01-10-18 全面解析Inception Score原理及其局限性 https://blog.csdn ...
- 使用BurpSuite和Hydra爆破相关的服务(9.25 第十一天)
使用BP和Hydra爆破相关的服务 Hydra:九头蛇,开源的功能强大的爆破工具,支持的服务有很多,使用Hydra爆破C/S架构的服务. 使用BurpSuite爆破web服务 DVWA:web应用程序 ...
- vue坑 - vue安装vue-cli报错coffee-script@1.12.7: CoffeeScript on NPM has moved to "coffeescript" (no hyphen)或者说不支
$ npm install -g vue-cli npm WARN deprecated coffee-script@1.12.7: CoffeeScript on NPM has moved to ...
- C++ #if 1
当注释掉大块代码时,使用"#if 0"比使用"/**/"要好,因为用"/**/"做大段的注释要防止被注释掉的代码中有嵌套的"/** ...
- jvm 优化
1.jvm 运行参数 1.1 jvm 参数类型 1.1.1 标准参数 -help -version 1.1.2 -X参数(非标准参数) -Xint -Xcomp 1.1.3-XX参数(使用率高) -X ...
- docker-compose 安装 mongodb
1. 修改 docker-compose.yml version: "2.1" services: php7.1: build: ./php image: php7.1-ext p ...