bzoj 1093 缩点+DP
首先比较明显的是如果存在一个半连通子图,我们将其中的环缩成点,那么该图仍为半连通子图,这样我们就可以先将整张图缩点,重新构图,新图为拓扑图,记录每个新的点表示的强连通分量中点的个数num[i],那么我们就可以DP了,新图中的每一条链都为原图的半连通子图,这样我们找到新图中的最长链就行了,找入度为0的点dfs做树上DP,这样我们可以知道每个点的len[i]代表从这个点开始的最长链的长度,len[i]=max(len[son of i])+num[i],然后我们求出来了第一问,对于第二问,我们需要找len[i]=ans1的点做dfs,然后设ans[i]为以i为根的子树的方案数,那么ans[i]+=ans[son of i] (len[i]=len[son of i]+num[i]),因为我们需要找最长链上的点来更新答案,这样最后再累计答案就好了。
反思:开始题中说没有重边,但是没有考虑到重构图之后的图是可能有重边的,这样第二问的答案就可能会被重复累加,所以我们DP的时候可以维护一个栈,和每个栈中元素的父亲,这样对于一个点枚举子节点,如果子节点没有在栈中出现过,那么就累加答案,该子节点进栈,dfs最后的时候再弹出所有栈中x的子节点。
/**************************************************************
Problem: 1093
User: BLADEVIL
Language: C++
Result: Accepted
Time:1992 ms
Memory:45028 kb
****************************************************************/
//By BLADEVIL
#include <cstdio>
#include <algorithm>
#define maxn 200020
#define maxm 4000040
using namespace std;
int n,m,d39,l,tot,time,size;
int last[maxn],other[maxm],pre[maxm],stack[maxn],low[maxn],dfn[maxn],flag[maxn],col[maxn],num[maxn],in[maxn];
int len[maxn],ans[maxn],father[maxn];
int p1,p2;
void connect(int x,int y){
pre[++l]=last[x];
last[x]=l;
other[l]=y;
//if (x>n||y>n) printf("%d %d\n",x,y);
}
void tarjan(int x){
//printf("%d %d\n",x,fa);
low[x]=dfn[x]=++time;
stack[++tot]=flag[x]=x;
//for (int i=1;i<=tot;i++) printf("%d ",stack[i]); printf("\n");
for (int p=last[x];p;p=pre[p]){
if (!dfn[other[p]]) tarjan(other[p]),low[x]=min(low[x],low[other[p]]); else
if (flag[other[p]]) low[x]=min(low[x],dfn[other[p]]);
}
if (low[x]==dfn[x]){
int cur=-;
while (cur!=x){
cur=stack[tot--];
flag[cur]=;
col[cur]=size;
num[size]++;
}
size++;
}
}
void dfs(int x){
int cur=;
for (int p=last[x];p;p=pre[p]){
if (!len[other[p]]) dfs(other[p]);
cur=max(cur,len[other[p]]);
}
len[x]=cur+num[x];
//printf(" %d %d\n",x,len[x]);
}
void work(int x){
for (int p=last[x];p;p=pre[p])
if (len[other[p]]+num[x]==len[x]) {
if (flag[other[p]]) continue;
if (!ans[other[p]]) work(other[p]);
ans[x]+=ans[other[p]];
stack[++tot]=other[p]; flag[other[p]]=; father[other[p]]=x;
}
if (!ans[x]) ans[x]=;
ans[x]%=d39;
while (father[stack[tot]]==x) flag[stack[tot--]]=;
}
int main(){
int x,y;
scanf("%d%d%d",&n,&m,&d39); size=n+;
for (int i=;i<=m;i++) scanf("%d%d",&x,&y),connect(x,y);
for (int i=;i<=n;i++) if (!low[i]) tarjan(i);
//for (int i=1;i<=n;i++) printf("%d %d %d\n",col[i],low[i],dfn[i]);
for (int i=;i<=n;i++)
for (int p=last[i];p;p=pre[p])
if (col[i]!=col[other[p]]) connect(col[i],col[other[p]]),in[col[other[p]]]++;
//for (int i=n+1;i<size;i++) printf("|%d %d\n",i,num[i]);
for (int i=n+;i<size;i++) if (!in[i]) dfs(i);
//for (int i=n+1;i<size;i++) printf("|%d %d\n",i,len[i]);
for (int i=n+;i<size;i++) p1=max(p1,len[i]);
for (int i=n+;i<size;i++) if (len[i]==p1) work(i);
for (int i=n+;i<size;i++) if (len[i]==p1) (p2+=ans[i])%=d39;
//for (int i=n+1;i<size;i++) printf("|%d %d\n",i,ans[i]);
printf("%d\n%d\n",p1,p2);
return ;
}
bzoj 1093 缩点+DP的更多相关文章
- BZOJ 1093 [ZJOI2007] 最大半连通子图(强联通缩点+DP)
题目大意 题目是图片形式的,就简要说下题意算了 一个有向图 G=(V, E) 称为半连通的(Semi-Connected),如果满足图中任意两点 u v,存在一条从 u 到 v 的路径或者从 v 到 ...
- bzoj 1093 [ZJOI2007]最大半连通子图——缩点+拓扑
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1093 缩点+拓扑,更新长度的时候维护方案数. 结果没想到处理缩点后的重边,这样的话方案数会算 ...
- BZOJ 1093 [ZJOI2007]最大半连通子图
1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 1986 Solved: 802[Submit][St ...
- POJ3160 Father Christmas flymouse[强连通分量 缩点 DP]
Father Christmas flymouse Time Limit: 1000MS Memory Limit: 131072K Total Submissions: 3241 Accep ...
- [bzoj 1093][ZJOI2007]最大半联通子图(强联通缩点+DP)
题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1093 分析: 首先肯定是先把强联通全部缩成一个点,然后成了一个DAG 下面要知道一点: ...
- BZOJ 1093 强连通缩点+DAG拓扑DP
缩点后在一个DAG上求最长点权链 和方案数 注意转移条件和转移状态 if (nowmaxn[x] > nowmaxn[v]) { ans[v] = ans[x]; nowmaxn[v] = no ...
- BZOJ 1093: [ZJOI2007]最大半连通子图( tarjan + dp )
WA了好多次... 先tarjan缩点, 然后题意就是求DAG上的一条最长链. dp(u) = max{dp(v)} + totu, edge(u,v)存在. totu是scc(u)的结点数. 其实就 ...
- bzoj 1093 [ ZJOI 2007 ] 最大半连通子图 —— 拓扑+DP
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1093 先缩点,然后就是找最长链,DP一下即可: 注意缩点后的重边!会导致重复计算答案. 代码 ...
- bzoj 1093 [ZJOI2007]最大半连通子图(scc+DP)
1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 2286 Solved: 897[Submit][St ...
随机推荐
- ubuntu 安装xdebug
Add XDebug to Ubuntu 14.04 Submitted by Wilbur on Tue, 06/17/2014 - 12:49pm It's pretty easy to add ...
- try catch finally 与continue的使用
using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.T ...
- MySQL中的条件语句
判断学生表中成绩是否小于60,将小于60的学生成绩列为不及格 学生表(student) 字段:姓名(name),学号(主键)(num),性别(sex),成绩(score) select *,if(sc ...
- Spring Bean注册和加载
Spring解密 - XML解析 与 Bean注册 Spring解密 - 默认标签的解析 Spring解密 - 自定义标签与解析 Spring解密 - Bean的加载流程
- java session特性
1.当前浏览器不关闭 则一直有效 servlet就能取到值(未设置过期时间情况下 或者在过期的时间范围内) 算成一次会话 再次会话内多个请求都能获得session 2.session保存在服务端,通 ...
- BZOJ4835 遗忘之树
点分树上的某个点和其某个子树在原树中的连接方式一般来说可以是由该点连向子树内任意一点,这样方案数即为所有子树大小之积.但有一种特殊情况是连接某点后导致编号最小的重心更换,只要去掉这种就行了,具体地可以 ...
- vue-cli项目打包出现空白页和路径错误问题
vue-cli项目打包: 1. 命令行输入:npm run build 打包出来后项目中就会多了一个文件夹dist,这就是我们打包过后的项目. 第一个问题,文件引用路径.我们直接运行打包后的文件夹 ...
- IO模式
二 IO模式 刚才说了,对于一次IO访问(以read举例),数据会先被拷贝到操作系统内核的缓冲区中,然后才会从操作系统内核的缓冲区拷贝到应用程序的地址空间.所以说,当一个read操作发生时,它会经历两 ...
- Mysql 基本语句练习
一.怎样查看数据库信息? desc 数据库名; 二.怎样查看数据表信息? desc 表名: //查看表的属性和属性值 或者用select语句: //查看表的行记录信息 select ...
- bzoj 3275: Number (最小割)
题目的意思是要选一些数,但是这些数如果满足两个条件的话就不能一起被选. type arr=record toward,next,cap:longint; end; const maxn=; maxm= ...