目录

20190810 西比拉先知系统——sibyl

20190810 西比拉先知系统——sibyl

思路:

考虑暴力:

很明显,每次维护与其相邻的点,将其的权值加上这个数

缺点:

复杂度过大,为$O(n^2)$,如遇到类似于菊花图的情况直接T掉

对于Subtask1:

显然是一个森林结构,并且有很多菊花图,此时我们可以考虑当修改一个点时,也顺便修改其父亲的点,并且打上一个tag,

当询问一个点时,只要把这个点的权值加上它父亲的tag即可

标算:

受到Subtask1的启示,我们可以当修改一个点时,只考虑修改其中一些相邻节点,再打上tag,查询点时,正好相反,只要把这

个点的权值加上其它未修改到这个节点的邻节点的tag

最优情况显然是度数大的由相邻节点修改时顺带维护,度数小的则加上相邻节点tag,考虑重新建图,所有点向度数大于$\sqrt n$连出边,统计修改时只要修改统计出边即可

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
// luogu-judger-enable-o2
#include<bits/stdc++.h>
using namespace std;
#define fre(x) freopen(x".in","r",stdin),freopen(x".out","w",stdout)
struct io{
    template<typename T>void read1(T &n)
    {
        n=0;char c;bool f=0;
        for(;!isdigit(c);c=getchar())if(c=='-')f=1;
        for(;isdigit(c);c=getchar())n=n*10+c-'0';
        if(f)n=-n;
    }
    template<typename T>io operator >> (T &n)
    {
        this->read1(n);return *this;
    }
}yin;
const int maxn=3*1e5+1000;
struct edge{
    int v,nxt;
}e[maxn<<1];
int in[maxn],tot,head[maxn],dd[maxn],sum[maxn],ad[maxn];
inline void __ADD(int u,int v)
{
    ++in[v];
    e[++tot].v=v;
    e[tot].nxt=head[u];
    head[u]=tot;
}
inline void add(int u,int v)
{
    __ADD(u,v);
    __ADD(v,u);
}
vector<int> a[maxn];
int main()
{
//     fre("sibyl");
    register int n,m,q;
    yin>>n>>m>>q;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        yin>>x>>y;
        add(x,y);
    }
    register int xx=sqrt(m);
    for(register int i=1;i<=n;++i)
    {
        if(in[i]>=xx)
        {
            dd[i]=1;
            for(register int j=head[i];j;j=e[j].nxt)
            {
                a[e[j].v].push_back(i);
            }
        }

    }
    while(q--)
    {
        register int t,x,y;
        yin>>t;
        if(t)
        {
            yin>>x>>y;
            sum[x]+=y;
            if(dd[x])ad[x]+=y;
            else
                for(register int i=head[x];i;i=e[i].nxt)
                {
                    sum[e[i].v]+=y;
                    // cout<<x<<endl;
                }
        }
        else
        {
            yin>>x;
            int ans=0;
            for(register int i=0;i<a[x].size();++i)
            {
                ans+=ad[a[x][i]];
                // cout<<*i<<endl;
            }
            printf("%d\n",(int)(ans+sum[x]));
        }
    }
    return 0;
}