hdu2544 最短路 Bellman_Ford Dijkstra Folyd

Bellman_Ford基本方法,和使用邻接表和队列的优化。由于没有负环而且最短路一定存在,所以不需要判负环那一步。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<string>
#include<map>
#include<set>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<sstream>
#define LL long long
#define OJ_DEBUG 0
#define READ_FILE 0
using namespace std;
const int NN_MAX = 110;
const int MM_MAX = NN_MAX*NN_MAX;
const int INF = 0x1fffffff;
struct Edge{
  int a,b,c;
  Edge (int a1=0,int b1=0,int c1=0):a(a1),b(b1),c(c1){}
  bool operator < (const Edge& chs) const {return c<chs.c;}
}theEdge[MM_MAX];
/**********************************************************/
int n,m,s,t,lenEdge;
vector<int>G[NN_MAX];
int d[NN_MAX];
/**********************************************************/
int min_2 (int x,int y) {return x<y?x:y;}
int max_2 (int x,int y) {return x>y?x:y;}
void Bellman_Ford_1();
void Bellman_Ford_2();
/**********************************************************/
int main()
{
    if (READ_FILE) freopen ("in.txt","r",stdin);
    while( scanf("%d %d",&n,&m) && (n||m) )
    {
        lenEdge=0;
        for(int i=0;i<=n;i++) G[i].clear();
        int t1,t2,t3;
        for(int i=0;i<m;i++){
            scanf("%d %d %d",&t1,&t2,&t3);
            theEdge[lenEdge++]=Edge(t1,t2,t3); G[t1].push_back(lenEdge-1);
            theEdge[lenEdge++]=Edge(t2,t1,t3); G[t2].push_back(lenEdge-1);
        }
        s=1,t=n;
        //Bellman_Ford_1();
        Bellman_Ford_2();
        printf("%d\n",d[t]);
    }
    return 0;
}
void Bellman_Ford_1()
{
    for(int i=0;i<=n;i++) d[i]=INF;
    d[s]=0;
    for(int k=0;k<n-1;k++){
        bool finish=true;
        for(int i=0;i<lenEdge;i++){
            int x=theEdge[i].a, y=theEdge[i].b;
            if(d[x]<INF) d[y]=min_2(d[y],d[x]+theEdge[i].c);
            finish=false;
        }
        if(finish) break;
    }
}
void Bellman_Ford_2()
{
    for(int i=0;i<=n;i++) d[i]=INF;
    d[s]=0;
    queue<int>qee;
    qee.push(s);
    while(!qee.empty())
    {
        int x=qee.front();qee.pop();
        for(int i=0;i<G[x].size();i++){
            Edge e=theEdge[G[x][i]];
            if(d[x]<INF && d[e.b]>d[x]+e.c){
                d[e.b]=d[x]+e.c;
                qee.push(e.b);
            }
        }
    }
}

Dijkstra基本方法,和使用邻接表和优先队列的优化。注意优先队列元素需要重载<

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<string>
#include<map>
#include<set>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<sstream>
#define LL long long
#define OJ_DEBUG 0
#define READ_FILE 1
using namespace std;
const int NN_MAX = 110;
const int MM_MAX = NN_MAX*NN_MAX;
const int INF = 0x1fffffff;
struct Edge{
  int a,b,c;
  Edge (int a1=0,int b1=0,int c1=0):a(a1),b(b1),c(c1){}
  bool operator < (const Edge& chs) const {return c<chs.c;}
}theEdge[MM_MAX];
struct CostAndNode{
    int d,b;
    CostAndNode(int d1=0,int b1=0):d(d1),b(b1){}
    bool operator < (const CostAndNode& chs) const {return d>chs.d;}
};
/**********************************************************/
int n,m,s,t,lenEdge;
vector<int>G[NN_MAX];
int d[NN_MAX],vis[NN_MAX],maps[NN_MAX][NN_MAX];
/**********************************************************/
int min_2 (int x,int y) {return x<y?x:y;}
int max_2 (int x,int y) {return x>y?x:y;}
void Dijkstra_1();
void Dijkstra_2();
/**********************************************************/
int main()
{
    if (READ_FILE) freopen ("in.txt","r",stdin);
    while( scanf("%d %d",&n,&m) && (n||m) )
    {
        lenEdge=0;
        for(int i=0;i<=n;i++) G[i].clear();
        int t1,t2,t3;
        /*for(int i=0;i<=n;i++)
            for(int j=i;j<=n;j++)
                maps[i][j]=maps[j][i]=(i==j?0:INF);*/
        for(int i=0;i<m;i++){
            scanf("%d %d %d",&t1,&t2,&t3);
            maps[t1][t2]=t3;
            maps[t2][t1]=t3;
            theEdge[lenEdge++]=Edge(t1,t2,t3); G[t1].push_back(lenEdge-1);
            theEdge[lenEdge++]=Edge(t2,t1,t3); G[t2].push_back(lenEdge-1);
        }
        s=1,t=n;
        //Dijkstra_1();
        Dijkstra_2();
        printf("%d\n",d[t]);
    }
    return 0;
}
void Dijkstra_1()
{
    memset(vis,0,sizeof(vis));
    for(int i=0;i<=n;i++) d[i]=INF;
    d[s]=0;
    for(int k=0;k<n;k++){
        int x=0;
        for(int i=1;i<=n;i++)
            if(!vis[i] && d[i]<d[x])
                x=i;
        vis[x]=1;
        for(int y=1;y<=n;y++)
            if(!vis[y])
                d[y]=min_2(d[y],d[x]+maps[x][y]);
    }
}
void Dijkstra_2()
{
    memset(vis,0,sizeof(vis));
    for(int i=0;i<=n;i++) d[i]=INF;
    d[s]=0;
    priority_queue<CostAndNode>pq;
    pq.push(CostAndNode(0,s));
    while(!pq.empty())
    {
        CostAndNode cad=pq.top();pq.pop();
        int x=cad.b;
        if(vis[x]) continue;
        vis[x]=1;
        for(int i=0;i<G[x].size();i++){
            Edge e=theEdge[G[x][i]];
            if(!vis[e.b] && d[e.b]>d[x]+e.c){
                d[e.b]=d[x]+e.c;
                pq.push( CostAndNode(d[e.b],e.b) );
            }
        }
    }
}
floyd也能过
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<string>
#include<map>
#include<set>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<sstream>
#define LL long long
#define OJ_DEBUG 0
#define READ_FILE 0
using namespace std;
const int NN_MAX = 110;
const int MM_MAX = NN_MAX*NN_MAX;
const int INF = 0x0fffffff;
/**********************************************************/
int n,m,s;
int maps[NN_MAX][NN_MAX];
/**********************************************************/
int min_2 (int x,int y) {return x<y?x:y;}
int max_2 (int x,int y) {return x>y?x:y;}
void floyd();
/**********************************************************/
int main()
{
    if (READ_FILE) freopen ("in.txt","r",stdin);
    while( scanf("%d %d",&n,&m) && (n||m) )
    {
        int t1,t2,t3;
        for(int i=0;i<=n;i++)
            for(int j=i;j<=n;j++)
                maps[i][j]=maps[j][i]=(i==j?0:INF);
        for(int i=0;i<m;i++){
            scanf("%d %d %d",&t1,&t2,&t3);
            maps[t1][t2]=t3;
            maps[t2][t1]=t3;
        }
        floyd();
        printf("%d\n",maps[1][n]);
    }
    return 0;
}
void floyd()
{
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(maps[i][k]<INF && maps[k][j]<INF)
                    maps[i][j]=min_2(maps[i][j],maps[i][k]+maps[k][j]);
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,关于 HDU4992 求所有原根的问题,我可以给出以下解答: 1. 首先,我们需要知道什么是原根。原根是指模 n 意义下,所有与 n 互质的数 a 的最小正整数 k 次幂对 n 取模的值覆盖了所有与 n 互质的数。简单来说,就是如果一个数 a 是模 n 意义下的原根,那么 a 的任何正整数次幂对 n 取模的值都不相同,且能覆盖所有与 n 互质的数。 2. 为了求模 n 意义下的所有原根,我们需要先求出与 n 互质的数的个数 phi(n)。phi(n) 可以使用欧拉函数求出。 3. 然后,我们需要枚举模 n 意义下的所有数,判断它是否是原根。具体来说,对于每个 a,我们需要判断 a 的每个小于 phi(n) 的正整数次幂对 n 取模的值是否都不相同,且能覆盖所有与 n 互质的数。如果是,那么 a 就是模 n 意义下的原根。 4. 代码实现可以参考以下 Java 代码: ``` import java.util.*; public class Main { static int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b); } static int phi(int n) { int res = n; for (int i = 2; i * i <= n; i++) { if (n % i == 0) { res = res / i * (i - 1); while (n % i == 0) { n /= i; } } } if (n > 1) { res = res / n * (n - 1); } return res; } static int pow(int a, int b, int mod) { int res = 1; while (b > 0) { if ((b & 1) != 0) { res = res * a % mod; } a = a * a % mod; b >>= 1; } return res; } static boolean check(int a, int n, int phi) { for (int i = 1, j = pow(a, i, n); i <= phi; i++, j = j * a % n) { if (j == 1) { return false; } } return true; } public static void main(String[] args) { Scanner scanner = new Scanner(System.in); while (scanner.hasNext()) { int n = scanner.nextInt(); int phi = phi(n); List<Integer> ans = new ArrayList<>(); for (int i = 1; i < n; i++) { if (gcd(i, n) == 1 && check(i, n, phi)) { ans.add(i); } } Collections.sort(ans); for (int x : ans) { System.out.print(x + " "); } System.out.println(); } } } ``` 其中,gcd 函数用于求最大公约数,phi 函数用于求欧拉函数,pow 函数用于快速幂求模,check 函数用于判断一个数是否是原根。在主函数中,我们依次读入每个 n,求出 phi(n),然后枚举模 n 意义下的所有数,判断它是否是原根,将所有原根存入一个 List 中,最后排序输出即可。 希望我的回答能够帮到你,如果你有任何问题,欢迎随时提出。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值