C++博客-寒月-随笔分类-图论

最大流模板【sap+gap优化+当前弧优化, 采用邻接表+反向弧指针,O(mn^2)】

wuxu — Tue, 10 May 2011 05:17:00 GMT

摘要: /**//* sap+gap优化+当前弧优化, 采用邻接表+反向弧指针时间复杂度:O(mn^2)*/#includeusing namespace std;const int maxnode = 1024... 阅读全文

wuxu 2011-05-10 13:17 发表评论

最大流模板【sap+gap优化,邻接矩正,O(mn^2)】

wuxu — Tue, 10 May 2011 05:12:00 GMT

摘要: /**//* sap+gap优化,采用邻接矩正时间复杂度:O(mn^2)*/#include#includeusing namespace std;const int msize=205; &nbs... 阅读全文

wuxu 2011-05-10 13:12 发表评论

最大流模板【EdmondsKarp算法，O(m^2n)】

wuxu — Tue, 10 May 2011 05:10:00 GMT

/*
   EdmondsKarp算法，O(m^2n)
*/
#include<iostream>
#include<queue>
using namespace std;
const int maxn=205;
const int inf=0x7fffffff;

int r[maxn][maxn]; //残留网络，初始化为原图
bool visit[maxn];
int pre[maxn];
int m,n;

bool bfs(int s,int t)  //寻找一条从s到t的增广路，若找到返回true
{
    int p;
    queue<int > q;
    memset(pre,-1,sizeof(pre));
    memset(visit,false,sizeof(visit));

    pre[s]=s;
    visit[s]=true;
    q.push(s);
    while(!q.empty())
    {
        p=q.front();
        q.pop();
        for(int i=1;i<=n;i++)
        {
            if(r[p][i]>0&&!visit[i])
            {
                pre[i]=p;
                visit[i]=true;
                if(i==t) return true;
                q.push(i);
            }
        }
    }
    return false;
}

int EdmondsKarp(int s,int t)
{
   int flow=0,d,i;
   while(bfs(s,t))
   {
       d=inf;
       for(i=t;i!=s;i=pre[i])
           d=d<r[pre[i]][i]? d:r[pre[i]][i];
       for(i=t;i!=s;i=pre[i])
       {
           r[pre[i]][i]-=d;
           r[i][pre[i]]+=d;
       }
       flow+=d;
   }
   return flow;
}

int main()
{
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        int u,v,w;
        memset(r,0,sizeof(r));///
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%d",&u,&v,&w);
            r[u][v]+=w;
        }
        printf("%d\n",EdmondsKarp(1,n));
    }
    return 0;
}

wuxu 2011-05-10 13:10 发表评论

最大流模板【dinic算法，边表存储（仿指针），实用于稀疏图，O(mn^2)】

wuxu — Tue, 10 May 2011 05:07:00 GMT

摘要: /**//* dinic算法，边表存储（仿指针），实用于稀疏图，O(mn^2)*/#includeusing namespace std;const int maxn=205; //最大定点数const int maxm=205... 阅读全文

wuxu 2011-05-10 13:07 发表评论

最大流模板【dinic算法,邻接矩正,O(mn^2)】

wuxu — Tue, 10 May 2011 05:05:00 GMT

摘要: /**//* dinic算法,邻接矩正,O(mn^2)*/#include#include#includeusing namespace std;const int maxnode=205;const int&n... 阅读全文

wuxu 2011-05-10 13:05 发表评论

hdu3560

wuxu — Wed, 25 Aug 2010 10:55:00 GMT

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3560
题意：求连通分量的个数和环图的个数。如果是环图那么每个点的度数应该为2。

#include<iostream>
using namespace std;
typedef struct
{
    int v;
    int next;
}edge;
edge e[600005];
int head[100005];
int deg[100005];
bool vis[100005];
int n,m,num1,num2;
void bfs(int u)
{
     num1++;
     vis[u]=1;
     bool flag=0;
     int front=0,rear=0;
     int q[100005];
     q[rear++]=u;
     if(deg[u]!=2) flag=1;
     while(front<rear)
     {
         int x=q[front++];
         for(int i=head[x];i!=-1;i=e[i].next)
         {
             if(!vis[e[i].v])
             {
                 vis[e[i].v]=1;
                 if(deg[e[i].v]!=2) flag=1;
                 q[rear++]=e[i].v;
             }
         }
     }
     if(!flag) num2++;
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF&&(n||m))
    {
        int i,j,uu,vv,tot=0;
        memset(head,-1,sizeof(head));
        memset(deg,0,sizeof(deg));
        memset(vis,0,sizeof(vis));
        for(i=0;i<m;i++)
        {
            scanf("%d%d",&uu,&vv);
            deg[uu]++;
            deg[vv]++;

            e[tot].v=vv;
            e[tot].next=head[uu];
            head[uu]=tot++;

            e[tot].v=uu;
            e[tot].next=head[vv];
            head[vv]=tot++;
        }

        num1=0,num2=0;
        for(i=0;i<n;i++)
        if(!vis[i])
        bfs(i);

        printf("%d %d\n",num1,num2);
    }
    //system("pause");
    return 0;
}

wuxu 2010-08-25 18:55 发表评论

pku1178

wuxu — Fri, 20 Aug 2010 08:34:00 GMT

floyd+枚举. 题目来源:http://acm.pku.edu.cn/JudgeOnline/problem?id=1178
我们将每一个格子看作一个结点,将该结点与附近的骑士能一步到达的结点连一条边.然后用floyd求出每对结点间的最段路径,这就是骑士从一个格子走到另外一个格子的最小步数.对于国王我们可以直接用坐标来求出他通往任意两个格子间的最小步数. 然后枚举他们的汇聚点和国王与骑士的汇合点,再枚举是哪个骑士和国王回合.取最小的那个值.详细见代码:

#include<iostream>
#include<cmath>
using namespace std;
const int inf=100000000;
int location[64],knight[64][64];
int kn_dir[8][2]={-2,-1,-2,1,-1,-2,-1,2,1,-2,1,2,2,-1,2,1};
char str[130];
void floyd()
{
     int i,j,k;
     for(k=0;k<64;k++)
     for(i=0;i<64;i++)
     for(j=0;j<64;j++)
     {
         if(knight[i][j]>knight[i][k]+knight[k][j])
         knight[i][j]=knight[i][k]+knight[k][j];
     }
}
void Init()
{
     int i,j;
     for(i=0;i<64;i++)
     for(j=0;j<64;j++)
     {
         if(i==j) knight[i][j]=0;
         else knight[i][j]=inf;
     }
     for(i=0;i<64;i++)
     {
         for(j=0;j<8;j++)
         {
             int x=i/8+kn_dir[j][0];
             int y=i%8+kn_dir[j][1];
             if(x>=0&&x<8&&y>=0&&y<8)
             {
                  knight[i][x*8+y]=1;
             }
         }
     }
     floyd();
}
int main()
{
    int i,j,k,num,ans=inf;
    Init();
    scanf("%s",str);
    for(i=0,num=0;str[i];i+=2)
    {
         location[num++]=8*(str[i+1]-'1')+(str[i]-'A');
    }
    for(i=0;i<64;i++)
    for(j=0;j<64;j++)
    {
         int xx=abs(location[0]/8-j/8);
         int yy=abs(location[0]%8-j%8);
         int king_len=xx>yy ? xx:yy;
         int knight_len=inf;
         for(k=1;k<num;k++)
         {
              int temp=0;
              for(int v=1;v<num;v++)
              if(v!=k)
              {
                   temp+=knight[location[v]][i];
              }
              temp+=knight[location[k]][j]+knight[j][i];
              if(knight_len>temp) knight_len=temp;
         }
         if(knight_len==inf&&king_len<ans) ans=king_len;
         else if(king_len+knight_len<ans) ans=king_len+knight_len;
    }
    printf("%d\n",ans);
    //system("pause");
    return 0;
}

wuxu 2010-08-20 16:34 发表评论

pku1062

wuxu — Wed, 11 Aug 2010 08:40:00 GMT

枚举+最短路
枚举最低等级,保证酋长在交易范围,然后求最短路。

#include<iostream>
using namespace std;

const int inf=100000000;
int m,n;
int map[105][105];
int dis[105],vis[105],lev[105];

int dijkstra()
{
    int ans=inf;
    int i,j,k;
    for(i=lev[1]-m;i<=lev[1];i++)
    {
        memset(vis,0,sizeof(vis));
        for(j=1;j<=n+1;j++)
            dis[j]=map[1][j];
        vis[1]=1;

           for(j=1;j<=n;j++)
        {
            int min=inf,cur=-1;
            for(k=1;k<=n+1;k++)
            {
                if(!vis[k]&&lev[k]>=i&&lev[k]<=i+m&&dis[k]<min)
                {
                    min=dis[k];
                    cur=k;
                }
            }
            if(cur==-1) break;
            vis[cur]=1;

            for(k=1;k<=n+1;k++)
            {
                if(!vis[k]&&lev[k]>=i&&lev[k]<=i+m&&map[cur][k]!=inf&&dis[k]>dis[cur]+map[cur][k])
                    dis[k]=dis[cur]+map[cur][k];
            }
        }
        if(dis[n+1]<ans)
            ans=dis[n+1];
    }
    return ans;
}

int main()
{
    int p,l,x,t,v;
    int i,j;
    scanf("%d%d",&m,&n);
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
        {
            if(i==j) map[i][j]=0;
            else map[i][j]=inf;
        }
    for(i=1;i<=n;i++)
    {
        scanf("%d%d%d",&p,&l,&x);
        map[i][n+1]=p;
        lev[i]=l;
        for(j=1;j<=x;j++)
        {
            scanf("%d%d",&t,&v);
            map[i][t]=v;
        }
    }
    lev[n+1]=lev[1];
    printf("%d\n",dijkstra());
    return 0;
}

wuxu 2010-08-11 16:40 发表评论