LINK：Decompose
看起来很难 实际上也很难 考验选手的dp 树链剖分 矩阵乘法的能力。
容易列出dp方程 暴力dp 期望得分28.
对于链的情况 容易发现dp方程可以转矩阵乘法 然后利用线段树维护矩阵即可。
这个矩阵很容易列出这里不再赘述。
对于100分 容易想到动态dp模型 LCT写动态dp是万万不能的。
而且这道题的dp方程和其他儿子也有些关系。
考虑树链剖分 然后分别计算轻儿子和重儿子的贡献。
让重儿子利用矩阵来进行转移 轻儿子当做常数.
这样每次修改的时候 修改的节点最多只有logn个.
用set维护需要维护的东西即可。
剩下的就是树链剖分型动态dp的套路 每次利用链顶的矩阵信息更新下一个节点即可。
一个细节：叶子节点可以直接列成(Lcdot L)的矩阵 只有第一个元素有值 这样更容易实现。
一个细节：可能矩阵乘法出来的值和原来的值不尽相同 此时考虑更新set的时候利用原来信息更新 然后更新原来信息即可。
思维难度：高 代码难度：极高。
const ll MAXN=100010;
ll n,Q,L,len,id;
ll top[MAXN],pos[MAXN],dfn[MAXN],fa[MAXN],c[MAXN],sum[MAXN];
ll a[MAXN][5],d[MAXN],son[MAXN],sz[MAXN],f[MAXN][5],w[MAXN];
ll lin[MAXN],ver[MAXN<<1],nex[MAXN<<1];
multiset<ll>s[MAXN][4];
struct wy
{
ll b[5][5];
ll l,r;
wy(){l=r=0;rep(1,L,i)rep(1,L,j)b[i][j]=-INF;}
wy friend operator +(wy a,wy b)
{
wy c;c.l=b.l;c.r=a.r;
rep(1,L,i)rep(1,L,j)rep(1,L,k)
c.b[i][j]=max(c.b[i][j],a.b[i][k]+b.b[k][j]);
return c;
}
}t[MAXN<<2];
inline void add(ll x,ll y)
{
ver[++len]=y;
nex[len]=lin[x];
lin[x]=len;
}
inline void dfs(ll x,ll father)
{
d[x]=d[father]+1;fa[x]=father;sz[x]=1;
rep(2,L,j)f[x][j]=-INF;w[x]=-INF;
f[x][1]=a[x][1];ll ans=0;
go(x)if(tn!=father)
{
dfs(tn,x);
if(sz[tn]>sz[son[x]])son[x]=tn;
sz[x]+=sz[tn];
rep(2,L,j)f[x][j]=max(f[x][j]+w[tn],f[tn][j-1]+a[x][j]+ans);
ans+=w[tn];
}
f[x][1]+=ans;
rep(1,L,j)w[x]=max(w[x],f[x][j]);
go(x)if(tn!=father&&tn!=son[x])
rep(1,L-1,j)s[x][j].insert(f[tn][j]-w[tn]);
sum[x]=ans-w[son[x]];
}
inline void dp(ll x,ll father)
{
top[x]=father;dfn[x]=++id;pos[id]=x;c[father]=x;
if(!son[x])return;
dp(son[x],father);
go(x)if(tn!=son[x]&&tn!=fa[x])dp(tn,tn);
}
inline void build(ll p,ll l,ll r)
{
l(p)=l;r(p)=r;
if(l==r)
{
ll x=pos[l];
if(sz[x]==1)
{
rep(1,L,i)rep(1,L,j)t[p].b[i][j]=-INF;
t[p].b[1][1]=a[x][1];
}
else
{
rep(1,L,i)t[p].b[i][1]=a[x][1]+sum[x];
ll flag=s[x][1].size();
rep(2,L,j)
{
ll ww=flag?(*--s[x][j-1].end()):-INF;
rep(1,L,i)t[p].b[i][j]=a[x][j]+sum[x]+ww;
t[p].b[j-1][j]-=ww;
}
}
return;
}
ll mid=(l+r)>>1;
build(zz,l,mid);
build(yy,mid+1,r);
t[p]=t[yy]+t[zz];
}
inline void change(ll p,ll x)
{
if(l(p)==r(p))
{
ll x=pos[l(p)];
if(sz[x]==1)
{
rep(1,L,i)rep(1,L,j)t[p].b[i][j]=-INF;
t[p].b[1][1]=a[x][1];
}
else
{
rep(1,L,i)t[p].b[i][1]=a[x][1]+sum[x];
ll flag=s[x][1].size();
rep(2,L,j)
{
//cout<<s[x][j-1].size()<<endl;
ll ww=flag?(*(--s[x][j-1].end())):-INF;
rep(1,L,i)t[p].b[i][j]=a[x][j]+sum[x]+ww;
t[p].b[j-1][j]-=ww;
//cout<<t[p].b[j-1][j]<<endl;
}
//rep(1,L,i){rep(1,L,j)cout<<t[p].b[i][j]<<' ';cout<<endl;}
}
return;
}
ll mid=(l(p)+r(p))>>1;
if(x<=mid)change(zz,x);
else change(yy,x);
t[p]=t[yy]+t[zz];
}
inline wy ask(ll p,ll l,ll r)
{
if(l<=l(p)&&r>=r(p))return t[p];
ll mid=(l(p)+r(p))>>1;
if(l>mid)return ask(yy,l,r);
if(r<=mid)return ask(zz,l,r);
return ask(yy,l,r)+ask(zz,l,r);
}
inline void Tchange(ll x)
{
ll fx=top[x];
while(fx!=1)
{
change(1,dfn[x]);
wy w1=ask(1,dfn[fx],dfn[c[fx]]);
x=fa[fx];//修改x.
ll cnt1=-INF;
rep(1,L,i)cnt1=max(cnt1,w1.b[1][i]);
sum[x]=sum[x]-w[fx]+cnt1;
fep(L-1,1,i)
{
s[x][i].erase(s[x][i].find(f[fx][i]-w[fx]));
s[x][i].insert(w1.b[1][i]-cnt1);
f[fx][i]=w1.b[1][i];
}
w[fx]=cnt1;fx=top[x];
}
change(1,dfn[x]);
wy ww=ask(1,dfn[1],dfn[c[1]]);ll ans=-INF;
rep(1,L,i)
{
ans=max(ans,ww.b[1][i]);
//cout<<ww.b[1][i]<<' ';
}
//puts("");
putl(ans);
}
signed main()
{
freopen("decompose.in","r",stdin);
freopen("decompose.out","w",stdout);
get(n);get(Q);get(L);
rep(2,n,i)add(read(),i);
rep(1,n,i)rep(1,L,j)get(a[i][j]);
dfs(1,0);dp(1,1);
build(1,1,n);
//wy ww=ask(1,dfn[1],dfn[c[1]]);ll ans=-INF;
//rep(1,L,i)ans=max(ans,ww.b[1][i]),cout<<ww.b[1][i]<<' ';
//puts("");putl(ans);
rep(1,Q,i)
{
ll get(x);
rep(1,L,j)get(a[x][j]);
Tchange(x);
}
return 0;
}