bzoj 3697

 

终于自己独立A了一道点分治。

点分治的题目和树规的题目很像,都是在树上乱搞,但点分治一般和路径更相关,如果用树规做会感觉冗余状态太多,内存和时间都无法承受,如果不用动态规划,直接在原树上运用这道题的方法,又无法保证时间复杂度(点分治让我们的算法对原树的形态依赖更小,可以保证时间复杂度)。

 

  1 /**************************************************************
  2     Problem: 3697
  3     User: idy002
  4     Language: C++
  5     Result: Accepted
  6     Time:2480 ms
  7     Memory:14012 kb
  8 ****************************************************************/
  9  
 10 #include <cstdio>
 11 #define N 100010
 12 #define fprintf(...)
 13  
 14 int n;
 15 int head[N], dest[N+N], wght[N+N], next[N+N], ntot;
 16 int fat[N], vis[N], siz[N], bac[N], dis[N], cur_root, cur_siz;
 17 int vv[N+N], svv[N+N][2], cvv[N+N][2];
 18 int stkc[N+N], tpc, stks[N+N], tps;
 19 long long ans;
 20  
 21 void insert( int u, int v, int w ) {
 22     ntot++;
 23     next[ntot] = head[u];
 24     wght[ntot] = w ? 1 : -1;
 25     dest[ntot] = v;
 26     head[u] = ntot;
 27 }
 28 void dfs_root( int u ) {
 29     siz[u] = 1, bac[u] = 0;
 30     for( int t=head[u]; t; t=next[t] ) {
 31         int v=dest[t];
 32         if( vis[v] || v==fat[u] ) continue;
 33         fat[v] = u;
 34         dfs_root(v);
 35         siz[u] += siz[v];
 36         if( siz[v]>bac[u] ) bac[u]=siz[v];
 37     }
 38     if( cur_siz-siz[u]>bac[u] ) bac[u]=cur_siz-siz[u];
 39     if( bac[u]<bac[cur_root] ) cur_root=u;
 40 }
 41 void dfs_dis( int u ) {
 42     if( dis[u]==0 && vv[N+0]>1 )
 43         ans++;
 44     bool ok = vv[N+dis[u]];
 45     if( ok ) {
 46         ans += svv[N-dis[u]][0]+svv[N-dis[u]][1];
 47         cvv[N+dis[u]][1]++;
 48     } else {
 49         ans += svv[N-dis[u]][1];
 50         cvv[N+dis[u]][0]++;
 51     }
 52     stkc[++tpc] = dis[u];
 53     fprintf( stderr, "Arrived %d cur_ans = %lld\n", u, ans );
 54     vv[N+dis[u]]++;
 55     for( int t=head[u]; t; t=next[t] ) {
 56         int v=dest[t], w=wght[t];
 57         if( vis[v] || v==fat[u] ) continue;
 58         dis[v] = dis[u]+w;
 59         fat[v] = u;
 60         dfs_dis( v );
 61     }
 62     vv[N+dis[u]]--;
 63 }
 64 void build_vdcp( int rt ) {
 65     /*   found the core   */
 66     cur_siz = siz[rt];
 67     cur_root = 0;
 68     fat[rt] = rt;
 69     dfs_root( rt );
 70     /*   collect the info  */
 71     rt = cur_root;
 72     vis[rt] = true;
 73     fprintf( stderr, "cur_root = %d\n", rt );
 74     for( int t=head[rt]; t; t=next[t] ) {
 75         int u=dest[t], w=wght[t];
 76         if( vis[u] ) continue;
 77         tpc = 0;
 78         vv[N+0]++;
 79         fat[u] = rt;
 80         dis[u] = w;
 81         dfs_dis( u );
 82         vv[N+0]--;
 83         while( tpc ) {
 84             int v = stkc[tpc--];
 85             if( !cvv[N+v][0] && !cvv[N+v][1] ) continue;
 86             svv[N+v][0] += cvv[N+v][0];
 87             svv[N+v][1] += cvv[N+v][1];
 88             stks[++tps] = v;
 89             cvv[N+v][0] = cvv[N+v][1] = 0;
 90         }
 91     }
 92     while( tps ) {
 93         int v = stks[tps--];
 94         svv[N+v][0] = svv[N+v][1] = 0;
 95     }
 96     for( int t=head[rt]; t; t=next[t] ) {
 97         int u=dest[t];
 98         if( vis[u] ) continue;
 99         build_vdcp( u );
100     }
101 }
102 int main() {
103     scanf( "%d", &n );
104     for( int i=1,u,v,w; i<n; i++ ) {
105         scanf( "%d%d%d", &u, &v, &w );
106         insert( u, v, w );
107         insert( v, u, w );
108     }
109     bac[0] = n;
110     siz[1] = n;
111     build_vdcp(1);
112     printf( "%lld\n", ans );
113 }
View Code

 

posted @ 2015-03-28 13:13  idy002  阅读(...)  评论(...编辑  收藏