【noip2015】【Day2】【二分】【倍增】运输计划

唉一个懒得写博客的我又来了…
懒得贴题目了，自己找找吧…
题目请戳这里orz（来源：vijos）

这道题呢，可以分两部分来做：
1.基础骗分…
2.进一步骗分….

千万要好好利用数据神奇的特点√
于是乎写链状嘛…我没写出来….但后来在修改程序的时候…写了个出来对拍orz…蠢死…..

---

思想：
1.二分答案，这个应该是很好理解的吧√
2.将不符合的路线整理出来，起点终点分别++，它们的lca-=2；
自下而上回溯“求和”求出整棵树的值，每个结点的值就表示这个结点往上的那一条边经过的次数√。貌似是叫拆分的思想吧….
3.求lca要用倍增来缩短时间√用tarjan不会超时，那个树规的那种倍增就会超时qwq
4.找出每个 值=不符合的路线数 的结点，并选出边权最大的，将所有路线的ans-max，若<=mid就成立咯
5.预处理要处理出每个方案的lca，每个方案的总长度。

---

链状代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#define LL long long
using namespace std;
const int N=300005;
int pa[N][2],zhi[N],len[N],sum[N],a,b,c,ma,n,m;
LL ans;
bool chek(int x)
{
    memset(sum,0,sizeof(sum));
    int tt=0,xu=0;
    bool yes=false;
    for(int i=1;i<=m;i++)
    if(len[i]>x){
        tt++;sum[pa[i][0]]++;sum[pa[i][1]]--;
    }
    for(int i=2;i<=n;i++)
    {
        sum[i]+=sum[i-1];
        if(sum[i]==tt&&zhi[i]>=ma-x){yes=true;break;}
    }
    return yes;
}
int main()
{
    freopen("transport.in","r",stdin);
    freopen("transport2.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n-1;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        if(a>b)swap(a,b);
        zhi[a]=c;
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&a,&b);
        if(a>b)swap(a,b);
        pa[i][0]=a;pa[i][0]=b;
        while(a!=b){len[i]+=zhi[a];a++;}
        ma=max(len[i],ma);
    }
    int le=0,ri=ma;
    while(le<=ri)
    {
        int mid=(le+ri)>>1;
        if(chek(mid))
        {ans=mid;ri=mid-1;}
        else le=mid+1;
    }
    cout<<ans<<endl;
}

---

下面贴上正确代码，因为不会tarjan版的倍增，只能过80分orz

#include<cstring>
#include<string>
#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
const int N=300005;
int n,m,fa[N][20],zh[N],dp[N],c;//2^19>1e5
int head[N*2],next[N*2],tov[N*2],zhi[N*2],tot,ma;
int poi[N][2],lca[N],len[N],tt,flag,sum[N],ans;
bool sym[N],yes;
void build(int a,int b,int c)
{
    tot++;next[tot]=head[a];tov[tot]=b;head[a]=tot;zhi[tot]=c;
    tot++;next[tot]=head[b];tov[tot]=a;head[b]=tot;zhi[tot]=c;
}
void dfs(int x)
{
    sym[x]=true;
    int u=head[x];
    while(u)
    {
        int v=tov[u];
        if(!sym[v])
        {
            zh[v]=zhi[u];
            dp[v]=dp[x]+1;
            fa[v][0]=x;
            dfs(v);
        }
        u=next[u];
    }
}
void dfs2(int x)
{
    if(yes)return;
    sym[x]=true;
    int u=head[x];
    while(u)
    {
        int v=tov[u];
        if(!sym[v])
        {
            dfs2(v);
            sum[x]+=sum[v];
        }
        u=next[u];
    }
    if(sum[x]==tt&&zh[x]>=flag)yes=true;
}
int LCA(int x)
{
    int a=poi[x][0],b=poi[x][1];
    int w=dp[b]-dp[a],i=0;
    while(w&&i<20)
    {
        if(w&1<<i)  {b=fa[b][i];w=w^1<<i;}
        i++;
    }
    if(a==b)return a;
    for(int i=19;i>=0;i--)
      if(fa[a][i]!=fa[b][i])
      {a=fa[a][i];b=fa[b][i];}
    a=fa[a][0];b=fa[b][0];
    return a;
}
void pre()
{
    scanf("%d%d",&n,&m);
    int a,b,c;
    for(int i=1;i<=n-1;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        build(a,b,c);
    }
    dfs(1);
    for(int i=1;i<20;i++)  //bz
      for(int j=1;j<=n;j++)
      fa[j][i]=fa[fa[j][i-1]][i-1];
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&a,&b);
        if(dp[a]>dp[b])swap(a,b);
        poi[i][0]=a;poi[i][1]=b;
        lca[i]=LCA(i);
        while(a!=lca[i]){len[i]+=zh[a];a=fa[a][0];}
        while(b!=lca[i]){len[i]+=zh[b];b=fa[b][0];}
        ma=max(ma,len[i]);
    }
}
bool chek(int x)
{
    tt=0;yes=false;
    memset(sum,0,sizeof(sum));
    memset(sym,false,sizeof(sym));
    for(int i=1;i<=m;i++)
        if(len[i]>x)
        {tt++;sum[poi[i][0]]++;sum[poi[i][1]]++;sum[lca[i]]-=2;}
    flag=ma-x;dfs2(1);
    return yes;
}
void workk()
{
    int le=0,ri=ma;
    while(le<=ri)
    {
        int mid=(le+ri)>>1;
        if(chek(mid))
        {ans=mid;ri=mid-1;}
        else le=mid+1;
    }
    if(!ans)ans=ma;
    printf("%d\n",ans);
}
int main()
{
    freopen("transport.in","r",stdin);
    freopen("transport.out","w",stdout);
    pre();
    workk();
}