「10.17-10.18」liu_runda’s模拟
暂咕
$day1$
A. 位运算
分类讨论，贡献分离。
B. 集合论
维护类似时间戳的东西
C. 连连看
我们可以看出两个方块能作出贡献，实际上是一个极大联通块(白块)所连节点，
然后我们就可以乱搞了
因为块会重复计算，容斥即可
有点懒，所以懒得细写了....
#include<bits/stdc++.h> #define int long long #define MAXN 1101 using namespace std; int read(){ char c=getchar();int ff=1;int x=0; while(c<'0'||c>'9'){if(c=='-')ff=-1;c=getchar();} while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+(c^48);c=getchar();} return x*ff; } int dex[5]={-1,1,0,0}; int dey[5]={0,0,-1,1}; int a[MAXN][MAXN]; //int kt[4][MAXN*MAXN*2]; const int mod=23333333; /*struct node{ int head[mod];int tot=0;int nxt[mod];int to[mod];int val[mod]; int find(int x){ int me=x%mod; for(int i=head[me];i;i=nxt[i]){ if(to[i]==x){return val[i];} } return 0; } void init(){ for(int i=1;i<=tot;++i)nxt[i]=0; tot=0; } void insert(int x){ int me=x%mod; for(int i=head[me];i;i=nxt[i]){ if(to[i]==x){val[x]++;return ;} } tot++; to[tot]=x;head[me]=tot;val[tot]++; } void erase(int x){ int me=x%mod; for(int i=head[me];i;i=nxt[i]){ to[i]=0;val[i]=0; } head[me]=0; return ; } }T1,T2,T3;*/ int f[MAXN][MAXN][5]; int vis[MAXN][MAXN]; vector<pair<int,int> >v[MAXN*MAXN*2];int ans=0; void work(int x){ if(!v[x].size())return ; //printf("x=%lld\n",x); map<int,int>T4,T1,T2,T3; //kt[1][0]=0;kt[2][0]=0;kt[3][0]=0; for(int i=0;i<v[x].size();++i){ int tox=v[x][i].first;int toy=v[x][i].second; int k=1;for(k=1;k<=5;++k)if(f[tox][toy][k]==0)break;k--; for(int j=1;j<=k;++j){ /* ans+=T1.find(f[tox][toy][j]); T1.insert(f[tox][toy][j]); */ ans+=T1[f[tox][toy][j]]; T1[f[tox][toy][j]]++; //kt[1][++kt[1][0]]=f[tox][toy][j]; } for(int j=1;j<=k;++j){ unsigned long long base=f[tox][toy][j]; for(int w=1;w<j;++w){ base=base*10000007ll+f[tox][toy][w]; /* ans-=T2.find(base); T2.insert(base); kt[2][++kt[2][0]]=base; */ ans-=T2[base]; T2[base]++; } } for(int j=1;j<=k;++j){ unsigned long long base=f[tox][toy][j]; for(int w=1;w<j;++w){ base=base*10000007ll+f[tox][toy][w]; for(int z=1;z<w;++z){ base=base*10000007ll+f[tox][toy][z]; /* ans+=T3.find(base); T3.insert(base); */ ans+=T3[base]; T3[base]++; // kt[3][++kt[3][0]]=base; } } } if(k==4){ unsigned long long base=f[tox][toy][4]*10000007ll; base=f[tox][toy][3]+base*10000007ll; base=f[tox][toy][2]+base*10000007ll; base=base*100000007ll+f[tox][toy][1]; ans-=T4[base]; T4[base]++; } } /* for(int i=1;i<=kt[1][0];++i){T1.erase(kt[1][i]);} for(int i=1;i<=kt[2][0];++i){T2.erase(kt[2][i]);} for(int i=1;i<=kt[3][0];++i){T3.erase(kt[3][i]);} T1.init();T2.init();T3.init(); */ } queue<pair<int,int> >q;int n,m; void puss(int tox,int toy,int sum){ int j=0; for(j=1;j<=4;++j){if(f[tox][toy][j]==0)break;} f[tox][toy][j]=sum; } int cnt=0; pair<int,int>st[MAXN*MAXN]; void BFS(int x,int y){ int sss=0; q.push(make_pair(x,y));vis[x][y]=1; while(!q.empty()){ int mx=q.front().first; int my=q.front().second; q.pop(); for(int i=0;i<=3;++i){ int tox=mx+dex[i];int toy=my+dey[i]; if(tox<1||tox>n)continue; if(toy<1||toy>m)continue; if(a[tox][toy]==0&&vis[tox][toy]==0){ vis[tox][toy]=1; q.push(make_pair(tox,toy)); } else if(a[tox][toy]!=0&&vis[tox][toy]==0){ vis[tox][toy]=1; st[++sss]=make_pair(tox,toy); puss(tox,toy,cnt); } } } for(int i=1;i<=sss;++i){ int tox=st[i].first;int toy=st[i].second; vis[tox][toy]=0; } sss=0; } int look(int x1,int y1,int x2,int y2){ for(int i=1;i<=4;++i){ if(f[x1][y1][i]==0)continue; for(int j=1;j<=4;++j){ if(f[x2][y2][j]==0)continue; if(f[x1][y1][i]==f[x2][y2][j])return 0; } } return 1; } void QAQ(){ int sum=0; for(int i=1;i<=n;++i){ for(int j=1;j<=m;++j){ if(a[i][j]==0)continue; for(int k=0;k<=3;++k){ int tox=i+dex[k];int toy=j+dey[k]; if(tox<=0||tox>n)continue; if(toy<=0||toy>m)continue; if(a[i][j]!=a[tox][toy])continue; sum+=look(i,j,tox,toy); } } } ans+=sum/2; } int mx=0;int k; signed main(){ //freopen("link6.in","r",stdin); n=read();m=read();k=read(); for(int i=1;i<=n;++i) for(int j=1;j<=m;++j){ a[i][j]=read(); if(a[i][j]==0)continue; v[a[i][j]].push_back(make_pair(i,j)); mx=max(a[i][j],mx); } for(int i=1;i<=n;++i){ for(int j=1;j<=m;++j){ if(!vis[i][j]&&a[i][j]==0){ cnt++; BFS(i,j); } } } for(int i=1;i<=mx;++i){work(i);} QAQ(); printf("%lld\n",ans); } /* 3 3 4 1 0 2 2 3 0 0 0 2 */
$day2$
A. 串串香
$hash$即可
B. 糊涂图
还没改出来QAQ....
C. 木叶下
树上直径，再用上$lca$即可
我们维护若干数组
$dp_{i}$表示从$x$的父亲不经过$x$的最大路径
$g_{i,4}$表示$x$的几个长链的儿子
$gofa_{i}$表示向$fa$的最长路径
$dis$是对于$dp$的倍增数组
然后根据数组定义打就没了
#include<bits/stdc++.h> #define int long long #define MAXN 210005 using namespace std; int read(){ char c=getchar();int ff=1;int x=0; while(c<'0'||c>'9'){if(c=='-')ff=-1;c=getchar();} while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+(c^48);c=getchar();} return x*ff; } int anss,diss[MAXN]; int dp[MAXN];//从x的父亲不经过x的最大路径 int fa[MAXN][25]; int dis[MAXN][25]; int g[MAXN][4];//从x到x的子树的最大路径 int son[MAXN][4]; int deep[MAXN]; int gofa[MAXN];//向fa走的最长路径 struct node{int to,n;}e[MAXN*2];int tot,head[MAXN]; void add(int u,int v){e[++tot].to=v;e[tot].n=head[u];head[u]=tot;} void pianfen(int x,int fa){ for(int i=head[x];i;i=e[i].n){ int to=e[i].to; if(to==fa)continue; pianfen(to,x); anss=max(anss,diss[x]+diss[to]+1); diss[x]=max(diss[to]+1,diss[x]); } } int n,m; void DP(int x,int fa){ priority_queue<pair<int,int> >q; g[x][1]=1; for(int i=head[x];i;i=e[i].n){ int to=e[i].to; if(to==fa)continue; DP(to,x); q.push(make_pair(g[to][1]+1,to)); //printf("psg=%lld x=%lld to=%lld\n",g[to][1]+1,x,to); } int th=1; while(!q.empty()){ g[x][th]=q.top().first; son[x][th]=q.top().second; q.pop(); //printf("g[%lld][%lld]=%lld %lld\n",x,th,g[x][th],son[x][th]); th++; if(th==4)break; } } void DFS2(int x,int fa){ for(int i=head[x];i;i=e[i].n){ int to=e[i].to; if(to==fa)continue; if(son[x][1]!=to)dp[to]=max(1ll,g[x][1]); else dp[to]=max(1ll,g[x][2]); //printf("dp[%lld]=%lld\n",to,dp[to]); DFS2(to,x); } } void jia_swaproot(int x,int fa){ for(int i=head[x];i;i=e[i].n){ int to=e[i].to; if(to==fa)continue; gofa[to]=max(gofa[x]+1,dp[to]+1); jia_swaproot(to,x); } } void DFS(int x,int faa){ //printf("xxx=%lld\n",x); for(int i=head[x];i;i=e[i].n){ int to=e[i].to; if(to==faa)continue; fa[to][0]=x; deep[to]=deep[x]+1; dis[to][0]=dp[to]; //printf("dis[%lld]=%lld\n",to,dis[to][0]); for(int j=1;j<=19;++j){ if(fa[to][j-1]==0)break; fa[to][j]=fa[fa[to][j-1]][j-1]; dis[to][j]=max(dis[to][j-1],dis[fa[to][j-1]][j-1]); } DFS(to,x); } } int Lca(int x,int y){ int maxn=0; if(deep[x]<deep[y])swap(x,y); int yuanx=x;int yuany=y; //printf("x=%lld y=%lld\n",x,y); for(int i=19;i>=0;--i){ if(deep[fa[x][i]]>=deep[y]){ maxn=max(dis[x][i]-1,maxn); //printf("maxn=%lld %lld\n",maxn,dis[x][i]); x=fa[x][i]; } } if(x==y){ maxn=max(maxn,g[yuanx][1]-1); maxn=max(gofa[x]-1,maxn); //printf("go=%lld %lld\n",gofa[x]-1,x); return maxn; } for(int i=19;i>=0;--i){ if(fa[x][i]!=fa[y][i]){ maxn=max(dis[x][i]-1,maxn); maxn=max(dis[y][i]-1,maxn); x=fa[x][i],y=fa[y][i]; } } int lca=fa[x][0]; for(int j=1;j<=3;++j){ if(son[lca][j]==x||son[lca][j]==y)continue; maxn=max(maxn,g[lca][j]-1); break; } maxn=max(maxn,gofa[lca]-1); maxn=max(maxn,g[yuanx][1]-1);maxn=max(maxn,g[yuany][1]-1); return maxn; } signed main(){ //freopen("myx.in","r",stdin); //freopen("wwb.out","w",stdout); n=read(); for(int i=1;i<=n-1;++i){ int u=read();int v=read(); add(u,v);add(v,u); } deep[1]=1; pianfen(1,0); DP(1,0); DFS2(1,0); jia_swaproot(1,0); DFS(1,0); anss++; if(anss%2==1)anss/=2,anss++; else anss/=2; m=read(); for(int i=1;i<=m;++i){ int u=read();int v=read(); if(u==v){ printf("%lld\n",anss); continue; } int ans=Lca(u,v); printf("%lld\n",ans); } }
