codevs 1835 魔法猪学院 A*寻k短路做了一个月卡死在spfa那了/(ㄒoㄒ)/~~
SPFA时点出队后一定要把在队内的标记置为false!SPFA时点出队后一定要把在队内的标记置为false!SPFA时点出队后一定要把在队内的标记置为false!
我因为这个卡了一个月大家信吗?测得时候总是60分,丢情况,总之我太弱了
一个月前傻傻的手写堆o(>﹏<)o,现在写Yveh教的我的stl大法的优先队列用于A*爽的飞起\(≧▽≦)/
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int point[5003],next[200003],v[200003],point2[5003],next2[200003],v2[200003],cnt=0,cnt2=0,N,M;
double c[200003],c2[200003],dist[5003],E;
bool p[5003];
struct datatype{
double g; int isk;
bool operator < (const datatype &a) const
{
return g+dist[isk]>a.g+dist[a.isk];
}
};
void insect(int x,int y,double z){next[cnt]=point[x];point[x]=cnt;v[cnt]=y;c[cnt]=z;cnt++;}
void insect2(int x,int y,double z){next2[cnt2]=point2[x];point2[x]=cnt2;v2[cnt2]=y;c2[cnt2]=z;cnt2++;}
void SPFA()
{
p[N]=1; queue<int>Q;Q.push(N); dist[N]=0.0; int i,now;
while (!Q.empty())
{
now=Q.front(); Q.pop(); p[now]=0;
for (i=point2[now];i!=-1;i=next2[i])
if (dist[v2[i]]>dist[now]+c2[i])
{dist[v2[i]]=dist[now]+c2[i];if (p[v2[i]]==0){p[v2[i]]=1;Q.push(v2[i]);}}
}
}
void Astar()
{
priority_queue<datatype>q; cnt=0;
int i,now; double nowg; datatype A; A.g=0.0; A.isk=1; q.push(A);
while ((E>0)&&(!q.empty()))
{
now=q.top().isk; nowg=q.top().g; q.pop();
if (now==N)
{E=E-nowg; if (E>=0) {cnt++;continue;} else break;}
for (i=point[now];i!=-1;i=next[i])
{A.isk=v[i];A.g=nowg+c[i];q.push(A);}
}printf("%d\n",cnt);
}
int main()
{
memset(point,-1,sizeof(point));
memset(point2,-1,sizeof(point2));
memset(next,-1,sizeof(next));
memset(next2,-1,sizeof(next2));
memset(v,0,sizeof(v));
memset(v2,0,sizeof(v2));
memset(c,0,sizeof(c));
memset(c2,0,sizeof(c2));
memset(dist,0x43,sizeof(dist));
memset(p,0,sizeof(p));
scanf("%d %d %lf\n",&N,&M,&E);
int i,x,y; double z;
for (i=1;i<=M;++i)
{
scanf("%d %d %lf\n",&x,&y,&z);
insect(x,y,z);
insect2(y,x,z);
}SPFA();
Astar();
return 0;
}
2016-03-14:改了改之前的手写堆,在codevs上比stl快多了,,,charge呵呵一笑:stl垃圾(内存我懒得改了,,,codevs上A了就行)
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int duin,N,M,next[200010],point[5010],v[200010],next2[200010],point2[5010],v2[200010];
double c[200010],c2[200010],dist[5010],duif[40000010],E,sum=0;
int cnt=0,cnt2=0,duiu[40000010];
bool vis[5010];
void insect(int x,int y,double z){cnt++;next[cnt]=point[x];point[x]=cnt;v[cnt]=y;c[cnt]=z;}
void insect2(int x,int y,double z){cnt2++;next2[cnt2]=point2[x];point2[x]=cnt2;v2[cnt2]=y;c2[cnt2]=z;}
void swapp(double &x,double &y){double z=x;x=y;y=z;}
void swapp2(int &x,int &y){int z=x;x=y;y=z;}
void SPFA()
{
int i,u; queue<int>q;
for (i=1;i<N;++i) dist[i]=1e150; memset(vis,0,sizeof(vis));
dist[N]=0; vis[N]=1; q.push(N);
while (!q.empty())
{
u=q.front(); q.pop();
for (i=point2[u];i!=0;i=next2[i])
if (dist[v2[i]]>dist[u]+c2[i])
{
dist[v2[i]]=dist[u]+c2[i];
if (vis[v2[i]]==0) {vis[v2[i]]=1;q.push(v2[i]);}
}
}
}
void weihu()
{
duiu[1]=duiu[duin];duif[1]=duif[duin];duin--;
int i=1;
while (i<duin)
{
if (((i*2+1>duin)||(duif[i]<duif[i*2+1]))&&((i*2>duin)||(duif[i]<duif[i*2])))
return;
if (i*2+1<=duin)
if (duif[i*2+1]<duif[i*2])
{
swapp(duif[i*2+1],duif[i]);
swapp2(duiu[i*2+1],duiu[i]); i=i*2+1;
}
else
{
swapp(duif[i*2],duif[i]);
swapp2(duiu[i*2],duiu[i]); i=i*2;
}
else
{
swapp(duif[i*2],duif[i]);
swapp2(duiu[i*2],duiu[i]); i=i*2;
}
}
}
void indui(int u,double f)
{
duin++;duiu[duin]=u;duif[duin]=f;
int i=duin;
while (i>1)
{
if (duif[i]>duif[i/2]) return;
else {swapp(duif[i],duif[i/2]);swapp2(duiu[i],duiu[i/2]); i=i/2;}
}
}
void astar()
{
int u,i,u2; double f,g,f2,g2; cnt=0;
duin=1; duiu[1]=1; duif[1]=dist[1];
while (duin>0)
{
u=duiu[1];f=duif[1];g=f-dist[u]; weihu();
if (u==N){sum+=f; if (sum<=E) {cnt++;continue;} else return;}
for (i=point[u];i!=0;i=next[i])
{
u2=v[i]; g2=g+c[i]; f2=g2+dist[u2];
indui(u2,f2);
}
}
}
int main()
{
scanf("%d %d %lf\n",&N,&M,&E);
int i,x,y;double z;
for (i=1;i<=M;++i)
{
scanf("%d %d %lf\n",&x,&y,&z);
insect(x,y,z); insect2(y,x,z);
}
SPFA();
astar();
printf("%d\n",cnt);
return 0;
}
这样就行啦
codevs 1835 魔法猪学院 A*寻k短路做了一个月卡死在spfa那了/(ㄒoㄒ)/~~的更多相关文章
- BZOJ 1975 魔法猪学院(A*求K短路)
显然每次贪心的走最少消耗的路径即可.那么也就是找出最短路,次短路,,,K短路之后消耗E的能量的最多的路径条数. 也就是裸的A*算法. #include <bits/stdc++.h> us ...
- k短路模板(洛谷P2483 [SDOI2010]魔法猪学院)(k短路,最短路,左偏树,priority_queue)
你谷数据够强了,以前的A*应该差不多死掉了. 所以,小伙伴们快来一起把YL顶上去把!戳这里! 俞鼎力的课件 需要掌握的内容: Dijkstra构建最短路径树. 可持久化堆(使用左偏树,因其有二叉树结构 ...
- BZOJ1975 SDOI2010魔法猪学院(启发式搜索+最短路+堆)
对反图跑最短路求出每个点到终点的最短路径,令其为估价函数大力A*,第k次到达某个点即是找到了到达该点的非严格第k短路,因为估价函数总是不大于实际值.bzoj可能需要手写堆.正解是可持久化可并堆,至今是 ...
- BZOJ-1975 魔法猪学院 K短路 (A*+SPFA)
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 1323 Solved: 433 [Submit][Statu ...
- Bzoj 1975: [Sdoi2010]魔法猪学院 dijkstra,堆,A*,K短路
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 1357 Solved: 446[Submit][Statu ...
- K短路 (A*算法) [Usaco2008 Mar]牛跑步&[Sdoi2010]魔法猪学院
A*属于搜索的一种,启发式搜索,即:每次搜索时加一个估价函数 这个算法可以用来解决K短路问题,常用的估价函数是:已经走过的距离+期望上最短的距离 通常和Dijkstra一起解决K短路 BZOJ1598 ...
- [BZOJ1975][SDOI2010]魔法猪学院(k短路,A*)
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 2748 Solved: 883[Submit][Statu ...
- bzoj 1975: [Sdoi2010]魔法猪学院 [k短路]
1975: [Sdoi2010]魔法猪学院 裸题... 被double坑死了 #include <iostream> #include <cstdio> #include &l ...
- bzoj1975: [Sdoi2010]魔法猪学院【k短路&A*算法】
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 2446 Solved: 770[Submit][Statu ...
随机推荐
- 【MVC 4】3.MVC 基本工具(创建示例项目、使用 Ninject)
作者:[美]Adam Freeman 来源:<精通ASP.NET MVC 4> 本次将考察三类工具,它们是每一位 MVC 程序员工具库的成员:DI容器.单元测试框架和模仿工具. ...
- 合工大OJ 1337 一加二减三
Description 题目描述:给一个串,形如一+二-三,求值 Input 第一行为一个正整数T,表示数据的组数,接下来有T行每行都是一个形如一+二-三的串,一,二,三均为正整数 Output 对于 ...
- jQuery对象与dom对象的区别与相互转换
什么是jQuery对象? 就是通过jQuery包装DOM对象后产生的对象.jQuery对象是jQuery独有的,其可以使用jQuery里的方法.例如: $("#test").htm ...
- UESTC 1227 & POJ 3667 Hotel
非常细腻的线段树题目啊,后来还是有个细节写错了,查了一个晚上..就不分析了. 代码: #include <iostream> #include <cstdio> #includ ...
- RMQ之ST算法模板
#include<stdio.h> #include<string.h> #include<iostream> using namespace std; ; ],M ...
- python 维吉尼亚
加密key='COMPUTER' plaintext='BLOCKCIPHERDESIGNPRINCIPLE' ascii='abcdefghijklmnopqrstuvwxyz'.upper() k ...
- 细说git merge & git rebase
git merge和git rebase两个都是用来合并两个分支用的,在使用过程中,这两个概念容易混淆. 在此,对这两个git技巧的用法进行详细描述,希望能帮助一些热爱git的朋友. -------- ...
- typicalapp.js
/** * 1.找出数字数组中最大的元素(使用Math.max函数) 2.转化一个数字数组为function数组(每个function都弹出相应的数字) 3.给object数组进行排序(排序条件是每个 ...
- final-----finalize----finally---区别
一.性质不同 (1)final为关键字: (2)finalize()为方法: (3)finally为为区块标志,用于try语句中: 二.作用 (1)final为用于标识常量的关键字,final标识的关 ...
- Node.js之事件events
Events a.EventEmitter支持多个事件监听,最大为10,也可以自定义最大数 //添加监听var EventEmitter = require('events').EventEmitte ...