SDTSC 2010 sotomon

SDTSC 2010 sotomon

题意：
给定N个点 有些节点可以通往同行 有些可以通往同列中的点 有些可以通往八连通的点 可以走过多次
问最多一次可以走过多少点

做法：
可以走过多次也就是说跟apio2009 atm一样 一个强连通分量内可以无限走
然后就是走过场 缩点 拓扑排序 dp

问题在于。。。离散化后每次要寻找(x,y)这个点是否存在  然后我的二分查找就因为小于大于号的问题挂了。。

（跟cqtsc2010 内部白点一样。。。挂在二分上了！！！）

  1  #include  < cstdio >
  2  #include  < algorithm >
  3  using   namespace  std;
  4  #define  n 300005
  5  #define  m 3000005
  6  struct  Tpnt
  7  {
  8       int  x,y,kind,o;
  9  }    T[ 100005 ];
 10  int  vtx[m],ne[m],L[n],tot,Tot,All,Sub,E;
 11  int  N,R,C,dfn[n],F[n],cnt[n],low[n],Stk[n],Deg[n],sub[n],p[n];
 12  int  x[ 100005 ],y[ 100005 ],cnt_x[ 1000005 ],place_x[ 1000005 ],place_y[ 1000005 ];
 13  int  X[m],Y[m];
 14  bool  vis[n];
 15  inline  bool  cmp_x( const  Tpnt  & a, const  Tpnt  & b)
 16  {
 17       return  a.x < b.x || a.x == b.x && a.y < b.y; 
 18  }
 19  inline  bool  cmp_y( const  Tpnt  & a, const  Tpnt  & b)
 20  {
 21       return  a.y < b.y || a.y == b.y && a.x < b.x;
 22  }
 23  inline  void  Ins( int  u, int  v)
 24  {
 25      vtx[ ++ tot] = v;ne[tot] = L[u];L[u] = tot;
 26  }
 27  inline  int  findx_x( int  X)
 28  {
 29       int  l = 0 ,r = N,mid;
 30       for  (;l + 1 < r;)
 31       if  (mid = (l + r) >> 1 ,T[mid].x < X)    l = mid;
 32       else     r = mid;
 33       if  (T[r].x == X)     return  r;
 34       return   - 1 ;
 35  }
 36  inline  int  findx_y( int  st, int  en, int  Y)
 37  {
 38       int  l = st - 1 ,r = en,mid;
 39       for  (;l + 1 < r;)
 40       if  (mid = (l + r) >> 1 ,T[mid].y < Y)    l = mid;
 41       else     r = mid;
 42       if  (T[r].y == Y)     return  r;
 43       return   - 1 ;
 44  }
 45  inline  void  Tarjan( int  u)
 46  {
 47      dfn[u] = low[u] =++ All;
 48      vis[Stk[ ++ Stk[ 0 ]] = u] = 1 ;
 49       for  (p[u] = L[u];p[u];p[u] = ne[p[u]])
 50       if  ( ! dfn[vtx[p[u]]])    Tarjan(vtx[p[u]]),low[u] = min(low[u],low[vtx[p[u]]]);
 51       else
 52       if  (vis[vtx[p[u]]])    low[u] = min(low[u],dfn[vtx[p[u]]]);
 53       if  (dfn[u] == low[u])
 54           for  ( ++ Sub;;)
 55          {
 56              vis[Stk[Stk[ 0 ]]] = 0 ;
 57              sub[Stk[Stk[ 0 ]]] = Sub;
 58               if  (Stk[Stk[ 0 ] -- ] == u)     break ;
 59          }
 60  }
 61  inline  void  Topo()
 62  {
 63      Stk[ 0 ] = 0 ;
 64       for  ( int  i = 1 ;i <= Sub; ++ i)
 65       if  ( ! Deg[i])    F[i] = cnt[i],Stk[ ++ Stk[ 0 ]] = i;
 66       for  (;Stk[ 0 ];)
 67      {
 68           int  u = Stk[Stk[ 0 ] -- ];
 69           for  ( int  p = L[u],v = vtx[p];p;v = vtx[p = ne[p]])
 70          {
 71              F[v] = max(F[v],F[u] + cnt[v]);
 72               -- Deg[v];
 73               if  ( ! Deg[v])    Stk[ ++ Stk[ 0 ]] = v;
 74          }
 75      }
 76  }
 77  int  main()
 78  {
 79      scanf( " %d%d%d " , & N, & R, & C);
 80       for  ( int  i = 1 ;i <= N; ++ i)
 81          scanf( " %d%d%d " , & T[i].x, & T[i].y, & T[i].kind),T[i].o = i;
 82      sort(T + 1 ,T + N + 1 ,cmp_x);
 83      E = 0 ;
 84       for  ( int  i = 1 ;i <= N; ++ i)
 85      {
 86           if  (T[i].x != T[i - 1 ].x)
 87          {
 88              x[ ++ x[ 0 ]] = T[i].x;
 89              place_x[T[i].x] = x[ 0 ];
 90          }
 91           ++ cnt_x[T[i].x];
 92          Ins(x[ 0 ] + N,T[i].o);
 93          X[ ++ E] = x[ 0 ] + N,Y[E] = T[i].o;
 94      }
 95      sort(T + 1 ,T + N + 1 ,cmp_y);
 96       for  ( int  i = 1 ;i <= N; ++ i)
 97      {
 98           if  (T[i].y != T[i - 1 ].y)
 99          {
100              y[ ++ y[ 0 ]] = T[i].y;
101              place_y[T[i].y] = y[ 0 ];
102          }
103          Ins(y[ 0 ] + N + x[ 0 ],T[i].o);
104          X[ ++ E] = y[ 0 ] + N + x[ 0 ],Y[E] = T[i].o;
105      }
106      Tot = N + x[ 0 ] + y[ 0 ];
107      sort(T + 1 ,T + N + 1 ,cmp_x);
108       for  ( int  i = 1 ;i <= N; ++ i)
109       if  (T[i].kind == 1 )
110      {
111          Ins(T[i].o,place_x[T[i].x] + N);
112          X[ ++ E] = T[i].o,Y[E] = place_x[T[i].x] + N;
113      }
114       else
115       if  (T[i].kind == 2 )
116      {
117          Ins(T[i].o,place_y[T[i].y] + N + x[ 0 ]);
118          X[ ++ E] = T[i].o,Y[E] = place_y[T[i].y] + N + x[ 0 ];
119      }
120       else
121           for  ( int  dx =- 1 ;dx <= 1 ; ++ dx)
122           if  (T[i].x + dx > 0 && T[i].x + dx <= R)
123          {
124               int  xst = findx_x(T[i].x + dx),xen;
125               if  (xst < 0 )     continue ;
126              xen = xst + cnt_x[T[i].x + dx] - 1 ;
127               for  ( int  dy =- 1 ;dy <= 1 ; ++ dy)
128               if  (T[i].y + dy > 0 && T[i].y + dy <= C)
129              {
130                   if  ( ! dx &&! dy)     continue ;
131                   int  pos = findx_y(xst,xen,T[i].y + dy);
132                   if  (pos < 0 || i == pos)     continue ;
133                  Ins(T[i].o,T[pos].o);
134                  X[ ++ E] = T[i].o,Y[E] = T[pos].o;
135              }
136          }
137       for  ( int  i = 1 ;i <= Tot; ++ i)
138       if  ( ! dfn[i])    Tarjan(i);
139       for  ( int  i = 1 ;i <= N; ++ i)
140           ++ cnt[sub[i]];
141      tot = 0 ;
142      memset(L, 0 , sizeof (L));
143       for  ( int  i = 1 ;i <= E; ++ i)
144       if  (sub[X[i]] != sub[Y[i]])
145      {
146          Ins(sub[X[i]],sub[Y[i]]);
147           ++ Deg[sub[Y[i]]];
148      }
149      Topo();
150       int  ret = 0 ;
151       for  ( int  i = 1 ;i <= Sub; ++ i)
152          ret = max(ret,F[i]);
153      printf( " %d\n " ,ret);
154       return   0 ;
155  }
156