BZOJ 2330 [SCOI2011]糖果 - 差分約束

發現篇博文講差分約束講的挺好
http://www.cppblog.com/menjitianya/archive/2015/11/19/212292.html

照着文章自學了一遍,發現這道題就是一道裸題嘛。
設一個超級源,表示基礎的0,要求所有的值均嚴格大於0,spfa一下,將距離數組相加就是答案,若存在負環,則不存在一個合理解。

關於差分約束的建圖,還想mark一下:
1.若求一個變量xi 滿足xxians 中ans的最小值,即xi 最小比x 大多少,變形爲xans+xi ,類似於單源最短路徑的更新,則將有向圖的邊從x連向xi ,求最短路;
2.若不等式形式爲xxians ,求ans最大值,則轉化爲最長路。
對於這道題的一個毒瘤數據:一條1->n的鏈,0向其他節點建邊的順序反向即可,從1更新到n,一次更新完畢,避免最壞退化成O(n2)

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<algorithm>

using namespace std;

const int maxn=100005;

struct edge
{
    int to,next,val;
}e[maxn<<2];

int n,m;
int head[maxn],dist[maxn],tim[maxn];
bool inq[maxn];

void insert(int a,int b,int c)
{
    static int cnt=0;
    e[++cnt].to=b;e[cnt].val=c;e[cnt].next=head[a];head[a]=cnt;
}
bool spfa()
{
    queue<int>q;
    q.push(0);inq[0]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        inq[u]=false;
        for(int i=head[u];i;i=e[i].next)
        {
            int v=e[i].to;
            if(dist[u]+e[i].val>dist[v])
            {
                dist[v]=dist[u]+e[i].val;
                tim[v]++;
                if(tim[v]>n)return false;
                if(inq[v])continue;
                inq[v]=true,q.push(v);
            }
        }
    }
    return true;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=n;i;i--)//倒着加邊 
        insert(0,i,1);
    for(int i=1,op,d1,d2;i<=m;i++)
    {
        scanf("%d%d%d",&op,&d1,&d2);
        if(op%2==0&&d1==d2){puts("-1");return 0;}
        if(op==1)insert(d1,d2,0),insert(d2,d1,0);
        else if(op==2)insert(d1,d2,1);
        else if(op==3)insert(d2,d1,0);
        else if(op==4)insert(d2,d1,1);
        else insert(d1,d2,0);
    }
    if(!spfa()){puts("-1");return 0;}
    long long ans=0;
    for(int i=1;i<=n;i++)
        ans+=dist[i];
    printf("%lld",ans);
    return 0;
} 
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章