• POJ 1741/1987 树的点分治


    树的点分治,主要思想是每次找子树的重心,计算经过根节点的情况数,再减去点对属于同一子树的情况。

      1 #include <iostream>
      2 #include <vector>
      3 #include <algorithm>
      4 #include <string>
      5 #include <string.h>
      6 #include <stdio.h>
      7 #include <stdlib.h>
      8 #include <queue>
      9 #include <stack>
     10 #include <map>
     11 #include <set>
     12 #include <cmath>
     13 #include <ctime>
     14 #include <cassert>
     15 #include <sstream>
     16 using namespace std;
     17 
     18 const int N=10010;
     19 struct Edge {
     20     int to,next;
     21     int w;
     22     Edge(){}
     23     Edge (int _t,int _n,int _w=0) {
     24         to=_t;
     25         next=_n;
     26         w=_w;
     27     }
     28 }edge[N<<1];
     29 int idx,head[N];
     30 void addEdge (int u,int v,int w) {
     31     ++idx;
     32     edge[idx]=Edge(v,head[u],w);
     33     head[u]=idx;
     34 }
     35 bool vis[N];
     36 
     37 int K;// 输入中的k
     38 
     39 int dis[N],disCnt;
     40 void getDis(int u,int f,int d) {
     41     dis[disCnt++]=d;
     42     for (int k=head[u];k;k=edge[k].next) {
     43         int v=edge[k].to;
     44         if (vis[v]||v==f) continue;
     45         getDis(v,u,d+edge[k].w);
     46     }
     47 }
     48 
     49 int bal,cmp;
     50 int getBal(int u,int f) {
     51     int son=1;
     52     int dp=0;
     53     for (int k=head[u];k;k=edge[k].next) {
     54         int v=edge[k].to;
     55         if (vis[v]||v==f) continue;
     56         int subSon=getBal(v,u);
     57         son+=subSon;
     58         dp=max(dp,subSon);
     59     }
     60     dp=max(dp,disCnt-son);
     61     if (dp<cmp){
     62         cmp=dp;
     63         bal=u;
     64     }
     65     return son;
     66 }
     67 
     68 int calc(int u,int initDis) {
     69     disCnt=0;
     70     getDis(u,-1,initDis);
     71     sort(dis,dis+disCnt);
     72     int ret=0;
     73     int l=0,r=disCnt-1;
     74     while (l<r) {
     75         if (dis[l]+dis[r]>K) r--;
     76         else ret+=r-l++;
     77     }
     78     return ret;
     79 }
     80 
     81 int ans=0;
     82 void solve(int u) {
     83     cmp=~0U>>1;
     84     getBal(u,-1);
     85     ans+=calc(bal,0);
     86     vis[bal]=true;
     87     for (int k=head[bal];k;k=edge[k].next) {
     88         int v=edge[k].to;
     89         if (vis[v]) continue;
     90         ans-=calc(v,edge[k].w);// 减去子树内的重复计算的情况
     91         solve(v);
     92     }
     93 }
     94 
     95 void init(int n) {
     96     idx=1;memset(head,0,sizeof head);
     97     memset(vis,false,sizeof vis);
     98     disCnt=n;// 由于在计算重心时,可以直接用之前算dis时的disCnt,所以这里初始化为n
     99     ans=0;
    100 }
    101 int main () {
    102     int n;
    103     while (scanf("%d %d",&n,&K)!=EOF) {
    104         if (n==0&&K==0) break;
    105         init(n);
    106         for (int i=1;i<n;i++) {
    107             int u,v,w;
    108             scanf("%d %d %d",&u,&v,&w);
    109             addEdge(u,v,w);
    110             addEdge(v,u,w);
    111         }
    112         solve(1);
    113         printf("%d
    ",ans);
    114     }
    115     return 0;
    116 }
  • 相关阅读:
    源码解析之–网络层YTKNetwork
    ARC和非ARC文件混编
    Xcode真机调试中"There was an internal API error"错误解决方法
    设置textView或者label的行间距方法
    iOS中__block 关键字的底层实现原理
    使用AVCaptureSession捕捉静态图片
    使用AVCaptureSession显示相机预览
    短小强悍的JavaScript异步调用库
    开源中国愚人节网页变模糊的js blur代码
    undefined与null的区别
  • 原文地址:https://www.cnblogs.com/micrari/p/4803021.html
Copyright © 2020-2023  润新知