P1330 封锁阳光大学,其实想明白了就是黑白点染色,分析的话,对于条件:”当两只河蟹封锁了相邻的两个点时,他们会发生冲突”,考虑对边的影响,则一条边,两个点选一个但只能选一个,所以不妨将模型转化为染色,不过有“不可能的情况”,这样的话其实判两端的颜色其实就行:(烦,没注意它可能不是强连通图)
using namespace std;
int n,m,pd=1;
int len=0,last[100001],col[100001],ans[3];
bool v[100001];
struct pp
{
int x,y,next;
};pp p[400001];
void ins(int x,int y)
{
int now=++len;
p[now]={x,y,last[x]};last[x]=now;
return ;
}
void dfs(int x,int lastc)
{
col[x]=lastc^1;ans[lastc^1]++;v[x]=false;
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;
if(v[y]==false&&col[y]==col[x]) pd=-1;
else
{
if(v[y]==true) dfs(y,lastc^1);
}
}
return ;
}
int main()
{
memset(last,-1,sizeof(last));memset(v,true,sizeof(v));
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
int x,y;scanf("%d%d",&x,&y);
ins(x,y);ins(y,x);
}
int kkk=0;
for(int i=1;i<=n;i++)
{
if(v[i])
{
ans[0]=ans[1]=0;
dfs(i,0),kkk+=min(ans[1],ans[0]);
}
}
if(pd==-1) printf("Impossible\n");
else printf("%d\n",kkk);
return 0;
}
今天也在努力做一题呢P2661 [NOIP2015 提高组] 信息传递这个大概是与环有关的,图论求最小环,尝试用并查集。也可以尝试用一下拓扑,来删除不在环中的节点后用dfs求解。不过我们这里讨论的是并查集的做法,并查集可以将用于判环我们是知道的,但怎么求最小路径呢,我们可以尝试采用路径压缩的方式(为什么不按序合并?每一个点的路径长及父亲都是给定了的,所以不行)如何压缩,直接从上一个节点继承即可~
假如说信息由A传递给B,那么就连一条由A指向B的边,同时更新A的父节点,A到它的父节点的路径长也就是B到它的父节点的路径长+1。如果有两个点祖先节点相同,那么就可以构成一个环,长度为两个点到祖先节点长度之和+1。
虽然但是dfs真的很简单好吗
//口是心非~~~~ no...
//这个处理方法极妙地方在于其对路径的理解之深,想到了以并查集维护路径的同时,以父亲维护长度
//具体看两句:dep[x]+=dep[last];若这一点需要更新,则直接取其父亲的长
// ans=min(ans,dep[x]+dep[y]+1);这句更新的话指的是若祖先相同的话就让自己的深度+y的+1
//不过说真的其实x与y的路径长一定有一个是0,自己yy一会就懂啦!
//+1是因为要多走一轮!
#include<bits/stdc++.h>
using namespace std;
int n,m,dep[2000001],fa[2000001],ans=1e9;
int findfa(int x)
{
if(fa[x]==x) return x;
int last=fa[x];fa[x]=findfa(fa[x]);
dep[x]+=dep[last];
return fa[x];
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++) fa[i]=i;
for(int i=1;i<=n;i++)
{
int x=i,y;scanf("%d",&y);
int fx=findfa(x),fy=findfa(y);
if(fx==fy) ans=min(ans,dep[x]+dep[y]+1);
else fa[fx]=fy,dep[x]=dep[y]+1;
}
printf("%d",ans);
return 0;
}
欸欸欸新算法:最短路树~我真棒!就是做一遍最短路之后,因为每一个的上一个点是确定了的所以直接for扫一遍看看能不能连上就行。若是要求字典序最小也不难,只要再多比较一个字典序就行。
#include<bits/stdc++.h>
#define int long long
using namespace std;
int n,m,t,len=0,last[200001],sum[200001];
int a[200001],ans=0,dis[200001],c[200001];
bool v[200001];
struct pp
{
int x,y,c,next;
};pp p[400001];
struct node
{
int x,dis;
friend bool operator < (const node &x,const node &y)
{
return x.dis > y.dis;
};
};priority_queue<node> q;
void ins(int x,int y,int c)
{
int now=++len;
p[now]={x,y,c,last[x]};last[x]=now;
return ;
}
void dj(int ST)
{
memset(dis,63,sizeof(dis));dis[ST]=0;
node e;e.x=ST;e.dis=dis[ST];q.push(e);
while(q.size()!=0)
{
node x=q.top();q.pop();
if(v[x.x]==false) continue;v[x.x]=false;
for(int i=last[x.x];i!=-1;i=p[i].next)
{
int y=p[i].y;
if(dis[y]>dis[x.x]+p[i].c)
{
dis[y]=dis[x.x]+p[i].c;
e.x=y,e.dis=dis[y];q.push(e);
}
}
}
return ;
}
void dfs(int x,int fa)
{
sum[x]=a[x];
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;if(y==fa) continue ;
dfs(y,x);sum[x]+=sum[y];
}
ans=max(ans,sum[x]*(dis[x]-t));
return ;
}
main()
{
memset(last,-1,sizeof(last));memset(v,true,sizeof(v));
scanf("%lld%lld%lld",&n,&m,&t);
for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
for(int i=1;i<=m;i++)
{
int x,y,c;scanf("%lld%lld%lld",&x,&y,&c);
ins(x,y,c);ins(y,x,c);
}
dj(1);
for(int i=1;i<=n;i++)
{
int x=i;
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y; if(dis[y]<dis[x]) continue ;
if(dis[x]+p[i].c==dis[y])
{
if(c[y]==0) c[y]=x;
if(c[y]>x) c[y]=x;
}
}
}
memset(last,-1,sizeof(last));len=0;
for(int i=2;i<=n;i++) ins(i,c[i],0),ins(c[i],i,0);
dfs(1,0);printf("%lld",ans);
return 0;
}
P1265 公路修建,用的是prim,用另一个会炸的。其实就是裸的最小生成树。不过对这玩意的板子还是不算熟。
#include<bits/stdc++.h>
using namespace std;
int n,m;
bool v[100001];
double dis[100001],x[100001],y[100001],ans=0;
double getlong(double x1,double y1,double x2,double y2)
{
return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
}
int main()
{
memset(dis,0x7f,sizeof(dis));memset(v,true,sizeof(v));
scanf("%d",&n);dis[1]=0;
for(int i=1;i<=n;i++) scanf("%lf%lf",&x[i],&y[i]);
for(int i=1;i<=n;i++)
{
int k=0;
for(int j=1;j<=n;j++)
{
if(v[j]&&dis[j]<dis[k]) k=j;
}
v[k]=false;
for(int j=1;j<=n;j++)
{
if(v[j]&&getlong(x[j],y[j],x[k],y[k])<dis[j]) dis[j]=getlong(x[j],y[j],x[k],y[k]);
}
}
for(int i=1;i<=n;i++) ans+=sqrt(dis[i]);
printf("%.2lf\n",ans);
return 0;
}
有一个欧拉路径的学习:
呃啊不写了吧刷点别的Number of Components神仙题目与大佬爆杀。看到题目大概没什么思路就开题解了,一篇就说,看到连通块直接考虑从转化入手,连通块数量=点数-边数(好神仙的想法!),然后将两个分开统计。上面那个式子可以自己手推一下,考虑剩下两个怎么统计,不妨单看一个点,若是能对结果造成贡献,那么 l<=a[ i ]&&a[ i ]<=r 所以就可以得出一个结论:i 这个点能造成的贡献是 a[ i ] * (n-a[ i ]+1) 也就是l能选择的数量-r能选择的数量。再来看边的统计,若是一条i到i+1的边被统计,那么肯定两者中的最小值>=l ,最大值<=r ,那么这将是绝杀!太强啦!绝杀!!!
#include<bits/stdc++.h>
using namespace std;
long long n,m,a[200001],ans=0;
int main()
{
scanf("%lld",&n);
for(int i=1;i<=n;i++)
{
scanf("%lld",&a[i]);ans+=a[i]*(n-a[i]+1);
}
for(int i=1;i<=n-1;i++) ans-=min(a[i],a[i+1])*(n-max(a[i],a[i+1])+1);
printf("%lld",ans);
return 0;
}
?刚刚打的板子没了?nmd
P1983 [NOIP2013 普及组] 车站分级拓扑的题目,其实没看出来,不过,呃分层之后就懂了。
#include<bits/stdc++.h>
using namespace std;
int n,m,len=0,last[4000001],dep[4000001],sta[4000001],a[4000001],ans=0;
bool v[4000001];
struct pp
{
int x,y,next;
};pp p[4000001];
void ins(int x,int y)
{
int now=++len;
p[now]={x,y,last[x]};last[x]=now;
return ;
}
void dfs(int x)
{
if(v[x]) return ;v[x]=true;
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;dfs(y);
dep[x]=max(dep[y]+1,dep[x]);
}
if(x>n) dep[x]--;ans=max(ans,dep[x]);
return ;
}
int main()
{
memset(last,-1,sizeof(last));memset(v,false,sizeof(v));
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
int s;scanf("%d",&s);
for(int j=1;j<=s;j++) scanf("%d",&a[j]),sta[a[j]]=i,ins(n+i,a[j]);
for(int j=a[1]+1;j<=a[s];j++)
{
if(sta[j]!=i) ins(j,n+i);
}
}
for(int i=1;i<=n;i++)
{
if(!v[i]) dfs(i);
}
printf("%d",ans+1);
return 0;
}
有一道带权并查集的题目:P1196 [NOI2002] 银河英雄传说,不妨维护一个区间的大小以及从一个点到根的距离,那么直接减去就行。
#include<bits/stdc++.h>
using namespace std;
int fa[300001],s[300001],siz[300001],n,m;
int findfa(int x)
{
if(fa[x]==x) return x;
int ff=fa[x];fa[x]=findfa(fa[x]);
s[x]+=s[ff];siz[x]=siz[fa[x]];//先转移,后找爹。
return fa[x];
}
int main()
{
for(int i=1;i<=300001;i++) fa[i]=i,siz[i]=1,s[i]=0;//初始化
scanf("%d",&m);
while(m--)
{
char op[10];scanf("%s",op+1);
if(op[1]=='M')
{
int x,y;scanf("%d%d",&x,&y);
int fx=findfa(x),fy=findfa(y);
fa[fx]=fy;s[fx]+=siz[fy];
siz[fx]+=siz[fy];siz[fy]=siz[fx];//考虑转移,先到根的距离,转移大小。
}
else
{
int x,y;scanf("%d%d",&x,&y);
int fx=findfa(x),fy=findfa(y);
if(fx!=fy) printf("-1\n");
else printf("%d\n",abs(s[x]-s[y])-1);
}
}
return 0;
}
慢慢做吧,后面会再见的。P1462 通往奥格瑞玛的道路
//最短路的皮子再套一个二分。
#include<bits/stdc++.h>
#define int long long
using namespace std;
int n,m,len=0,last[100001],b;//我们用边权来表示血量,那么扩展到下一个点的时候就要比较金钱了啦
int a[100001],dis[100001];
bool v[100001];
struct pp
{
int x,y,c,next;
};pp p[4000001];
struct node
{
int x,dis;
friend bool operator < (const node &x,const node &y)
{
return x.dis > y.dis;
}
};priority_queue<node > q;
void ins(int x,int y,int c)
{
int now=++len;
p[now]={x,y,c,last[x]};last[x]=now;
return ;
}
bool c_dj(int ST,int ED,int lim)
{
if(a[ST]>lim) return false;
memset(v,false,sizeof(v));memset(dis,63,sizeof(dis));dis[ST]=0;
node e;e.x=ST,e.dis=0;q.push(e);
while(q.size()!=0)
{
node x=q.top();q.pop();
if(v[x.x]) continue ; v[x.x]=true;
for(int i=last[x.x];i!=-1;i=p[i].next)
{
int y=p[i].y;
if(a[y]<=lim&&dis[y]>dis[x.x]+p[i].c)
{
dis[y]=dis[x.x]+p[i].c;
e.dis=dis[y],e.x=y;q.push(e);
}
}
}
return (dis[ED]<=b);
}
main()
{
memset(last,-1,sizeof(last));
scanf("%lld%lld%lld",&n,&m,&b);
for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
for(int i=1;i<=m;i++)
{
int x,y,c;scanf("%lld%lld%lld",&x,&y,&c);
ins(x,y,c);ins(y,x,c);
}
int l=0,r=1e9+1,ans=-1;
while(l<=r)
{
int mid=(l+r)/2;
if(c_dj(1,n,mid)) r=mid-1,ans=mid;
else l=mid+1;
}
if(ans==-1) printf("AFK");
else printf("%lld",ans);
return 0;
}
P2921 [USACO08DEC]Trick or Treat on the Farm G之前比赛的题目呢,不过那时候拿的是寄搜过的。那么这次尝试一下缩点?好吧最后用记忆化优化了一下不过有些细节还是难搞。
#include<bits/stdc++.h>
using namespace std;
int n,m,len=0,last[200001],dfsx[200001],low[200001],dfs_x=0;
int q[200001],tot=0,id[200001],point=0,siz[200001],ans[200001],f[200001];
bool v[200001];
struct pp
{
int x,y,next;
};pp p[200001];
void ins(int x,int y)
{
int now=++len;
p[now]={x,y,last[x]};last[x]=now;
return ;
}
void SH(int x)
{
low[x]=dfsx[x]=++dfs_x;
v[x]=true;q[++tot]=x;
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;
if(!dfsx[y]) SH(y),low[x]=min(low[x],low[y]);
else if(v[y]) low[x]=min(low[x],dfsx[y]);
}
if(dfsx[x]==low[x])
{
point++;
while(dfsx[x]<=dfsx[q[tot]])
{
id[q[tot]]=point;siz[point]++;
v[q[tot--]]=false;dfs_x--;
}
}
return ;
}
void dfs(int x)
{
if(f[x]) return ;//printf("%d ",x);
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;dfs(y);f[x]=f[y]+1;
}
return ;
}
int main()
{
memset(last,-1,sizeof(last));memset(v,false,sizeof(v));
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
int x;scanf("%d",&x);
ins(i,x);if(i==x) f[i]=ans[i]=1;
}
for(int i=1;i<=n;i++)
{
if(!dfsx[i]) SH(i);//printf("%d ",i);
}
for(int i=1;i<=n;i++)
{
if(siz[id[i]]!=1) f[i]=ans[i]=siz[id[i]];
}
for(int i=1;i<=n;i++)
{
if(!ans[i]) dfs(i),ans[i]=f[i];
}
for(int i=1;i<=n;i++) printf("%d\n",ans[i]);
return 0;
}
二分图水题,不过要黑白染色:P6268 [SHOI2002]舞会,太久没用二分图,手法生疏了调了一会。
#include<bits/stdc++.h>
using namespace std;
int n,m,len=0,last[100001],ans=0,color[100001],op[1000001];
bool v[100001];
struct pp
{
int x,y,next;
};pp p[100001];
void ins(int x,int y)
{
int now=++len;
p[now]={x,y,last[x]};last[x]=now;
return ;
}
void getcolor(int x,int c)
{
color[x]=c;//printf("%d ",c);
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;if(color[y]!=-1) continue ;
getcolor(y,c^1);
}
return ;
}
bool dfs(int x)
{
for(int i=last[x];i!=-1;i=p[i].next)
{
int y=p[i].y;if(v[y]) continue ;v[y]=true;
if((!op[y])||dfs(op[y]))
{
op[y]=x;
return true ;
}
}
return false;
}
int main()
{
memset(last,-1,sizeof(last));memset(color,-1,sizeof(color));memset(op,0,sizeof(op));
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
int x,y;scanf("%d%d",&x,&y);
x++,y++;ins(x,y),ins(y,x);
}
for(int i=1;i<=n;i++)
{
if(color[i]==-1) getcolor(i,0);
}
for(int i=1;i<=n;i++)
{
if(color[i]!=1) continue ;
memset(v,false,sizeof(v));//printf("%d ",i);
if(dfs(i)) ans--;//printf("%d ",i);
}
ans+=n;printf("%d",ans);
return 0;
}