bzoj 3364: [Usaco2004 Feb]Distance Queries 距离咨询 LCA
2017-10-12 11:22
525 查看
→题目链接←
跑最短路果然不可行...
#include<cstdio>
#include<queue>
#include<vector>
using namespace std;
struct node{
int to,len;
node(int x,int y){to=x,len=y;}
};
int n,m;
vector<node>v[40040];
int fa[40040];
int deep[40040];
int par[40040][20];
int dis[40040];
void dfs(int x){
for(int i=0; i<v[x].size(); i++){
int to=v[x][i].to;
int len=v[x][i].len;
if(fa[x]==to)continue;
fa[to]=x;
deep[to]=deep[x]+1;
par[to][0]=x;
dis[to]=dis[x]+len;
dfs(to);
}
}
void cal(){
for(int j=1; (1<<j)<=n; j++){
for(int i=1; i<=n; i++){
par[i][j]=par[par[i][j-1]][j-1];
}
}
}
int lca(int x,int y){
if(deep[x]>deep[y])swap(x,y);
int num=deep[y]-deep[x];
for(int i=18; i>=0; i--){
if(num>=(1<<i)){
y=par[y][i];
num-=(1<<i);
}
}
if(x!=y){
for(int i=18; i>=0; i--){
if(par[x][i]!=par[y][i]){
x=par[x][i];
y=par[y][i];
}
}
x=par[x][0];
}
return x;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=0; i<m; i++){
int x,y,z;
char c[2];
scanf("%d%d%d%s",&x,&y,&z,c);
v[x].push_back(node(y,z));
v[y].push_back(node(x,z));
}
dfs(1);
cal();
int x;
scanf("%d",&x);
for(int i=0; i<x; i++){
int y,z,xx;
scanf("%d%d",&y,&z);
xx=lca(y,z);
printf("%d\n",dis[y]+dis[z]-2*dis[xx]);
}
return 0;
}
【想说的话】
TLE了3遍后我深刻的认识到了....跑最短路果然不可行...
【题解】
显然是裸的LCA【代码】
#include<iostream>#include<cstdio>
#include<queue>
#include<vector>
using namespace std;
struct node{
int to,len;
node(int x,int y){to=x,len=y;}
};
int n,m;
vector<node>v[40040];
int fa[40040];
int deep[40040];
int par[40040][20];
int dis[40040];
void dfs(int x){
for(int i=0; i<v[x].size(); i++){
int to=v[x][i].to;
int len=v[x][i].len;
if(fa[x]==to)continue;
fa[to]=x;
deep[to]=deep[x]+1;
par[to][0]=x;
dis[to]=dis[x]+len;
dfs(to);
}
}
void cal(){
for(int j=1; (1<<j)<=n; j++){
for(int i=1; i<=n; i++){
par[i][j]=par[par[i][j-1]][j-1];
}
}
}
int lca(int x,int y){
if(deep[x]>deep[y])swap(x,y);
int num=deep[y]-deep[x];
for(int i=18; i>=0; i--){
if(num>=(1<<i)){
y=par[y][i];
num-=(1<<i);
}
}
if(x!=y){
for(int i=18; i>=0; i--){
if(par[x][i]!=par[y][i]){
x=par[x][i];
y=par[y][i];
}
}
x=par[x][0];
}
return x;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=0; i<m; i++){
int x,y,z;
char c[2];
scanf("%d%d%d%s",&x,&y,&z,c);
v[x].push_back(node(y,z));
v[y].push_back(node(x,z));
}
dfs(1);
cal();
int x;
scanf("%d",&x);
for(int i=0; i<x; i++){
int y,z,xx;
scanf("%d%d",&y,&z);
xx=lca(y,z);
printf("%d\n",dis[y]+dis[z]-2*dis[xx]);
}
return 0;
}
相关文章推荐
- BZOJ3364: [Usaco2004 Feb]Distance Queries 距离咨询 LCA
- LCA【bzoj3364】 [Usaco2004 Feb]Distance Queries 距离咨询
- LCA【bzoj3364】: [Usaco2004 Feb]Distance Queries 距离咨询
- BZOJ 3364: [Usaco2004 Feb]Distance Queries 距离咨询
- 【BZOJ】3364: [Usaco2004 Feb]Distance Queries 距离咨询
- BZOJ3364: [Usaco2004 Feb]Distance Queries 距离咨询
- 【bzoj3362/3363/3364/3365】[Usaco2004 Feb]树上问题杂烩 并查集/树的直径/LCA/树的点分治
- [BZOJ2684][USACO 2004 FEB]距离询问
- BZOJ_3365_[Usaco2004 Feb]Distance Statistics 路程统计&&POJ_1741_Tree_点分治
- bzoj 3372: [Usaco2004 Feb]Moo University -- Financial Aid 财政补助(set+贪心)
- BZOJ 3367: [Usaco2004 Feb]The Big Game 球赛( dp )
- bzoj3362[Usaco2004 Feb]Navigation Nightmare 导航噩梦
- bzoj 3366: [Usaco2004 Feb]Breeding 奶牛饲育
- BZOJ 3365: [Usaco2004 Feb]Distance Statistics 路程统计 点分治
- BZOJ 3365 [Usaco2004 Feb]Distance Statistics 路程统计 树的点分治
- 树的直径【bzoj3363】 [Usaco2004 Feb]Cow Marathon 奶牛马拉松
- bzoj 3369: [Usaco2004 Feb]Pesky Parentheses 烦人括号(思维题)
- [BZOJ3362][Usaco2004 Feb]Navigation Nightmare 导航噩梦(加权并查集)
- bzoj 3361: [Usaco2004 Jan]培根距离
- bzoj 3362: [Usaco2004 Feb]Navigation Nightmare 导航噩梦 带权并查集