欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

【模板】图论

程序员文章站 2022-05-22 14:35:56
...

基础图论

链式前向星

带权值

int head[400001],ver[2000001],nxt[2000001],val[2000001],tot=0;
void add(int x,int y,int z) {
    ver[++tot]=y,val[tot]=z,nxt[tot]=head[x],head[x]=tot;
    ver[++tot]=x,val[tot]=z,nxt[tot]=head[x],head[x]=tot;
}

不带权值

int head[100001]={0},nxt[400001],ver[400001],tot=0;
void add(int x,int y) {
    ver[++tot]=y,nxt[tot]=head[x],head[x]=tot;
    ver[++tot]=x,nxt[tot]=head[y],head[y]=tot;
}

最短路算法

SPFA

int d[100001];
bool iq[100001] = {0};
queue<int>q;
void spfa(int s) {
    memset(d, 0x7f, sizeof(d));
    d[s] = 0;
    q.push(s);
    while(! q.empty()) {
        s = q.front(), q.pop();
        iq[s] = 0;
        for(int i = head[s]; i; i = nxt[i])
            if(d[ver[i]] > d[s] + val[i]){
                d[ver[i]] = d[s] + val[i];
                if(! iq[ver[i]])iq[ver[i]] = 1, q.push(ver[i]);
        }
    }
}

堆优化迪杰斯特拉

int d[100001];
bool v[100001];
priority_queue<pair<int,int> >q;
void dij(int s) {
    memset(d, 0x7f, sizeof(d));
    d[s] = 0;
    q.push(make_pair(0, s));
    while(! q.empty()) {
        s = q.top().second, q.pop();
        if(v[s])continue;
        v[s] = 1;
        for(int i = head[s]; i; i = nxt[i])
            if(d[ver[i]] > d[s] + val[i])
                q.push(make_pair(-(d[ver[i]] = d[s] + val[i]), ver[i]));
    }
}

高级图论

线段树优化连边

(非递归线段树)

int N=1;
void build(int n) {
    for(;N <= n + 1; N <<= 1);
    for(int i = N - 1; i >= 1; i --)
        add(i, i << 1 | 1, 0), add(i, i << 1, 0);
}
int find(int x) {
    return x + N;
}
void add_tree(int s, int t, int x, int v) {
    x = find(x);
  for(s = N + s - 1, t = N + t + 1; s ^ t ^ 1; s >>= 1, t >>= 1) {
      if(~ s & 1)add(x, s ^ 1, v);
      if(t & 1)add(x, t ^ 1, v);
  }
}