BZOJ NOI十连测第二测 T1

出题人居然是个哲学家。。

26%的程序，太SB了。。。本来我的想法也是二分+贪心，但是贪心是个怪怪的SX贪心。。

 #include<algorithm>

 #include<cstdio>

 #include<cmath>

 #include<cstring>

 #include<iostream>

 #include<time.h>

 #define inf 0x7fffffff

 struct edge{

     int u,v,id;

 }e[];

 int tot,go[],first[],next[];

 int Tot,Go[],First[],Next[];

 int c[],b[],son[],ans,vis[],pd;

 int pass[],op[],C[],Cnt,block,sum,B[];

 int id1,id2;

 int n,k;

 int read(){

     char ch=getchar();int t=,f=;

     while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}

     while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}

     return t*f;

 }

 void Dfs(int x,int fa){

     if (vis[x]) return;

     c[x]=b[x];son[x]=;

     for (int i=first[x];i;i=next[i]){

         int pur=go[i];

         if (pur==fa) continue;

         if (pass[i]) continue;

         Dfs(pur,x);

         son[x]+=son[pur];

         c[x]+=c[pur];

     }

     if (c[x]&&(n-c[x]>)) ans=std::min(ans,std::max(son[x],n-son[x]));

 }

 void sbpianfen(){

     ans=inf;

     Dfs(,);

     printf("%d\n",ans);

 }

 void insert(int x,int y){

     tot++;

     go[tot]=y;

     next[tot]=first[x];

     first[x]=tot;

     pass[tot]=;

 }

 void add(int x,int y){

     insert(x,y);op[tot]=tot+;insert(y,x);op[tot]=tot-;

 }

 void Insert(int x,int y){

     Tot++;

     Go[Tot]=y;

     Next[Tot]=First[x];

     First[x]=Tot;

 }

 void Add(int x,int y){

     Insert(x,y);Insert(y,x);

 }

 int bfs(int x){

     int h=,t=;

     bool ck=;

     c[]=x;vis[x]=;if (b[x]) ck=;

     while (h<=t){

         int now=c[h++];

         for (int i=First[now];i;i=Next[i]){

             int pur=Go[i];

             if (vis[pur]) continue;

             c[++t]=pur;

             vis[pur]=;

             if (b[pur]) ck=;

         }

     }

     if (!ck) pd=;

     return t;

 }

 void dFs(int x){

     if (x==n){

         pd=;

         int bl=;

         for (int i=;i<=n;i++)

          First[i]=;

         Tot=;

         for (int i=;i<n;i++)

          if (!e[i].id) Add(e[i].u,e[i].v);

         for (int i=;i<=n;i++)

           vis[i]=;

         for (int i=;i<=n;i++)

          if (!vis[i])

           bl=std::max(bl,bfs(i));

         if (pd==) ans=std::min(ans,bl);

         return;

     }

     e[x].id=;

     dFs(x+);

     e[x].id=;

     dFs(x+);

 }

 void sxpianfen(){

     for (int i=;i<n;i++) e[i].id=;

     dFs();

 }

 void clear(int x){

     int h=,t=;

     C[]=x;

     bool ppp=;

     if (b[x]) Cnt++,ppp=;

     vis[x]=;b[x]=;

     while (h<=t){

         int now=C[h++];

         for (int i=first[now];i;i=next[i]){

             int pur=go[i];

             if (pass[i]) continue;

             if (vis[pur]) continue;

             C[++t]=pur;

             vis[pur]=;

             if (b[pur]) Cnt++,ppp=;

             b[pur]=;

         }

     }

     if (ppp) {

         pd=;

         return;

     }

     if (t>block) pd=;

     sum-=t;

 }

 void work(int x,int mid,int fa,int Edge){

     if (c[x]==&&b[fa]==&&son[x]<=block) {

         pass[Edge]=;

         pass[op[Edge]]=;

         clear(x);

         id1=x;id2=fa;

         if (pd) return;

         return;

     }

     else

     if (k-Cnt-c[x]==&&b[x]==&&sum-son[x]<=block){

         pass[Edge]=;

         pass[op[Edge]]=;

         clear(fa);

         id1=x;id2=fa;

         if (pd) return;

         return;

     }

     else

     if (c[x]==&&son[x]<=block&&son[fa]>=block){

         pass[Edge]=;

         pass[op[Edge]]=;

         clear(x);

         id1=x;id2=fa;

         if (pd) return;

         return;

     }else

     if (k-Cnt-c[x]==&&sum-son[x]==block){

         pass[Edge]=;

         pass[op[Edge]]=;

         clear(fa);

         id1=x;id2=fa;

         if (pd) return;

         return;

     }

     for (int i=first[x];i;i=next[i]){

         int pur=go[i];

         if (pur==fa) continue;

         if (vis[pur]) continue;

         if (pass[i]) continue;

         work(pur,mid,x,i);

     }

 }

 bool check(int mid){

     for (int i=;i<n;i++) std::swap(e[rand()%(n-)+],e[rand()%(n-)+]);

     for (int i=;i<=n;i++) first[i]=;

     tot=;

     for (int i=;i<n;i++)

      add(e[i].u,e[i].v);

     block=mid;pd=;

     for (int i=;i<=tot;i++)

      pass[i]=;

     for (int i=;i<=n;i++)

      b[i]=B[i],vis[i]=;

     id1=;id2=;vis[]=;

     sum=n;Cnt=;

     for (int i=;i<=k;i++){

       if (vis[id1]&&vis[id2]) break;

       if (k-Cnt==) break;

       if (!vis[id1]){

        Dfs(id1,);

        work(id1,mid,,);

       }else{

        Dfs(id2,);

        work(id2,mid,,);

       }

       if (pd) return ;

     }

     for (int i=;i<=n;i++)

      if (!vis[i]) clear(i);

     if (pd) return ;

     return ;

 }

 void solve(){

     int l=n/k,r=n,Ans=n;

     while (l<=r){

         int mid=(l+r)/;

         if (check(mid)) r=mid-,Ans=mid;

         else l=mid+;

     }

     printf("%d\n",Ans);

 }

 void solve2(){

     int l=n/k,r=n,Ans=n;

     while (l<=r){

         int mid=(l+r)/;

         if (check(mid)) r=mid-,Ans=mid;

         else l=mid+;

     }

     int Tn=Ans;

     l=n/k,r=n,Ans=n;

     while (l<=r){

         int mid=(l+r)/;

         if (check(mid)) r=mid-,Ans=mid;

         else l=mid+;

     }

     Tn=std::min(Tn,Ans);

     printf("%d\n",Tn);

 }

 int main(){

     //♂

     freopen("deep.in","r",stdin);

     freopen("deep.out","w",stdout);

     srand(time(NULL));

     ans=inf;

     n=read();k=read();

     for (int i=;i<n;i++){

         int x=read(),y=read();

         add(x,y);

         e[i].u=x;e[i].v=y;

     }

     for (int i=;i<=k;i++){

         int x=read();

         b[x]=;B[x]=;

     }

     if (n<=) {sxpianfen();printf("%d\n",ans);return ;}

     if (k==) {printf("%d\n",n);return ;}

     if (k==) {sbpianfen();return ;}

     if (n<=) {solve();return ;}

     if (k==) {solve2();return ;}

     if (k==n) {printf("%d\n",);return ;}

     if (k==n-) {printf("%d\n",);return ;}

     solve();

 }

AC程序不知道短到哪里去了

 #include<algorithm>

 #include<cstdio>

 #include<cmath>

 #include<cstring>

 #include<iostream>

 #define N 200005

 int tot,go[N*],next[N*],first[N];

 int n,k,b[N],fa[N],mx[N],sum[N],vis[N],c[N];

 int read(){

     char ch=getchar();int t=,f=;

     while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}

     while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}

     return t*f;

 }

 void insert(int x,int y){

     tot++;

     go[tot]=y;

     next[tot]=first[x];

     first[x]=tot;

 }

 void add(int x,int y){

     insert(x,y);insert(y,x);

 }

 void bfs(){

     int h=,t=;for (int i=;i<=n;i++) vis[i]=;

     c[]=;

     vis[]=;

     while (h<=t){

         int now=c[h++];

         for (int i=first[now];i;i=next[i]){

             int pur=go[i];

             if (vis[pur]) continue;

             vis[pur]=;

             c[++t]=pur;

             fa[pur]=now;

         }

     }

 }

 bool check(int mid){

     int f;

     memset(sum,,sizeof sum);

     memset(mx,,sizeof mx);

     for (int i=n;i>=;i--){

         int u=c[i];

         if (b[u]) mx[u]=mid-;

         else sum[u]++;

         f=(mx[u]>=sum[u])?mx[u]-sum[u]:-sum[u];

         if (u==&&f<) return ;

         if (f<) sum[fa[u]]-=f;

         else mx[fa[u]]=std::max(mx[fa[u]],f);

     }

     return ;

 }

 int main(){

     n=read();k=read();

     for (int i=;i<n;i++){

         int x=read(),y=read();

         add(x,y);

     }

     for (int i=;i<=k;i++){

         int x=read();

         b[x]=;

     }

     bfs();

     int l=n/k,r=n,ans;

     while (l<=r){

         int mid=(l+r)>>;

         if (check(mid)) r=mid-,ans=mid;

         else l=mid+;

     }

     printf("%d\n",ans);

 }

BZOJ NOI十连测第二测 T1的更多相关文章

BZOJ NOI十连测第二测 T2
思路:20%可以搜索.. #include<algorithm> #include<cstdio> #include<cmath> #include<cstr ...
BZOJ NOI十连测第一测 T1
思路:首先考虑t=1的情况,t等于1,那么所有位置的颜色相同,我们不用考虑概率的问题,那么,k+d*x在模d下都相等,我们考虑预处理一个数组s[i][j],代表d为i,起始位置为j的等差数列的和,这个 ...
BZOJ NOI十连测第一测 T2
思路:看到这题,就感觉是一道很熟悉的题目: http://www.cnblogs.com/qzqzgfy/p/5535821.html 只不过这题的K最多可以到N,而且边权不再只是1,考试的时候yy了 ...
痞子衡嵌入式：测一测i.MXRT1170 Raw NAND启动时间(从POR到进App的Reset_Handler)
大家好,我是痞子衡,是正经搞技术的痞子.今天痞子衡给大家介绍的是恩智浦i.MX RT1170 Raw NAND启动时间. 关于i.MXRT1170这颗划时代的MCU,痞子衡去年10月在其刚发布的时候, ...
「NOI十联测」深邃
「NOI十联测」深邃要使得最大的连通块最小,显然先二分答案. 先固定1结点为根. 对于一个果实,显然是先处理子树中未分配的点,再向外延伸. 每个结点记录一个\(si[]\),表示子树中未分配的点数, ...
「NOI十联测」奥义商店
「NOI十联测」奥义商店若lzz想花费最少的钱,那么显然要选择数目较少的颜色. 先考虑暴力的写法. 每次向两边统计,每个物品要求被买的概率可以由上一个物品推出. now=1;//now 被买概率 M ...
「NOI十联测」黑暗
「NOI十联测」黑暗 \(n\) 个点的无向图,每条边都可能存在,一个图的权值是连通块个数的 \(m\) 次方,求所有可能的图的权值和.(n≤30000,m≤15) 令\(ans[n][m]\)为n个 ...
NOI十连测第六测 T1
思路: 用treap动态维护,记一个sum1,sum2,注意!,写treap如果有删除操作,千万不能把权值相同的分开来..,这在删除的时候会进入死循环,这是一个惨痛的教训... #include< ...
NOI十连测第五测 T1
#include<cstdio> #include<cstring> #include<cmath> #include<iostream> #inclu ...

随机推荐

js django 数据互动
后台传输到客户端的数据怎么跟js进行互动 eval("("+{{data|safe}}+")") 获取到的数据可以直接使用 eval 有时候会存在编码问题所以推 ...
Codeforces 543A Writing Code
http://codeforces.com/problemset/problem/543/A 题目大意:n个人,一共要写m行程序,每个程序员每行出现的bug数为ai,要求整个程序出现的bug数不超过b ...
《Programming WPF》翻译第3章 4.我们进行到哪里了？
原文:<Programming WPF>翻译第3章 4.我们进行到哪里了? 控件是由应用程序创建的块.它们描述了用户用来交互的界面特征.控件提供了行为,依赖样式和模板来表示一个外观.输入 ...
【HDU1875】畅通工程再续（MST基础题）
更改成实形数即可.第一次敲完直接交,CE了一次.晕. #include <iostream> #include <cstring> #include <cstdio> ...
Java配置文件Properties的读取、写入与更新操作
/** * 实现对Java配置文件Properties的读取.写入与更新操作 */ package test; import java.io.BufferedInputStream; import j ...
debian下编译libev库
系统为Linux debian 2.6.32-5-686.这是裸系统,连xwindows都没有.帐户为root,不是的注意一下权限.这里想说明安装过程及出现的问题,故打印的信息较多,以供出现错误的读者 ...
cStringIO模块例子
# Vorbis comment support for Mutagen # Copyright 2005-2006 Joe Wreschnig # # This program is free so ...
Ruby新手教程和技巧
Ruby真的比Java更好? Ruby On Rails 创始人:对Java 说再见这两周以来环绕Java发生的两件大事:EclipseCon 和TheServerSide Java Sympo ...
向html某个元素中添加信息
<!DOCTYPE html> <html> <head> <meta charset="utf-8"> <meta http ...
网络基础知识HTTP(1) --转载
为什么要写网络? 作为网站开发人员,你所开发的软件产品最终是要在网络上运行的.这就像一个生产商,要生产供给东北地区的产品,而生产商对东北的天气.地理.人文毫无了解.生产商的产品肯定是不可用的,或者低端 ...

BZOJ NOI十连测 第二测 T1

BZOJ NOI十连测 第二测 T1的更多相关文章

随机推荐

热门专题

BZOJ NOI十连测第二测 T1

BZOJ NOI十连测第二测 T1的更多相关文章