食物链 bzoj-4562 Haoi-2016

题目大意:给你n个点,m条边的DAG,求所有的满足条件的链,使得每条链的起点是一个入度为0的点,中点是一条出度为0的点。

注释:$1\le n\le 10^5$,$1\le m\le 2*10^5$。

想法:考试T2,全场切

动态规划

状态:dp[i]表示从这个点到出度为0的点的方案数。

转移:dp[i]+=dp[to[i]]

然后用记忆化爆搜即可。

最后,附上丑陋的代码... ...

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <map>
#define MOD 1000000007
#define N 100010
#define M 200010
using namespace std;
typedef long long ll;
ll dp[N];
int num[N];
char s[N][100];
map<string,int>mp;
// ll stack[M<<1];
// ll val(int a)
// {
// ll ans1=1,ans2=1,ans3=1;
// int k=strlen(s[a]+1);
// for(int i=1;i<=k/3;i++) ans1=ans1*(s[i]-'0')%mod;
// for(int i=k/3+1;i<=k/3*2;i++) ans2=ans2*(s[i]-'0')%mod;
// for(int i=k/3*2+1;i<=k;i++) ans3=ans3*(s[i]-'0')%mod;
// }
struct Node
{
int val,id,now;
}qqq[N];
// int map[1000000];
int head[N],to[M],nxt[M],tot;
inline void add(int x,int y)
{
to[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
int dfs(int pos,int fa)
{
if(dp[pos]) return dp[pos];
for(int i=head[pos];i;i=nxt[i])
{
if(to[i]==fa) continue;
dp[pos]+=dfs(to[i],pos);
dp[pos]%=MOD;
}
if(dp[pos]==0) dp[pos]=1;
return dp[pos]%MOD;
}
inline bool cmp1(Node o,Node oo)
{
return o.val<oo.val;
}
inline bool cmp2(Node o,Node oo)
{
return o.id<oo.id;
}
ll sum[N];
int main()
{
// freopen("chain.in","r",stdin);
// freopen("chain.out","w",stdout);
ll ans=0,/* top=0, */cnt=0;
int m;
scanf("%*d%d",&m);
// puts("Fuck");
for(int i=1;i<=m;i++)
{
string s1,s2;
cin >> s1 >> s2 ;
if(!mp[s1]) mp[s1]=++cnt;
if(!mp[s2]) mp[s2]=++cnt;
add(mp[s1],mp[s2]);
num[mp[s1]]++;
sum[mp[s2]]++;
// scanf("%s%s",s[1]+1,s[2]+1);
// int x=val(1),y=val(2);
// printf("Bitch %d %d\n",x,y);
// top+=2;
// stack[++top]=x,stack[++top]=y;
// qqq[top-1].val=x,qqq[top].val=y;
// qqq[top-1].id=top-1,qqq[top].id=top;
// add(x,y);
}
// sort(qqq+1,qqq+top+1,cmp1);
// qqq[0].val=qqq[1].val-1;
// for(int i=1;i<=top;i++)
// {
// if(qqq[i-1].val!=qqq[i].val) cnt++;
// qqq[i].now=cnt;
// }
// for(int i=1;i<=top;i++)
// {
// printf("Woc %d %d\n",qqq[i].val,qqq[i].now);
// }
// sort(qqq+1,qqq+top+1,cmp2);
// for(int i=1;i<=top;i++)
// {
// printf("Fuck %d\n",qqq[i].now);
// }
// for(int i=2;i<=top;i+=2)
// {
// add(qqq[i-1].now,qqq[i].now);
// sum[qqq[i].now]++;
// maxnow=max(maxnow,sum[qqq[i].now]);
// maxnow=max(maxno)
// }
for(int i=1;i<=cnt;i++)
{
if(!sum[i]&&!num[i]) continue;
if(!sum[i]) ans+=dfs(i,0);
// printf("CaoNiMa %d\n",i);
}
// for(int i=1;i<=cnt;i++)
// {
// printf("EYWR %lld\n",dp[i]);
// }
printf("%lld\n",ans);
return 0;
}

小结:先想到dp是关键。学长写的全是什么toposort..都是O(n)的...

[bzoj4562][Haoi2016]食物链_记忆化搜索_动态规划的更多相关文章

  1. P3183 [HAOI2016]食物链[拓扑/记忆化搜索]

    题目来源:洛谷 题目描述 如图所示为某生态系统的食物网示意图,据图回答第1小题现在给你n个物种和m条能量流动关系,求其中的食物链条数.物种的名称为从1到n编号M条能量流动关系形如a1 b1a2 b2a ...

  2. CF894E Ralph and Mushrooms_强连通分量_记忆化搜索_缩点

    Code: #include<cstdio> #include<stack> #include<cstring> using namespace std; cons ...

  3. 【BZOJ1048】分割矩阵(记忆化搜索,动态规划)

    [BZOJ1048]分割矩阵(记忆化搜索,动态规划) 题面 BZOJ 洛谷 题解 一个很简单的\(dp\),写成记忆化搜索的形式的挺不错的. #include<iostream> #inc ...

  4. [ACM_动态规划] 数字三角形(数塔)_递推_记忆化搜索

    1.直接用递归函数计算状态转移方程,效率十分低下,可以考虑用递推方法,其实就是“正着推导,逆着计算” #include<iostream> #include<algorithm> ...

  5. poj 1088 滑雪_记忆化搜索

    题意:略 直接用记忆化搜索就行了 #include<cstdio> #include<iostream> using namespace std; int n,m; int m ...

  6. 洛谷 p2618 数字工程 记忆化搜索_ 线性筛

    我们在线筛的同时处理出每个数的所有质因子,记忆化搜索的时候直接枚举质因子即可. 时间复杂度为 O(nlogn)O(nlogn)O(nlogn) Code: #include<cstdio> ...

  7. BZOJ4513 SDOI2016 储能表 记忆化搜索(动态规划)

    题意: 题面中文,不予翻译:SDOI2016储能表 分析: 据说有大爷用一些奇怪的方法切掉了这道题%%%%% 这里用的是大众方法——动态规划. 其实这是一道类似于二进制数位dp的动态规划题,(但是实际 ...

  8. 阿牛的EOF牛肉串-记忆化搜索或动态规划

    C - 阿牛的EOF牛肉串 Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u Submi ...

  9. UVa 10599【lis dp,记忆化搜索】

    UVa 10599 题意: 给出r*c的网格,其中有些格子里面有垃圾,机器人从左上角移动到右下角,只能向右或向下移动.问机器人能清扫最多多少个含有垃圾的格子,有多少中方案,输出其中一种方案的格子编号. ...

随机推荐

  1. Webkit内核开源爬虫蜘蛛引擎

    C#开发的基于Webkit内核开源爬虫蜘蛛引擎 https://www.cnblogs.com/micro-chen/p/9075590.html 概述 在各个电商平台发展日渐成熟的今天.很多时候,我 ...

  2. Java执行定时任务

    一.用java.util.Timer 使用JAVA类Timer可实现简单的延迟和周期性任务,其中的任务使用java.util.TimerTask表示.任务的执行方式有两种: 按固定速率执行:即sche ...

  3. Flink之DataStreamAPI入门

    目录 Types Transformations Defining UDFs 本文API基于Flink 1.4 def main(args: Array[String]) { // 第一种会自动判断用 ...

  4. idea UML快捷键

  5. Oracle配置说明

    当Oracle安装完成后,为后续能够顺利得导出空表,特做一下配置(重点关注2.1) 1.1.查询空表select table_name from user_tables where NUM_ROWS= ...

  6. TCP/IP详解(三)

    超时与重传: TCP在发送一个包时,启动一个定时器,如果在定时器溢出之前没有收到ACK,则认为发出的包丢失了,此时会重传丢失的包.这就是超时重传. 其中定时器的时间不是一个固定值,它是根据RTT计算的 ...

  7. POJ 1010 题目翻译+题解

    题目实在是太难懂了,我翻译了一下... STAMPS Description Have you done any Philately lately? 你最近有没有集邮? You have been h ...

  8. POJ 2342 Anniversiry Party(TYVJ1052 没有上司的舞会)

    题意: P1052 没有上司的舞会 描述 Ural大学有N个职员,编号为1~N.他们有从属关系,也就是说他们的关系就像一棵以校长为根的树,父结点就是子结点的直接上司.每个职员有一个快乐指数.现在有个周 ...

  9. Spring的AOP机制---- AOP的注解配置---- AOP的注解配置

    3333隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约隐隐约约噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢噢 ...

  10. .Net Core添加分布式Session

    一.Session HTTP是一个无状态协议,Web服务器将每一个请求都视为独立请求.并且不保存之前请求中用户的值. Session 状态是ASP.NET Core提供的一个功能,它可以在用户通应用访 ...