地理课(geography)

题目描述

地理课上，老师给出了一个巨大的地图，由于世界日新月异，会有一些道路在某一时刻被删除，也会有一些道路在某一时刻被修建。这里的道路均为双向的。

老师认为，有一些城市被分在了一个连通块中可以相互到达，而有一些城市不能够相互到达。而他想知道，每个时刻所有连通块大小的乘积是多少？

wzy看到这个地图的时候就蒙了，还好那只上天的喵及时帮助了他。现在他把这个毒瘤的地图拿过来给你，想试试看你能不能求出来。由于答案可能很大，输出乘积mod109+7mod109+7即可。

输入

第一行两个数n,mn,m，表示有nn个点，mm个时刻。接下来mm行每行三个数，要么是1uv1uv，要么是2uv2uv，分别表示添加一条无向边和删除一条无向边。

输出

共mm，每行一个数表示连通块大小乘积mod1,000,000,007mod1,000,000,007

上面是每个时刻操作后的图。乘积分别为：

2*1*1*1=2,3*1*1=3,3*1*1=3,3*2=6,5,3*2=6

数据范围及约定

subtask1：30pts，n≤1,000,m≤2,000n≤1,000,m≤2,000

subtask2：20pts，满足没有删除操作。

subtask3：50pts，n,m≤100,000n,m≤100,000保证没有重边自环，不会删除不存在的边。

来源

noip2018南外-模拟

solution

好题好题

考虑离线，按时间建一棵线段树

把每一个边的编号加进线段树对应的区间中

我们现在希望得到每一个叶子节点的值

那么每次暴力求是不行了

我们记下每一层的边连接的并查集的根与大小

返回时把边拆开即可

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<map>
#include<vector>
#define mod 1000000007
#define maxn 100005
#define ll long long
using namespace std;
int n,m,f[maxn],sz[maxn],t1,t2,op;
ll now,ans[maxn],ny[maxn];
struct node{int u,v,st,ed;
}e[maxn];
map<ll,int>ls;
struct no{int l,r;vector<int>x;
}tree[maxn*8];
struct Node{int id,x,y,size;
};
vector<Node>p;
void build(int k,int L,int R){tree[k].l=L,tree[k].r=R;if(L==R)return;int mid=L+R>>1;build(k*2,L,mid);build(k*2+1,mid+1,R);
}
void add(int k,int L,int R,int v){//cout<<k<<' '<<L<<' '<<R<<''<<endl;if(tree[k].l>=L&&tree[k].r<=R){tree[k].x.push_back(v);return;}int mid=tree[k].l+tree[k].r>>1;if(L<=mid)add(k*2,L,R,v);if(R>mid)add(k*2+1,L,R,v);
}
ll work(ll a,int num){ll an=1;while(num){if(num&1)an=an*a;a=a*a;a%=mod;an%=mod;num>>=1;}return an;
}
int getf(int k){if(f[k]==k)return k;return getf(f[k]);
}
void hb(int k){//cout<<"k: "<<k<<' '<<tree[k].l<<' '<<tree[k].r<<endl;for(int i=0;i<tree[k].x.size();i++){int t=tree[k].x[i];//cout<<t<<endl;int f1=getf(e[t].u),f2=getf(e[t].v);if(f1!=f2){now=now*ny[sz[f1]]%mod*ny[sz[f2]]%mod;now=now*(sz[f1]+sz[f2])%mod;if(sz[f1]<sz[f2]){p.push_back((Node){k,f1,f2,sz[f1]});f[f1]=f2;sz[f2]+=sz[f1];sz[f1]=0;}else {p.push_back((Node){k,f2,f1,sz[f2]});f[f2]=f1;sz[f1]+=sz[f2];sz[f2]=0;}}}
}
void turn(int k){int Size=p.size()-1;if(Size<0)return;//cout<<"Size "<<Size<<' '<<tree[k].l<<' '<<tree[k].r<<endl;for(int i=Size;p[i].id==k&&i>=0;i--){//cout<<p[i].x<<' '<<p[i].y<<' '<<sz[p[i].x]<<' '<<sz[p[i].y]<<endl;now=now*ny[sz[p[i].y]]%mod;f[p[i].x]=p[i].x;sz[p[i].x]=p[i].size;sz[p[i].y]-=sz[p[i].x];now=now*sz[p[i].x]%mod*sz[p[i].y]%mod;p.pop_back();}
}
void dfs(int k){hb(k);if(tree[k].l==tree[k].r){ans[tree[k].l]=now;turn(k);return;}dfs(k*2);dfs(k*2+1);turn(k);
}
int main()
{cin>>n>>m;for(int i=1;i<=m;i++){scanf("%d%d%d",&op,&e[i].u,&e[i].v);if(e[i].u>e[i].v)swap(e[i].u,e[i].v);ll tmp=(ll)e[i].u*100000+e[i].v;if(op==1){ls[tmp]=i;e[i].st=i;}else {int la=ls[tmp];e[la].ed=i-1;//ls[tmp]=0;}}build(1,1,m);for(int i=1;i<=m;i++){if(e[i].ed==0)e[i].ed=m;if(e[i].st&&e[i].ed)add(1,e[i].st,e[i].ed,i);}for(int i=1;i<=n;i++)ny[i]=work(i,mod-2);for(int i=1;i<=n;i++)f[i]=i,sz[i]=1;now=1;dfs(1);for(int i=1;i<=m;i++)printf("%lld\n",ans[i]);return 0;
}