POJ 1860

      做惯了spfa的邻接表算法，转换成邻接矩阵就SB了，直接wa了好几天。最终在看着别人的代码敲完邻接矩阵后才发现自己的基础真的很差唉，题目做完快半个月了，今天再整理时又看到了这题，不禁愕然！唉，最短路学的还是不怎么扎实！

 

先上代码，日后慢慢看（代码中残留着邻接表的部分）。

 

#include<cstdio> #include<cstring> #include<queue> #include<algorithm> using namespace std; #define init(a,what) memset(a,what,sizeof(a)) #define read freopen("zx.in","r",stdin) #define write freopen("zx.out","w",stdout) #define INF 0x3f3f3f3f #define eps 1e-8 const int MAXN=100+10; //const int MAXM=1000+10; queue<int>q; int kmoney,mnum,haven,ans,a,b,circle[MAXN];//u[MAXM],v[MAXM],next[MAXN],first[MAXN]; double curmoney,dist[MAXN],r[MAXN][MAXN],c[MAXN][MAXN];//rab,cab,rba,cba,w[MAXN],m[MAXN] bool vis[MAXN]; /*void read_G(int a,int b,double cminte,double cexcept) { ans++; u[ans]=a, v[ans]=b, w[ans]=cminte, m[ans]=cexcept; next[ans]=first[u[ans]]; first[u[ans]]=ans; }*/ bool spfa(int start) { init(vis,false); init(circle,0); for(int i=1;i<=kmoney;i++) dist[i]=(i==start ? curmoney : 0); q.push(start); circle[start]++; while(!q.empty()) { int xx=q.front(); q.pop(); vis[xx]=false; for(int i=1;i<=kmoney;i++) { if(r[xx][i]>eps && (dist[i]+eps)<(dist[xx]-c[xx][i])*r[xx][i]) { dist[i]=(dist[xx]-c[xx][i])*r[xx][i]; if(!vis[i]) { vis[i]=true; q.push(i); } if(++circle[i]>kmoney) return true;//优化！ } } /*for(int e=first[xx];e!=-1;e=next[e]) { if(dist[v[e]]<(dist[xx]-m[xx])*w[e]) { dist[v[e]]=(dist[xx]-m[xx])*w[e]; if(!vis[v[e]]) { vis[v[e]]=true; q.push(v[e]); } } }*/ } return dist[start]>curmoney+eps; } int main() { //read, write; while(scanf("%d%d%d%lf",&kmoney,&mnum,&haven,&curmoney)!=EOF) { for(int i=1;i<=mnum;i++) { //ans=0; scanf("%d%d",&a,&b); scanf("%lf%lf%lf%lf",&r[a][b],&c[a][b],&r[b][a],&c[b][a]); //scanf("%d%d%lf%lf%lf%lf",&a,&b,&rab,cab,&rba,&cba); //read_G(a,b,rab,cab); //read_G(b,a,rba,cba); } if(spfa(haven)) printf("YES/n"); else printf("NO/n"); } return 0; }