最短路随笔

刷到的时候再继续更新 QAQ

反向建图&&Dij
emmm 正向建图TLE了 ，就是给了我们n个点m条单向边和终点S
问我们最快到达终点的时间是多久 或者无法到达 此外还有q个起始点
（本来我是建了一个超级原点的 结果WA了） emm 后来想着还是反向建图吧，枚举起点跑一跑取最小 就过了 挺简单的

 #include<bits/stdc++.h>
#include<stdlib.h>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<queue>
#include<time.h>
#include <cstdio>
#include <iostream>
#include <vector>
#define ll int
//#define int long long
#define inf 0x3f3f3f3f
#define mods 1000000007
#define modd 998244353
#define PI acos(-1)
#define fi first
#define se second
#define lowbit(x) (x&(-x))
#define mp make_pair
#define pb push_back
#define si size()
#define E exp(1.0)
#define fixed cout.setf(ios::fixed)
#define fixeds(x) setprecision(x)
#define IOS ios::sync_with_stdio(false);cin.tie(0)
 using namespace std;
 ll gcd(ll a,ll b){if(a<0)a=-a;if(b<0)b=-b;return b==0?a:gcd(b,a%b);}
template<typename T>void read(T &res){bool flag=false;char ch;while(!isdigit(ch=getchar()))(ch=='-')&&(flag=true);
for(res=ch-48;isdigit(ch=getchar());res=(res<<1)+(res<<3)+ch - 48);flag&&(res=-res);}

const int manx=1e5+5;
const int mamx=3*1e5+5;
ll head[manx],d[manx];
bool vis[manx];
ll n,m,k=0,s,kk,ans;
struct node{
    ll v,next,w;
}a[mamx];
void add(ll u,ll v,ll w)
{
    a[++k].next=head[u];
    a[k].w=w;
    a[k].v=v;
    head[u]=k;
}
void dij()
{
    memset(d,inf,sizeof(d));
    memset(vis,0,sizeof(vis));
    d[s]=0;
    priority_queue<pair<ll,ll> >q;
    q.push(mp(0,s));
    while(q.size()){
        ll u=q.top().se;
        q.pop();
        if(vis[u]) continue;
        vis[u]=1;
        for(int i=head[u];i;i=a[i].next){
            ll v=a[i].v,w=a[i].w;
           if(d[v]>d[u]+w){
          //printf("qwq %lld\n",d[v]);      //松弛操作，更新距离
            d[v]=d[u]+w;
           // printf("sad %lld %lld %lld %lld\n",d[u],u,res,d[v]);
            q.push(mp(-d[v],v)); //把更新的距离和点入队，这里距离取负变成小根堆
            }

        }
    }
}
ll cnt[1200];
ll e[1010][1010];
signed main(){
ll ed;
while(scanf("%d%d%d",&n,&m,&ed)!=EOF){
 if(n<=0||m<=0||ed<=0){
    break;
 }

 k=0;
 memset(head,0,sizeof(head));
memset(e,inf,sizeof(e));
for(int i=1;i<=m;i++){
ll u,v,w;
read(u);
read(v);
read(w);
e[u][v]=min(e[u][v],w);
//printf("Q  %lld\n",e[u][v]);
}
 ll q;
 ll op=0;
 read(q);
 while(q--){
 ll x;
 read(x);
 cnt[++op]=x;
 }

 for(int i=1;i<=n;i++){
    for(int j=1;j<=n;j++){
        if(i==j||e[i][j]==inf){continue;}

         add(j,i,e[i][j]);
    }
 }
s=ed;
dij();
ll mi=inf;
for(int i=1;i<=op;i++){
 mi=min(d[cnt[i]],mi);

}
if(mi>=inf){
    printf("-1\n");
    continue;
}
else{
    printf("%d\n",mi);
    continue;
}

}
return 0;


}

放一个小插曲吧 数论 拓展欧几里得判断是否解存在
牛客&&数学
给出a k b c 问是否存在 a*n 在k进制下 数位和 模b等于c
等价于 ax+by=c 有解则c%gcd(a,b)=0 而且由于进位影响 因此判定
ax+by+z(1-k)=c 是否有解
判断c%gcd(b,gcd(a,1-k)) 是否为0