C++博客-Smile-文章分类-多校联赛题解http://www.cppblog.com/Smile3/category/17318.htmlSmilezh-cnThu, 17 Nov 2011 03:00:21 GMTThu, 17 Nov 2011 03:00:21 GMT60BIT 第三场 The King’s Problemhttp://www.cppblog.com/Smile3/articles/151498.htmlSmile3Smile3Wed, 20 Jul 2011 10:47:00 GMThttp://www.cppblog.com/Smile3/articles/151498.htmlhttp://www.cppblog.com/Smile3/comments/151498.htmlhttp://www.cppblog.com/Smile3/articles/151498.html#Feedback0http://www.cppblog.com/Smile3/comments/commentRss/151498.htmlhttp://www.cppblog.com/Smile3/services/trackbacks/151498.html知识点：强连通分量、缩点、再加用最大二分匹配求最小路径覆盖。
bug处：突然发现没有模板自己敲的时候巨多的错误，比如说求强连通分量时连反向表都没有建，再比如说求二分匹配的时候在一个图里表示两个图的时候应是flag标记数组写不好，注意是在if(0== flag[b[i][j]])之后写flag[b[i][j]] = 1;因为此时是j点在交错路上面。最后是没有想清楚那个最小路径覆盖的数目应该是缩点之后的有向无环图的点的个数，即强连通分量的个数减去最大匹配数。
思维：觉得还是知识点懂得太少了，并且现场有点怕，不敢去啃这个题目，我觉得比较好的策略是我们把简单题目做完之后就选定一道题目，啃，我觉得自己这方面没有做好，随着我们知识系统的完善，我觉得我们应该每一次都努力的啃一道难题出来。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>

using namespace std;

#define max 5010
int n,m;
vector<int> a[max],b[max],adjact[max];
int flag[5010];
int id[max];
int result[max];
int pipe[max];

void initone()//初始化注意建反向图
{
     memset(a,0,sizeof(a));
     memset(adjact,0,sizeof(adjact));
     scanf("%d%d",&n,&m);
     for( int i=1; i<=m; i++ )
     {
          int aa,bb;
          scanf("%d%d",&aa,&bb);
          a[aa].push_back(bb);
          adjact[bb].push_back(aa);
     }
}

void vistone(int i,int &sig)//第一次深度优先搜索，注意其中的id为改点的最后离开时间，所以id的赋值应该要放到递归的后面，注意sig的&符号
{
     flag[i] = 1;
     int num = a[i].size();
     for( int j=0; j<num; ++j )
     {
          if( 0 == flag[a[i][j]] )
            vistone(a[i][j],sig);
     }
     id[++sig] = i;
     return;
}

void visttwo(int i,int sig)//注意此时用到的图为反向之后的图
{
     result[i] = sig;
     flag[i] = 1;
     int num = adjact[i].size();
     for( int j=0; j<num; ++j )
     {
          if( 0 == flag[adjact[i][j]] )
             visttwo(adjact[i][j],sig);
     }
     return;
}
                  
int f()//用来求强连通分量，首先对原图dfs一次，记住每个点的离开时间，然后按离开时间从大到小再对反向图dfs一次，注意第二次的sig要自增。
{
     memset(id,0,sizeof(id));
     memset(flag,0,sizeof(flag));
     int sig = 0;
     for( int i=1; i<=n; ++i )
     {
         if( 0 == flag[i] )
           vistone(i,sig);
     }
     sig = 0;
     memset(flag,0,sizeof(flag));
     for( int i=n; i>=1; --i )
     {
          if( 0 == flag[ id[i] ] )
          {
              visttwo(id[i],++sig);
          }
     }
     return sig;
}

void inittwo( )//建缩点之后的图，注意此时建成之后的图中的点应该是强连通分量。所以注意建图过程中的赋值。
{
     memset(b,0,sizeof(b));
     for( int i=1; i<=n; i++ )
     {
          int num = a[i].size();
          for( int j=0; j<num; ++j )
          {
               if( result[i] != result[a[i][j]] )
                  b[result[i]].push_back(result[a[i][j]]);
          }
     }
     return;
}

bool find(int i)//注意find函数中的flag赋值处，pipe数组的赋值，只要单向赋值。
{
     int num = b[i].size();
     for( int j=0; j<num; ++j )
     {
          if( 0 == flag[b[i][j]] )
          {
              flag[b[i][j]] = 1;
              if( 0 == pipe[b[i][j]] || find( pipe[b[i][j]] ) )
              {
                  pipe[b[i][j]] = i;
                  return true;
              }
          }
     }
     return  false;
}
 
int hungary(int ss)//注意每一次调用find之前都要将flag数组中的元素赋为0.
{
    int sum = 0;
    memset(pipe,0,sizeof(pipe));
    for( int i=1; i<=ss; i++ )
    {
         memset(flag,0,sizeof(flag));
         if( find(i) )
           sum++;
    }
    return sum;
}    

int main()
{
   int t;
   scanf("%d",&t);
   while( t-- )
   {
          initone();
          int ss = f();
          inittwo();
          int nn = hungary(ss);
          printf("%d\n",ss-nn);
   }
   return 0;
}

Smile3 2011-07-20 18:47 发表评论
]]>hnu第一场 Cat vs. Dog《最大二分匹配 匈牙利算法》http://www.cppblog.com/Smile3/articles/151294.htmlSmile3Smile3Mon, 18 Jul 2011 04:06:00 GMThttp://www.cppblog.com/Smile3/articles/151294.htmlhttp://www.cppblog.com/Smile3/comments/151294.htmlhttp://www.cppblog.com/Smile3/articles/151294.html#Feedback0http://www.cppblog.com/Smile3/comments/commentRss/151294.htmlhttp://www.cppblog.com/Smile3/services/trackbacks/151294.html重点：在于建图，建好图之后就一个纯裸的二分匹配。
知识：匈牙利算法 注意 最大独立集是集合的总元素个数减去最大二分匹配数
bug：此题无bug处，数据如果要符合实际情况，无法出特殊数据，要卡只能卡时间。

#include<cstdio>//此次建图是建的二分图，上面一篇文章是用的一个图求二分匹配。
#include<cstring>
#include<vector>
using namespace std;

#define max 510
int ch1[max][2];
int ch2[max][2];
int flag[max];
int vy[max];
int vx[max];
vector<int> p1[max];
int num1,num2;
int p,m,n;

void init()
{
     char a,b;
     int d,c;
     memset(p1,0,sizeof(p1));
     num1 = 0;
     num2 = 0;
     scanf("%d%d%d",&m,&n,&p);
     for( int i=1; i<=p; ++i )
     {
          getchar();
          scanf("%c%d %c%d",&a,&c,&b,&d);
          if( 'C' == a )
          {
             ch1[++num1][0] = c;
             ch1[num1][1] = d;
          }
          else
          {
              ch2[++num2][0] = c;
              ch2[num2][1] = d;
          }
     }
     for( int i=1; i<=num1; ++i )
         for( int j=1; j<=num2; ++j )
         {
              if( ch1[i][0] == ch2[j][1] || ch1[i][1] == ch2[j][0] )
              p1[i].push_back(j);
         }
         return ;
}

bool find(int i)//注意flag数组的赋值
{
     int num = p1[i].size();
     for( int j=0; j<num; ++j )
     {
          if( flag[ p1[i][j] ] == 0 )
          {
              flag[ p1[i][j] ] = 1;
              if( vy[ p1[i][j] ] == 0 || find( vy[ p1[i][j] ] ) )
              {
                  vx[i] = p1[i][j];
                  vy[ p1[i][j] ] = i;
                  return true;
              }
          }
   }
   return  false;
}

int main()
{
    int t;
    scanf("%d",&t);
    while( t-- )
    {
       init();
       int sig = 0;
       memset(vx,0,sizeof(vx));
       memset(vy,0,sizeof(vy));
       for( int i=1; i<=num1; ++i )
       {
            memset(flag,0,sizeof(flag));
            if( find(i) )
            {
                sig ++;
            }
       }
       printf("%d\n",p-sig);
    }
    return 0;
}

Smile3 2011-07-18 12:06 发表评论
]]>hnu第一场 Equivalent Sets 强连通分量 添加边使其成为一个强连通分量http://www.cppblog.com/Smile3/articles/151071.htmlSmile3Smile3Fri, 15 Jul 2011 09:22:00 GMThttp://www.cppblog.com/Smile3/articles/151071.htmlhttp://www.cppblog.com/Smile3/comments/151071.htmlhttp://www.cppblog.com/Smile3/articles/151071.html#Feedback0http://www.cppblog.com/Smile3/comments/commentRss/151071.htmlhttp://www.cppblog.com/Smile3/services/trackbacks/151071.html考察知识点：强连通分量的的求法
bug：注意强连通分量个数只有一个的情况
思维：这个转化倒是知道，但是知识点不足，还有一个更重要的是邻接表的写法和邻接矩阵的写法有很多的不同，特别是在第二次dfs的时候注意元素的表示不是i了而应该是
for( int i= 0; i<num; ++i)
       if( flag[ b[cur][i] ] == 0  )
          visttwo(b[cur][i],sig);
再贴一个tyvj的纯裸的强连通分量：tvvj1111
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define  max 210
int a[max][max],b[max][max];//a为图G的邻接矩阵，b为图G的反图的邻接矩阵
int  bleg[max];
int flag[max];
int numb[max];
int n;
void vistone(int cur,int &sig)
{
     flag[cur] = 1;
     for( int i=1; i<=n; ++i)
         if( flag[i] == 0 && a[cur][i] == 1 )
          vistone(i,sig );
     numb[++sig] = cur;
}
void visttwo( int cur ,int sig)
{
     flag[cur] = 1;
     bleg[cur] =sig;
     for( int i= 1; i<=n; ++i)
       if( flag[i] == 0 && b[cur][i] == 1 )
          visttwo(i,sig);
}     
int korareju()
{
    int sig =0;
    memset(flag,0,sizeof(flag));
    for( int i=1;i<=n; ++i)
       if( flag[i] == 0)
           vistone(i,sig);
    memset(flag,0,sizeof(flag));
    sig=0;
    for(int i=n; i>=1; --i )
       if( flag[ numb [i] ] == 0 )
       {
           visttwo(numb[i], (++sig) );
       }
    return sig;
}
int main()
{
    scanf("%d",&n);
    int i = 0;
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    while( (i++)<n )
    {
       while( true )
       {
           int x;
           scanf("%d",&x);
           if( x == 0)
              break;
           a[i][x] = 1;
           b[x][i] = 1;
       }
    }
    int w = korareju();
    printf("%d\n",w);
    return 0;
}
   


#include<cstdio>
#include<cstring>
#include<iostream>
#include<vector>
using namespace std;
#define  max 20010
vector< int > a[max],b[max];//a为图G的邻接矩阵，b为图G的反图的邻接矩阵
int  bleg[max];
int flag[max];
int numb[max];
int out[max];
int in[max];
int n;
void vistone(int cur,int &sig)
{
     flag[cur] = 1;
     int num = a[cur].size();
     for( int i=0; i<num; ++i)
         if( flag[ a[cur][i] ] == 0 )
              vistone(a[cur][i],sig );
     numb[++sig] = cur;
}
void visttwo( int cur ,int sig)
{
     flag[cur] = 1;
     bleg[cur] =sig;
     int num = b[cur].size();
     for( int i= 0; i<num; ++i)
       if( flag[ b[cur][i] ] == 0  )
          visttwo(b[cur][i],sig);
}     
int korareju()
{
    int sig =0;
    memset(flag,0,sizeof(flag));
    for( int i=1;i<=n; ++i)
       if( flag[i] == 0)
           {
             vistone(i,sig);
           }
    memset(flag,0,sizeof(flag));
    sig=0;
    for(int i=n; i>=1; --i )
       if( flag[ numb[i] ] == 0 )
       {
           visttwo( numb[i], (++sig) );
       }
    return sig;
}
int main()
{
    int m;
    while( scanf("%d%d",&n,&m) == 2 )
   {
    int k = 0;
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(numb,0,sizeof(numb));
    memset(out,0,sizeof(out));
    memset(in,0,sizeof(in));
    memset(bleg,0,sizeof(bleg));
    while( (k++)<m )
    {
           int x,y;
           scanf("%d%d",&x,&y);
           a[x].push_back(y);
           b[y].push_back(x);
    }
    if( n <= 1 )
    {
        printf("0\n");
        continue;
    }
    int w = korareju();
    if( w == 1 )
    {
        printf("0\n");
        continue;
    }
    for( int i=1; i<=n; ++i )
    {
      int num = a[i].size();
      for(int j=0; j<num; ++j)
         if( bleg[i] != bleg[ a[i][j] ] )
         {
             out[bleg[i]]++;
             in[ bleg[ a[i][j] ] ]++;
         }
    }
    int od=0;
    int id=0;
    for( int i= 1; i<=w; ++i)
    {
       if( out[i] == 0 )
          od++;
       if( in[i] == 0 )  
          id++;
    }
    int max0=od;
    if(max0<id)
      max0=id;   
    printf("%d\n",max0);
   }
    return 0;
}

Smile3 2011-07-15 17:22 发表评论
]]>hnu第一场 Earth Hour 《单源最短路》http://www.cppblog.com/Smile3/articles/151069.htmlSmile3Smile3Fri, 15 Jul 2011 08:51:00 GMThttp://www.cppblog.com/Smile3/articles/151069.htmlhttp://www.cppblog.com/Smile3/comments/151069.htmlhttp://www.cppblog.com/Smile3/articles/151069.html#Feedback0http://www.cppblog.com/Smile3/comments/commentRss/151069.htmlhttp://www.cppblog.com/Smile3/services/trackbacks/151069.html这是一个单源最短路的变形，枚举的是关节点；
小技巧：在开始用三次dijstra算法求出三个关键点到所有其他点的距离，保存下来，枚举关节点时就不要每一次都求解，以节省时间，否则会超时。
思维：这里卡的是转换过程，另外一个是知识点，对于知识点现在已经学过了，以后要加油，对于思维的转化，还是要靠经验的积累和考场的思考。
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define INF 10000000
int a[210][3];
int l[210][210];
int d1[210],d2[210],d3[210];
void initiate(int n)
{
   for( int i=0; i<n; ++i )
     for(int j=0; j<n; ++j )
     {
            int r = a[i][2]+a[j][2];
            int len = ( a[i][0] - a[j][0] )*( a[i][0] - a[j][0] ) + ( a[i][1] - a[j][1] )*( a[i][1] - a[j][1] );
            if( r*r >= len )
            {
                l[i][j] = 1;
                l[j][i] = 1;
            }
     }
     for(int i=0; i<n; ++i)
       l[i][i]=0;
}
void dijstra(int i,int n,int d[])//每一次选取已有集合的最小d值,更新其他d值
{
   d[i] = 0;
   int flag[210];
   memset(flag,0,sizeof(flag));
   int tt=n;
   while( tt-- )
   {
          int max =INF ;
          int num ;
          for(int j=0; j< n; j++ )
            if( d[j] < max && (flag[j] == 0) )
                 {
                     max = d[j];
                     num = j;                     
                 }
          flag[num] = 1;
          for( int t=0; t<n; ++t )
          {
            if( flag[t] == 0 && l[num][t] == 1 )
            {
                 if( d[t] > (d[num]+1) )
                 {
                     d[t] = d[num]+1;
                 }
            }
          }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while( t-- )
    {
      int n;
      scanf("%d",&n);
      memset(l,0,sizeof(l));
      //printf("n:%d\n",n);
      for( int i=0; i<n; ++i )
        scanf("%d%d%d",&a[i][0],&a[i][1],&a[i][2]);
      int w=n;
      initiate(w);
      for( int k=0; k<210; ++k )
      {
              d1[k] = INF;
              d2[k] = INF;
              d3[k] = INF;
      }
      dijstra(0,n,d1);
      dijstra(1,n,d2);
      dijstra(2,n,d3);
      int min =INF;
      for( int j=0; j<n; ++j)
           if( min > d1[j]+d2[j]+d3[j] )
             min = d1[j]+d2[j]+d3[j];
      if( min == INF )
         printf("-1\n");
      else
         printf("%d\n",(n-min-1));
    }
    return 0;
}

Smile3 2011-07-15 16:51 发表评论
]]>