題目鏈接 ：http://acm.hdu.edu.cn/showproblem.php?pid=6670

Mindis
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 548 Accepted Submission(s): 119

Problem Description
平面上有 n 個矩形，矩形的邊平行于坐標軸，現在度度熊需要操控一名角色從 A 點走到 B 點。
該角色可以上下左右移動，在恰被 k 個矩形覆蓋的區域，該角色的速率為 k+1 個距離/秒（矩形覆蓋區域包括邊界）。

請求出 A 移動到 B 最快需要多少秒。

Input
第一行一個整數 T (1≤T≤5) 表示數據組數。
對于每組數據，第一行輸入一個整數 n (1≤n≤200)。
接下來 n 行每行 4 個整數 x1,y1,x2,y2 (0≤x1<x2≤1000000000,0≤y1<y2≤1000000000)，分別表示矩形的左下角和右上角的坐標。
最后一行四個整數 xa,ya,xb,yb ((0≤xa,xb,ya,yb≤1000000000) 代表 A 和 B 的坐標。

Output
對于每組數據，輸出一個小數表示答案。答案保留 5 位小數。

Sample Input
1
1
5 5 6 6
7 7 8 8

Sample Output
2.00000

Source
2019 年百度之星·程序設計大賽 - 初賽一

思路：

先把題目中給出的所有點離散化一下，得到一個n*m 的網格，n和m最大都是400

我們枚舉每一個矩形，同時枚舉離散化之后的矩形坐標中包括哪些點，

這樣能維護出每一個點被多少個矩形包括，

cnt[x][y][k] 表示網格上橫坐標為x，縱坐標為y，方向為k，（0，1，2，3代表上下左右）這條離散化后長度為1的邊被矩形覆蓋的次數。

然后相鄰的點建邊（邊的代價為距離/速度，即時間，距離我們可以通過坐標減去得出，速度由該點被多少個矩形包括得出），跑圖的最短路徑即可，

代碼寫起來很長，其實不難寫，因為寫很多部分是不需要動腦的機械操作。，注意細節即可，

細節見代碼：

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> #include <queue> #include <stack> #include <map> #include <set> #include <vector> #include <iomanip> #define ALL(x) (x).begin(), (x).end() #define sz(a) int(a.size()) #define all(a) a.begin(), a.end() #define rep(i,x,n) for(int i=x;i<n;i++) #define repd(i,x,n) for(int i=x;i<=n;i++) #define pii pair<int,int> #define pll pair<long long ,long long> #define gbtb ios::sync_with_stdio(false),cin.tie(0),cout.tie(0) #define MS0(X) memset((X), 0, sizeof((X))) #define MSC0(X) memset((X), '\0', sizeof((X))) #define pb push_back #define mp make_pair #define fi first #define se second #define eps 1e-6 #define gg(x) getInt(&x) #define chu(x) cout<<"["<<#x<<" "<<(x)<<"]"<<endl using namespace std; typedef long long ll; ll gcd(ll a, ll b) {return b ? gcd(b, a % b) : a;} ll lcm(ll a, ll b) {return a / gcd(a, b) * b;} ll powmod(ll a, ll b, ll MOD) {ll ans = 1; while (b) {if (b % 2)ans = ans * a % MOD; a = a * a % MOD; b /= 2;} return ans;} inline void getInt(int* p); const int maxn = 400010; const ll inf = 1e18; /*** TEMPLATE CODE * * STARTS HERE ***/ int t; int n; struct juzhen {int x1,x2,y1,y2; }a[500]; std::vector<int> vx,vy; int sn,sm; int getid(int x,int y) {return (x-1)*sn+y; } double dis[maxn]; struct node {int to;double val;node (){}node(int tt,double vv){to=tt;val=vv;}bool operator <(const node &b)const{return val>b.val;} }; std::vector<node> v[maxn]; int cnt[500][500][5]; bool check(int x,int y) {return x>=1&&x<=sn&&y>=1&&y<=sm; } juzhen aid; priority_queue<node> heap; void dij() {while(heap.size()){heap.pop();}dis[getid(aid.x1,aid.y1)]=0;heap.push(node(getid(aid.x1,aid.y1),0));node temp;while(!heap.empty()){temp=heap.top();heap.pop();for(auto x:v[temp.to]){if(dis[x.to]>x.val+dis[temp.to]){dis[x.to]=x.val+dis[temp.to];heap.push(node(x.to,dis[x.to]));}}} } int main() {//freopen("D:\\code\\text\\input.txt","r",stdin);//freopen("D:\\code\\text\\output.txt","w",stdout);gbtb;cin>>t;while(t--){cin>>n;vx.clear();vy.clear();repd(i,1,n){cin>>a[i].x1>>a[i].y1>>a[i].x2>>a[i].y2;vx.push_back(a[i].x1);vx.push_back(a[i].x2);vy.push_back(a[i].y1);vy.push_back(a[i].y2);}cin>>aid.x1>>aid.y1>>aid.x2>>aid.y2;vx.push_back(aid.x1);vx.push_back(aid.x2);vy.push_back(aid.y1);vy.push_back(aid.y2);sort(ALL(vx));sort(ALL(vy));vx.erase(unique(ALL(vx)),vx.end());vy.erase(unique(ALL(vy)),vy.end());sn=sz(vx);sm=sz(vy);repd(i,1,sn){repd(j,1,sm){dis[getid(i,j)]=inf;v[getid(i,j)].clear();repd(z,0,3)cnt[i][j][z]=1;}}aid.x1=lower_bound(ALL(vx),aid.x1)-vx.begin()+1;aid.x2=lower_bound(ALL(vx),aid.x2)-vx.begin()+1;aid.y1=lower_bound(ALL(vy),aid.y1)-vy.begin()+1;aid.y2=lower_bound(ALL(vy),aid.y2)-vy.begin()+1;repd(i,1,n){int x1=lower_bound(ALL(vx),a[i].x1)-vx.begin()+1;int y1=lower_bound(ALL(vy),a[i].y1)-vy.begin()+1;int x2=lower_bound(ALL(vx),a[i].x2)-vx.begin()+1;int y2=lower_bound(ALL(vy),a[i].y2)-vy.begin()+1;repd(j,x1+1,x2-1){repd(k,y1+1,y2-1){repd(z,0,3)cnt[j][k][z]++;}}// 0 1 2 3// 上 下 左 右repd(j,x1+1,x2-1){cnt[j][y2][1]++;cnt[j][y2][2]++;cnt[j][y2][3]++;cnt[j][y1][0]++;cnt[j][y1][2]++;cnt[j][y1][3]++;}repd(j,y1+1,y2-1){cnt[x1][j][0]++;cnt[x1][j][1]++;cnt[x1][j][3]++;cnt[x2][j][0]++;cnt[x2][j][1]++;cnt[x2][j][2]++;}cnt[x1][y1][0]++;cnt[x1][y1][3]++;cnt[x2][y1][0]++;cnt[x2][y1][2]++;cnt[x1][y2][3]++;cnt[x1][y2][1]++;cnt[x2][y2][1]++;cnt[x2][y2][2]++;}repd(i,1,sn){repd(j,1,sm){if(check(i-1,j)){v[getid(i,j)].push_back(node(getid(i-1,j),1.00*(vx[i-1]-vx[i-2])/cnt[i][j][2]));}if(check(i+1,j)){v[getid(i,j)].push_back(node(getid(i+1,j),1.00*(vx[i]-vx[i-1])/cnt[i][j][3]));}if(check(i,j-1)){v[getid(i,j)].push_back(node(getid(i,j-1),1.00*(vy[j-1]-vy[j-2])/cnt[i][j][1]));}if(check(i,j+1)){v[getid(i,j)].push_back(node(getid(i,j+1),1.00*(vy[j]-vy[j-1])/cnt[i][j][0]));}}}dij(); // cout<<getid(aid.x2,aid.y2)<<endl;cout<<fixed<<setprecision(5)<<dis[getid(aid.x2,aid.y2)]<<endl;}return 0; }inline void getInt(int* p) {char ch;do {ch = getchar();} while (ch == ' ' || ch == '\n');if (ch == '-') {*p = -(getchar() - '0');while ((ch = getchar()) >= '0' && ch <= '9') {*p = *p * 10 - ch + '0';}}else {*p = ch - '0';while ((ch = getchar()) >= '0' && ch <= '9') {*p = *p * 10 + ch - '0';}} }

轉載于:https://www.cnblogs.com/qieqiemin/p/11391631.html

總結

以上是生活随笔為你收集整理的2019 年百度之星·程序设计大赛 - 初赛一 C. HDU 6670 Mindis 离散化+dijkstra的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。