• bzoj1023: [SHOI2008]cactus仙人掌图


    传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1023

    思路:类似树形DP记录一个f[i]表示最远(因为有环,所以这个定义是有一些限制条件的)

    先用点双缩点,每个环的信息可以挂到最高点上

    树上的差不多

    对于环上的,从一边扫过去,因为dis(i,j)有单调性,用单调队列搞一搞即可

    更详细的题解:http://ydcydcy1.blog.163.com/blog/static/21608904020131493113160/


    <span style="font-size:14px;">#include<cstdio>
    #include<vector>
    #include<cstring>
    #include<iostream>
    #include<algorithm>
    #define abs(a) (a<0?(-(a)):a)
    const int maxn=100010,maxm=500010;
    typedef unsigned int uint;
    using namespace std;
    int n,m,low[maxn],dfn[maxn],tim,sta[maxn],top,bnm[maxn],bcnt,ans,q[maxn],head,tail,a[maxn],f[maxn];bool vis[maxm];
    vector<int> bel[maxn],scc[maxn];
    
    struct Tgraph{
    	int pre[maxm],now[maxn],son[maxm],tot;
    	void add(int a,int b){pre[++tot]=now[a],now[a]=tot,son[tot]=b;}
    	void ins(int a,int b){add(a,b),add(b,a);}
    	void dp1(int x,int fa){
    		int id=x-n,cnt=scc[id].size(),top=0,siz=0;
    		for (int i=0;i<cnt;i++) a[++siz]=scc[id][i];
    		for (int i=0;i<cnt;i++) a[++siz]=scc[id][i];
    		head=1,tail=0;
    		for (int i=1;i<=siz;i++){
    			while (head<=tail&&i-q[head]>cnt/2) head++;
    			if (head<=tail) ans=max(ans,f[a[i]]+f[a[q[head]]]+i-q[head]);
    			while (head<=tail&&f[a[q[tail]]]-q[tail]<f[a[i]]-i) tail--;
    			q[++tail]=i;if (a[i]==fa&&!top) top=i-1;
    		}
    		for (int i=0;i<cnt;i++){
    			int d=abs(i-top);if (d>cnt/2) d=cnt-d;
    			f[x]=max(f[x],f[scc[id][i]]+d);
    		}
    		//printf("ans%d
    ",ans);
    	}
    	void dp2(int x,int fa){
    		int max1=0,max2=0;
    		for (int y=now[x];y;y=pre[y]) if (son[y]!=fa){
    			max2=max(max2,f[son[y]]+(son[y]<=n));
    			if (max2>max1) swap(max1,max2);
    		}
    		f[x]=max1;ans=max(ans,max1+max2);
    	}
    	void tree_dp(int x,int fa){
    		for (int y=now[x];y;y=pre[y]) if (fa!=son[y]) tree_dp(son[y],x);
    		if (x>n) dp1(x,fa);else dp2(x,fa);
    	}
    }g1,g2;
    
    void tarjan(int x,int fa){
    	dfn[x]=low[x]=++tim,sta[++top]=x;
    	for (int y=g1.now[x];y;y=g1.pre[y]) if (!vis[y]){
    		int v=g1.son[y];vis[y]=vis[y^1]=1;
    		if (!dfn[v]) tarjan(v,x),low[x]=min(low[x],low[v]);
    		else low[x]=min(low[x],dfn[v]);
    	}
    	if (low[x]==dfn[x]&&fa) top--,g2.ins(x,fa),bnm[x]++,bnm[fa]++;
    	if (low[x]==dfn[fa]){
    		int xx;bcnt++;
    		bel[fa].push_back(bcnt),scc[bcnt].push_back(fa);
    		do{xx=sta[top--],bel[xx].push_back(bcnt),scc[bcnt].push_back(xx);}while (x!=xx);
    	}
    }
    
    void rebuild(){
    	/*puts("");for (int i=1;i<=bcnt;i++){printf("i%d ",i);
    		for (uint j=0;j<scc[i].size();j++) printf("%d ",scc[i][j]);puts("");
    	}
    	puts("fuckpp
    
    
    
    ");*/
    	for (int i=1;i<=n;i++)
    		if (bnm[i]+bel[i].size()>=2)
    			for (uint j=0;j<bel[i].size();j++)
    				g2.ins(i,bel[i][j]+n);//,printf("%d %d
    ",i,bel[i][j]+n)
    }
    
    int main(){
    	scanf("%d%d",&n,&m),g1.tot=1;
    	for (int i=1,x,y,k;i<=m;i++){
    		scanf("%d%d",&k,&x);
    		for (int j=1;j<k;j++) scanf("%d",&y),g1.ins(x,y),x=y;
    	}
    	tarjan(1,0),rebuild(),g2.tree_dp(bcnt?n+1:1,0),printf("%d
    ",ans);
    	return 0;
    }
    </span>

  • 相关阅读:
    启动 YARN 并运行 MapReduce 程序(伪分布式模式)
    启动 HDFS 并运行 MapReduce 程序(伪分布式模式)
    简单计算器(stack)
    Linux定时发邮件脚本
    HttpClient接口调用-客户端
    获取时间字符串
    Visual Assist代码高亮突然失效
    批量快速生成员工文件夹工具
    日语学习笔记整理(汉译日)
    有关使用PLSQL Developer时出现报错ora-12514解决的方法
  • 原文地址:https://www.cnblogs.com/thythy/p/5493471.html
Copyright © 2020-2023  润新知