Equation

题目描述

有一棵\(n\)个点的以\(1\)为根的树,以及\(n\)个整数变量\(x_i\)。树上\(i\)的父亲是\(f_i\),每条边\((i,f_i)\)有一个权值\(w_i\),表示一个方程\(x_i+x_{f_i}=w_i\),这\(n-1\)个方程构成了一个方程组。

现在给出\(q\)个操作,有两种类型:

•\(1\ u\ v\ s\),表示询问加上\(x_u+x_v=s\)这个方程后,整个方程组的解的情况。具体来说,如果方程有唯一解,输出此时\(x_1\)的值;如果有无限多个解,输出inf;如果无解,输出none. 注意每个询问是独立的.

•\(2\ u\ w\),表示将\(w_u\)修改为\(w\).

输入输出格式

输入格式

从文件equation.in中读入数据.

第一行两个整数\(n,q\)。

接下来\(n-1\)行,第\(i\)行有两个整数\(f_{i+1}\)和\(w_{i+1}\)。

接下来\(q\)行,每行表示一个操作,格式见问题描述。

输出格式

输出到文件equation.out中.

对于每个询问输出一行表示答案.

说明

对于所有数据,有\(1\le n,q\le 10^6,1\le f_i\le i -1,1\le u,v\le n; -10^3\le w,w_i\le 10^3,-10^9\le s\le 10^9\).

• \(\text{Subtask1}(3\%),n\le 10,q=0\).

• \(\text{Subtask2}(18\%), n=2\).

• \(\text{Subtask3}(32\%), n,q\le 10^3\).

• \(\text{Subtask4}(33\%), n,q\le 10^5\).

• \(\text{Subtask5}(14\%)\),没有特殊的约束.


给个1e5,给个1e6真的坑。

1e6就认为是\(O(n)\)做法好吗(反正我大常数这辈子过不了1e6的\(nlogn\)了

不过为了卡两个\(log\)的还可以理解了


发现连上一条边后,把路径的边权正负相加可以消去很多项。

当路径长度为奇数,判断无解or多解

偶数,解出来判断是否有整数解

然后求上去就行了

发现我们要维护路径求和和单点修改

可以无脑树剖但是不能拿满分

路径求和可以根据lca分成两个分别求

维护一个到根节点的路径长度

跑一下dfs序,就变成了维护区间求和和单点加,树状数组就可以了

然而我常数真的大。。

细节处理起来听麻烦的


Code:

#include <cstdio>
#define ll long long
const int N=1e6+10;
int read()
{
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9') {if(c=='-') f=-f;c=getchar();}
while(c>='0'&&c<='9') {x=x*10+c-'0';c=getchar();}
return x*f;
}
int n,q;
int head[N],to[N],Next[N],cnt;
void add(int u,int v)
{
to[++cnt]=v,Next[cnt]=head[u],head[u]=cnt;
}
int dep[N],f[N][21],w[N],dfn[N],siz[N],dfsclock,tmp;
ll s[N];
void swap(int &x,int &y){tmp=x,x=y,y=tmp;}
void modify(int x,ll de)
{
for(int i=x;i<=n;i+=i&-i)
s[i]+=de;
}
ll query(int x)
{
ll sum=0;
for(int i=x;i;i-=i&-i)
sum+=s[i];
return sum;
}
void dfs(int now)
{
dfn[now]=++dfsclock,siz[now]=1;
modify(dfn[now],w[now]*(dep[now]&1?1:-1));
for(int i=1;f[now][i-1];i++) f[now][i]=f[f[now][i-1]][i-1];
for(int i=head[now];i;i=Next[i])
dep[to[i]]=dep[now]+1,dfs(to[i]),siz[now]+=siz[to[i]];
modify(dfn[now]+siz[now],w[now]*(dep[now]&1?-1:1));
}
int LCA(int u,int v)
{
for(int i=20;~i;i--)
if(dep[f[u][i]]>=dep[v])
u=f[u][i];
if(u==v) return u;
for(int i=20;~i;i--)
if(f[u][i]!=f[v][i])
u=f[u][i],v=f[v][i];
return f[u][0];
}
int main()
{
//freopen("data.in","r",stdin);
//freopen("data.out","w",stdout);
n=read(),q=read();
for(int i=2;i<=n;i++)
f[i][0]=read(),w[i]=read(),add(f[i][0],i);
dfs(1);
for(int p,op,u,v,s,i=1;i<=q;i++)
{
op=read();
if(op==1)
{
u=read(),v=read(),s=read();
if(dep[u]<dep[v]) swap(u,v);
int lca=LCA(u,v);p=u;
if(dep[u]-dep[v]&1)
{
ll k=(query(dfn[u])-query(dfn[lca]))*(dep[u]&1?1:-1);
k+=(query(dfn[v])-query(dfn[lca]))*(dep[v]&1?1:-1);
if(s==k) printf("inf\n");
else printf("none\n");
}
else
{
ll k=1ll*(query(dfn[u])-query(dfn[lca]))*(dep[u]&1?1:-1);
k+=1ll*(query(dfn[v])-query(dfn[lca]))*(dep[v]&1?-1:1);
k=k+s;
if(k&1) printf("none\n");
else
{
k>>=1;
printf("%lld\n",query(dfn[p])-k*(dep[p]&1?1:-1));
}
}
}
else
{
u=read();
p=read();
modify(dfn[u],(p-w[u])*(dep[u]&1?1:-1));
modify(dfn[u]+siz[u],(w[u]-p)*(dep[u]&1?1:-1));
w[u]=p;
}
}
return 0;
}

2018.10.6

雅礼集训 Day6 T2 Equation 解题报告的更多相关文章

  1. 雅礼集训 Day7 T1 Equation 解题报告

    Reverse 题目背景 小\(\text{G}\)有一个长度为\(n\)的\(01\)串\(T\),其中只有\(T_S=1\),其余位置都是\(0\).现在小\(\text{G}\)可以进行若干次以 ...

  2. 雅礼集训 Day6 T1 Merchant 解题报告

    Merchant 题目描述 有\(n\)个物品,第\(i\)个物品有两个属性\(k_i,b_i\),表示它在时刻\(x\)的价值为\(k_i\times x+b_i\). 当前处于时刻\(0\),你可 ...

  3. 雅礼集训 Day3 T2 u 解题报告

    u 题目背景 \(\frac 14\) 遇到了一道水题,完全不会做,于是去请教小\(\text{D}\).小\(\text{D}\)看了一眼就切掉了这题,嘲讽了\(\frac 14\)一番就离开了. ...

  4. 雅礼集训 Day3 T2 v 解题报告

    v 题目背景 \(\frac 14\)遇到了一道水题,又完全不会做,于是去请教小\(\text{D}\).小\(\text{D}\)看了\(0.607\)眼就切掉了这题,嘲讽了\(\frac 14\) ...

  5. 「雅礼集训 2017 Day2」解题报告

    「雅礼集训 2017 Day2」水箱 我怎么知道这种题目都能构造树形结构. 根据高度构造一棵树,在树上倍增找到最大的小于约束条件高度的隔板,开一个 \(vector\) 记录一下,然后对于每个 \(v ...

  6. 「雅礼集训 2017 Day1」 解题报告

    「雅礼集训 2017 Day1」市场 挺神仙的一题.涉及区间加.区间除.区间最小值和区间和.虽然标算就是暴力,但是复杂度是有保证的. 我们知道如果线段树上的一个结点,\(max=min\) 或者 \( ...

  7. 雅礼集训 Day3 T3 w 解题报告

    w 题目背景 \(\frac 14\)遇到了一道水题,双完全不会做,于是去请教小\(\text{D}\).小\(\text{D}\)看了\(0.607^2\)眼就切掉了这题,嘲讽了\(\frac 14 ...

  8. 雅礼集训 Day1 T3 画作 解题报告

    画作 题目描述 小\(\mathrm{G}\)的喜欢作画,尤其喜欢仅使用黑白两色作画. 画作可以抽象成一个\(r\times c\)大小的\(01\)矩阵.现在小\(\mathrm{G}\)构思好了他 ...

  9. 雅礼集训 Day5 T3 题 解题报告

    题 题目背景 由于出题人赶时间所以没办法编故事来作为背景. 题目描述 一开始有\(n\)个苹果,\(m\)个人依次来吃苹果,第\(i\)个人会尝试吃\(u_i\)或\(v_i\)号苹果,具体来说分三种 ...

随机推荐

  1. noip2018 洛谷 P5020 货币系统

    关键: 要使m最小,(m,b)中的数不能用(n,a)中的数表示出来 对于 3  19  10  6 19=10+3+3+3 6=3+3 只有3 和 10 不能被(n,a)中的数表示 所以m=2 只需要 ...

  2. MySql错误1045 Access denied for user 'root'@'localhost' (using password:YES)

    1.先停止mysql服务 2.进入mysql的安装路径,找到并打开my.ini文件,找到[mysqld],在该行下面添加 skip_grant_tables,也就是通知mysql,在登陆的时候跳过密码 ...

  3. oracle之bitmap索引

    oracle常见的索引是BTree索引和Bitmap索引. BTree索引特点: 默认索引 适合大量增删改查 不能用or操作符 适合高基数的列(即唯一值多) 创建sql:create index li ...

  4. java 反射bean

    描述:新bean和旧bean拥有相同的字段属性.将旧bean的属性值,赋给新bean. Java代码如下: /** * 反射Bean * @param bean 反射的对象 * @param newB ...

  5. Shell学习——终端打印

    1.echo1.1.默认情况下,echo在每次调用后会添加一个换行符1.2.待打印的内容,可以用单引号.双引号或者直接打印,不同的方式,有各自的限制1.2.1.使用不带引号的echo时,没法打印分好( ...

  6. Python 对象(type/object/class) 作用域 一等函数 (慕课--Python高级,IO并发 第二章)

    在python中一共有两种作用域:全局作用域和函数作用域全局作用域:在全局都有效,全局作用域在程序执行时创建,在程序执行结束时销毁:所有函数以外的区域都是全局作用域:在全局作用域中定义的变量,都属于全 ...

  7. keil5的安装及问题

    win8+keil 注意,在进行破解的时候首先要打开一个工程,而且keil要用管理员的身份进行运行, 才可以破解完成 发现打开之后,出现这样的错误. 原因是因为在创建工程的时候在下图中点了是,要点否才 ...

  8. LoadRunner使用代理远程执行提示找不到“pre_cci.c”文件

    好久没有使用LoadRunner了,工作需要使用一下,执行总是提示找不到“pre_cci.c”文件,找问题花了很长时间终于找到问题了.万事还是需要找到原因: cci 会将 pre_cci.c 文件作为 ...

  9. PCB工艺要求

    项目 加工能力 工艺详解   层数 1~6层 层数,是指PCB中的电气层数(敷铜层数).目前嘉立创只接受1~6层板.   板材类型 FR-4板材 板材类型:纸板.半玻纤.全玻纤(FR-4).铝基板,目 ...

  10. Java密码学综述---密码学基本功能

    机密性 只有发送方与指定接收方才能理解报文的内容,监听者只能截取到加密后的报文信息,但不能还原报文的内容,做到了保密 鉴别 发送方与接收方都应该能证实通信过程所涉及的另一方,通信的另一方确实具有它们所 ...