差分约束 链式前向星 (学习

    在学习差分约束的时候先学了下链式前向星存图,因为这存的方便,也优化了最短路算法,具体代码如下:

struct xxx
{
    int v,w;
    int next;
}edge[maxn];
int head[maxn];
int cnt=0;

void add(int u,int v,int w)
{
    edge[cnt].v=v;
    edge[cnt].w=w;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}

    cnt 相当于每条路存的位置

    head[u] 表示以u为开头的最后一条路的存储位置

    .next 存的是上一条同样以u开头的路存的位置。

比如有两条路 :

                        1->2==3   

                        1->3==4 

    那么edge[1].next=head[1]=0,head[1]=1。这样存可以轻松遍历出以任意u开头的所有路径:

/******head内初值为-1******/

    int t=u;
    for(int i=head[t]; i!=-1; i=edge[i].next)
    {
        
    }
    

    差分约束其实就是把有关系的不等式转换成最短路来求 

比如: A-B<=10

            B-C<=10

            求A-C最大为多少?

对于这种情况 我们可以将两个不等式转换为路径 即:B->A==10 , C->B==10 ,然后求出C->A的最短路径即为结果

再比如: A-B>=10

               B-C>=10

               求A-C最小为多少

那我们就可以建立:B->A==10 , C->B==10 ,但此时要求C->A 的最长路

可以得出个规律:求最大值, 不等式为 <= ,求最短路

                             求最小值,不等式为 >= ,   求最长路

若是不等式同时出现 <=和>=那我们只需按所需求的最大或最小把不等式方向改变再建立路径即可;

但差分约束一般麻烦在隐含条件,看例题 POJ - 1716

题目意思就是:给出n个区间,现在要你找出一个点集,使得这n个区间都至少有2个元素在这个点集里面,问这个点集最少有几个点。

第一行一个整数n。 接下来n行,每行两个整数a,b,表示区间的左端点和右端点,被空格隔开。

题目就这么多 ,那我们首先要发现的是 区间[a,b]内的点必定超过两个 那我们以区间内点的数量 D[ ]来建立不等式就是:

D[b+1]-D[a]>=2,但这一个不等式又显然不能得到最终结果,所以我们还要加上 0<=D[i+1]-D[i]<=1,再根据题意找最少的 那就是建路求最长路, 所以说差分约束烦就在要找隐含条件。附上代码:
 

#include<stdio.h>
#include<string>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
#define maxn  10005
#define LL long long
int inf=1000000007;

struct xxx
{
    int v,w;
    int next;
}edge[4*maxn];
int head[maxn];
int dis[maxn];
bool vis[maxn];
int cnt=0;
int mi=maxn,ma;

void add(int u,int v,int w)
{
    edge[cnt].v=v;
    edge[cnt].w=w;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}
void spfa()
{
    dis[mi]=0;
    queue<int> q;
    q.push(mi);
    vis[mi]=1;
    while(!q.empty())
    {
        int t=q.front();
        q.pop();
        vis[t]=0;
        for(int i=head[t];i!=-1;i=edge[i].next)
        {
            int j=edge[i].v;
            if(dis[j]<dis[t]+edge[i].w)
            {
                dis[j]=dis[t]+edge[i].w;
                if(!vis[j])
                {
                    q.push(j);
                    vis[j]=1;
                }
            }
        }
    }
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        mi=maxn;
        ma=0;
        cnt=0;
        memset(head,-1,sizeof(head));
        memset(edge,0,sizeof(edge));
        for(int i=0;i<n;i++)
        {
            int a,b;
            scanf("%d %d",&a,&b);
            mi=min(a,mi);
            ma=max(b+1,ma);
            add(a,b+1,2);
        }
        for(int i=mi;i<=ma;i++){
            vis[i]=0;
            dis[i]=-inf;
        }
        for(int i=mi;i<ma;i++)
        {
            add(i,i+1,0);
            add(i+1,i,-1);
        }
        spfa();
        printf("%d\n",dis[ma]);
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值