wannalfy 挑战赛8 F 白云的树(树形dp)
链接:https://www.nowcoder.com/acm/contest/57/F
时间限制:C/C++ 2秒,其他语言4秒
64bit IO Format: %lld
题目描
白云会定期对树做一些修改。
输入描述:
第一行1个整数n,Q,表示树的结点个数和事件个数。
第二行n个整数表示val
1 ... n
。
第三行n-1个整数表示2...n号结点的父亲。
输出描述:
对于op=1,每行一个数表示答案。答案对1e9+7
取模。
输入例子:
15 15
6 4 8 6 8 9 10 9 2 9 8 3 3 6 2
1 1 1 3 1 3 6 5 3 10 9 7 9 13
1 13 8
1 2 4
1 12 8
1 11 2
0 2 9
0 11 4
1 3 6
0 4 5
1 13 7
1 8 2
1 7 6
1 7 8
0 5 5
0 1 7
1 1 6
输出例子:
128592000
11304
35520768
72
10402608
16325280
81
6030720
359079840
8686888
-->
输入
15 15
6 4 8 6 8 9 10 9 2 9 8 3 3 6 2
1 1 1 3 1 3 6 5 3 10 9 7 9 13
1 13 8
1 2 4
1 12 8
1 11 2
0 2 9
0 11 4
1 3 6
0 4 5
1 13 7
1 8 2
1 7 6
1 7 8
0 5 5
0 1 7
1 1 6
输出
128592000
11304
35520768
72
10402608
16325280
81
6030720
359079840
8686888
备注:
//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
#include <bits/stdc++.h>
#define mst(a,b) memset((a),(b), sizeof a)
#define lowbit(a) ((a)&(-a))
#define IOS ios::sync_with_stdio(0);cin.tie(0);
#define MP make_pair
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int mod=1e9+;
const int maxn=1e5+;
ll dp[maxn][];
int val[maxn],fa[maxn];
vector<int>son[maxn];
int qpow(ll a,int b){
ll res=;
while(b){
if(b&)res=res*a%mod;
b>>=;
a=a*a%mod;
}
return res;
}
void add(int a,int b){
for(int i=;i;--i)
for(int j=;j<i;++j)
dp[a][i]=(dp[a][i]+dp[b][j]*dp[a][i-j])%mod; }
void del(int a,int b){
for(int i=;i<=;++i)
for(int j=;j<i;++j)
dp[a][i]=(dp[a][i] - dp[b][j]*dp[a][i-j])%mod;
}
void dfs(int pos){
dp[pos][]=;dp[pos][]=val[pos];
for(int i=;i<son[pos].size();++i){
int to=son[pos][i];
dfs(to);
add(pos,to);
}
}
void update_1(int pos){
if(pos==)return;
update_1(fa[pos]);
del(fa[pos],pos);
}
void update_2(int pos){
if(pos==)return;
add(fa[pos],pos);
update_2(fa[pos]);
}
void update(int k,int c){
update_1(k);
int cc=qpow(val[k],mod-);
for(int i=;i<=;++i)
dp[k][i]=dp[k][i]*cc%mod*c%mod;
val[k]=c;
update_2(k);
}
ll dd[][];
int now;
void qq(int pos,int s){
if(pos!=)qq(fa[pos],pos);
now^=;
memcpy(dd[now],dp[pos],sizeof dp[pos]);
for(int i=;i<=;++i)
for(int j=;j<i;++j)
dd[now][i]=(dd[now][i] - dp[s][j]*dd[now][i-j])%mod;
if(pos!=){
for(int i=;i;--i)
for(int j=;j<i;++j)
dd[now][i]=(dd[now][i]+dd[now^][j]*dd[now][i-j])%mod;
}
}
int query(int k,int s){
if(k==)return (dp[k][s]+mod)%mod;
int f=fa[k];
qq(f,k);
ll ret=;
for(int i=;i<s;++i)
ret=(ret+dd[now][i]*dp[k][s-i])%mod;
return (ret+mod)%mod;
}
int main(){
#ifdef local
freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
#endif
int n,q;scanf("%d%d",&n,&q);
for(int i=;i<=n;++i)scanf("%d",&val[i]);
for(int i=;i<=n;++i){
scanf("%d",&fa[i]);
son[fa[i]].push_back(i);
}
dfs();
while(q--){
int op,a,b;scanf("%d%d%d",&op,&a,&b);
if(op)printf("%d\n",query(a,b));
else update(a,b);
}
return ;
}
首先是随机树的树高是logn这点是不知道的,所以看到题就是蒙圈的,这点是学到了
然后在码的过程中update那个函数,里面对于k点dp值的更新我一开始是这样的
for(int i=;i<=;++i){
dp[k][i]=(dp[k][i]-val[k]*dp[k][i-])%mod;
dp[k][i]=(dp[k][i]+mod)%mod;
}
val[k]=c;
for(int i=;i;--i){
dp[k][i]=(dp[k][i]+dp[k][i-]*val[k])%mod;
}
错的很离谱,引以为鉴
wannalfy 挑战赛8 F 白云的树(树形dp)的更多相关文章
- 牛客挑战赛30 小G砍树 树形dp
小G砍树 dfs两次, dp出每个点作为最后一个点的方案数. #include<bits/stdc++.h> #define LL long long #define fi first # ...
- wannalfy 挑战赛7 F Masha与老鼠(贪心+dp)
链接:https://www.nowcoder.net/acm/contest/56/F 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524288K 6 ...
- 【BZOJ-3572】世界树 虚树 + 树形DP
3572: [Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1084 Solved: 611[Submit][Status ...
- 【BZOJ-2286】消耗战 虚树 + 树形DP
2286: [Sdoi2011消耗战 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 2120 Solved: 752[Submit][Status] ...
- bzoj 2286(虚树+树形dp) 虚树模板
树链求并又不会写,学了一发虚树,再也不虚啦~ 2286: [Sdoi2011]消耗战 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 5002 Sol ...
- 洛谷 P1453 城市环路 ( 基环树树形dp )
题目链接 题目背景 一座城市,往往会被人们划分为几个区域,例如住宅区.商业区.工业区等等.B市就被分为了以下的两个区域--城市中心和城市郊区.在着这两个区域的中间是一条围绕B市的环路,环路之内便是B市 ...
- BZOJ_2286_[Sdoi2011]消耗战_虚树+树形DP+树剖lca
BZOJ_2286_[Sdoi2011]消耗战_虚树+树形DP Description 在一场战争中,战场由n个岛屿和n-1个桥梁组成,保证每两个岛屿间有且仅有一条路径可达.现在,我军已经侦查到敌军的 ...
- 51nod 1353 树 | 树形DP经典题!
51nod 1353 树 | 树形DP好题! 题面 切断一棵树的任意条边,这棵树会变成一棵森林. 现要求森林中每棵树的节点个数不小于k,求有多少种切法. 数据范围:\(n \le 2000\). 题解 ...
- BZOJ 2286 消耗战 (虚树+树形DP)
给出一个n节点的无向树,每条边都有一个边权,给出m个询问,每个询问询问ki个点,问切掉一些边后使得这些顶点无法与顶点1连接.最少的边权和是多少.(n<=250000,sigma(ki)<= ...
随机推荐
- 使用svn未响应卡死的几个原因,commit时checkout时
1.commit 时 很可能是:检索文件内容过多导致,解决:不要在最外层文件夹目录下commit 2.checkout时 很可能是:地址错误
- MATLAB中产生高斯白噪声的两个函数
MATLAB中产生高斯白噪声非常方便,可以直接应用两个函数,一个是WGN,另一个是AWGN.WGN用于产生高斯白噪声,AWGN则用于在某一信号中加入高斯白噪声.1.WGN:产生高斯白噪声 y = wg ...
- IDeajCommunity 配置smart tomcat插件
ntelliJ IDEA社区版没有自带tomcat. 装插件--smart tomcat. IntelliJ IDEA>>Preferences>>Plugins>> ...
- 用python实现简单的计算器(加减乘除小括号等)
需求:实现能计算类似 1 - 2 * ( (60-30 +(-40/5) * (9-2*5/3 + 7 /3*99/4*2998 +10 * 568/14 )) - (-4*3)/ (16-3*2) ...
- [USACO12FEB]牛券Cow Coupons(堆,贪心)
[USACO12FEB]牛券Cow Coupons(堆,贪心) 题目描述 Farmer John needs new cows! There are N cows for sale (1 <= ...
- 清楚webView的缓存
+(void)clearCache{ NSHTTPCookie *cookie; NSHTTPCookieStorage *storage = [NSHTTPCookieStorage s ...
- python基础练习题4
题目:现有一个数据库记录文件(0005.txt)保证了学生课程签到的数据记录('2017-03-13 11:50:09',271,131),('2017-03-14 11:52:19',273,131 ...
- samba服务和client挂载
服务端 1.安装samba服务 yum -y install samba 2.创建系统用户 因为Samba 服务程序的数据库要求账户必须在当前系统中已经存在,否则日后创建文件时将导致文件的权限属性混乱 ...
- 微信小程序-自制弹出框禁止页面上下滑动
弹出 fixed 弹窗后,在弹窗上滑动会导致下层的页面一起跟着滚动. 解决方法: 在弹出层加上 catchtouchmove 事件 两种方法:(在电脑上测试是没有用的,这是触摸事件.因此,需要在手机端 ...
- h5页面弹窗时页面固定(弹窗下面的页面不滑动)
页面出现弹窗时,底部页面不能随之滑动怎么解决? 只需将页面的body增加一个样式 overflow:hidden;就能解决 jq: //开启弹窗 $('body').attr('style','ove ...