【2019.9.16】Za

其他 2019-09-19 22:31:10 阅读次数: 0

天天爱跑步

树上差分

咕==

#include<bits/stdc++.h>
using namespace std;
const int N=3e5+10,M=50000+10,inf=0x3f3f3f3f;
int n,m,w[N],ans[N],c1[N<<1],c2[N<<1];
vector<int>a1[N],b1[N],a2[N],b2[N];
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int head[N],tot=0;
struct edge{int v,nxt;}e[N<<1];
void add(int u,int v){e[++tot]=(edge){v,head[u]},head[u]=tot;}
int dep[N],f[N][25];
void bfs(int u,int fa){
    dep[u]=dep[fa]+1,f[u][0]=fa;
    for(int i=head[u],v;i;i=e[i].nxt)
    if((v=e[i].v)!=fa) bfs(v,u);
}
void prepare(){
    bfs(1,0);
    for(int i=1;i<=20;++i)
        for(int j=1;j<=n;++j)
            f[j][i]=f[f[j][i-1]][i-1];
}
int LCA(int x,int y){
    if(dep[x]<dep[y]) swap(x,y);
    for(int i=20;i>=0;--i)
        if(dep[f[x][i]]>=dep[y]) x=f[x][i];
    if(x==y) return x;
    for(int i=20;i>=0;--i)
    if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
    return f[x][0];
}

bool vis[N];
void dfs(int x){
    int w1=c1[dep[x]+w[x]],w2=c2[w[x]-dep[x]+n];
    vis[x]=1;
    for(int i=head[x],v;i;i=e[i].nxt)
    if(!vis[v=e[i].v]) dfs(v);
    for(int i=0;i<a1[x].size();++i) ++c1[a1[x][i]];
    for(int i=0;i<b1[x].size();++i) --c1[b1[x][i]];
    for(int i=0;i<a2[x].size();++i) ++c2[a2[x][i]+n];
    for(int i=0;i<b2[x].size();++i) --c2[b2[x][i]+n];
    ans[x]+=c1[dep[x]+w[x]]-w1+c2[w[x]-dep[x]+n]-w2;
}


int main(){
    freopen("in.txt","r",stdin);
    rd(n),rd(m);
    for(int i=1,u,v;i<n;++i) rd(u),rd(v),add(u,v),add(v,u);
    for(int i=1;i<=n;++i) rd(w[i]);
    prepare();
    for(int i=1,x,y,z;i<=m;++i){
        rd(x),rd(y),z=LCA(x,y);
        a1[x].push_back(dep[x]),
        a2[y].push_back(dep[x]-(dep[z]<<1)),
        b1[f[z][0]].push_back(dep[x]),
        b2[z].push_back(dep[x]-(dep[z]<<1));
    }
    dfs(1);
    for(int i=1;i<=n;++i) printf("%d ",ans[i]);
    return 0;
}

求桥

void tarjan(int u,int ine){
    dfn[u]=low[u]=++idx;
    for(int i=head[u],v;i;i=e[i].nxt){
        if(!dfn[v=e[i].v]) {
            tarjan(v,i),low[u]=min(low[u],low[v]);
            if(low[v]>dfn[u]) bridge[i]=bridge[i^1]=1;
        }
        else if(i!=(ine^1)) low[u]=min(low[u],dfn[v]);
    }
}

求割点

void tarjan(int u){
    dfn[u]=low[u]=++idx;int kid=0;
    for(int i=head[u],v;i;i=e[i].nxt){
        if(!dfn[v=e[i].v]){
            tarjan(v),low[u]=min(low[u],low[v]);
            if(low[v]>=dfn[u])
                if(++kid>1||u!=rt) cut[u]=1;
        }
        else low[u]=min(low[u],dfn[v]);
    }
}

边双连通e-DCC

void tarjan(int u,int ine){
    dfn[u]=low[u]=++idx;
    for(int i=head[u],v;i;i=e[i].nxt){
        if(!dfn[v=e[i].v]) {
            tarjan(v,i),low[u]=min(low[u],low[v]);
            if(low[v]>dfn[u]) bridge[i]=bridge[i^1]=1;
        }
        else if(i!=(ine^1)) low[u]=min(low[u],dfn[v]);
    }
}
int bl[N],dcc;
void dfs(int x){
    bl[x]=dcc;
    for(int i=head[u],v;i;i=e[i].nxt)
        if(!bl[v=e[i].v]&&!bridge[i]) dfs(v);
}
int main(){
    for(int i=1;i<=n;++i)
        if(!bl[i]) ++dcc,dfs(i);
    return 0;
}

点双连通v-DCC

void tarjan(int u){
    dfn[u]=low[u]=++idx,s.push(u);
    if(u==rt&&!head[u]) dcc[++cnt].push_back(u);
    int kid=0;
    for(int i=head[u],v,;i;i=e[i].nxt){
        if(!dfn[v=e[i].v]){
            tarjan(v),low[u]=min(low[u],low[v]);
            if(low[v]>=low[u]){
                if(++kid>1||u!=rt) cut[u]=1;
                ++cnt;int z;
                do{
                    z=s.top(),s.pop();
                    dcc[cnt].push_back(z);
                }while(z!=v)
                    dcc[cnt].push_back(u);
            }
            else low[u]=min(low[u],dfn[v]);
        }
    }
}

[JSOI2010] 连通数

bzoj2208 luogu4306

缩点完后用bitset来统计==

按拓扑序来

#include<bits/stdc++.h>
using namespace std;
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)<(y)?(x):(y))
const int N=2000+10,M=4000000+10,inf=0x3f3f3f3f;
int n;
char ss[N];
bitset<N>f[N];
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int head[N],hd[N],tot=0,tt=0;
struct edge{int v,nxt;}e[M],E[M];
void Add(int u,int v){E[++tt]=(edge){v,hd[u]},hd[u]=tt;}
void add(int u,int v){e[++tot]=(edge){v,head[u]},head[u]=tot;}

int idx=0,Bcnt=0,dfn[N],low[N],bl[N],sz[N];
stack<int>s;bool inst[N];
void tarjan(int u){
    dfn[u]=low[u]=++idx,s.push(u),inst[u]=1;
    for(int i=head[u],v;i;i=e[i].nxt)
        if(!dfn[v=e[i].v]) tarjan(v),low[u]=Min(low[u],low[v]);
        else if(inst[v]&&dfn[v]<low[u]) low[u]=dfn[v];
    if(dfn[u]==low[u]){
        ++Bcnt;int v;
        do{
            v=s.top(),s.pop(),inst[v]=0,f[Bcnt].set(v);
            bl[v]=Bcnt,++sz[Bcnt];
        }while(u!=v);
    }
}

int Q[N],in[N];
queue<int>q;
void topsort(){
    for(int i=1;i<=Bcnt;++i) if(!in[i]) q.push(i);
    while(!q.empty()){
        int u=q.front();
        q.pop(),Q[++Q[0]]=u;
        for(int i=hd[u],v;i;i=E[i].nxt)
            if(!--in[v=E[i].v]) q.push(v);
    }
    for(int i=Q[0];i;--i)
        for(int j=hd[Q[i]];j;j=E[j].nxt)
            f[Q[i]]|=f[E[j].v];
}

int main(){
    freopen("in.txt","r",stdin);
    rd(n);
    for(int i=1;i<=n;++i){
        scanf("%s",ss+1);
        for(int j=1;j<=n;++j)
            if(ss[j]=='1') add(i,j);
    }
    for(int i=1;i<=n;++i)
    if(!dfn[i]) tarjan(i);
    for(int u=1,v;u<=n;++u)
        for(int i=head[u];i;i=e[i].nxt)
            if(bl[u]!=bl[v=e[i].v]) Add(bl[u],bl[v]),++in[bl[v]];
    int ans=0;
    topsort();
    for(int i=1;i<=Bcnt;++i) ans+=f[i].count()*sz[i];
    printf("%d",ans);
    return 0;
}

缩完点后topsort时DP

#include<bits/stdc++.h>
using namespace std;
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)<(y)?(x):(y))
const int N=100000+10,M=100000 +10,inf=0x3f3f3f3f;
int n,m,in[N];
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int head[N],hd[N],tot=0,tt=0;
struct edge{int v,w,nxt;}e[M],E[M<<1];
void Add(int u,int v,int w){E[++tt]=(edge){v,w,hd[u]},hd[u]=tt;}
void add(int u,int v,int w){e[++tot]=(edge){v,w,head[u]},head[u]=tot;}
int idx=0,Bcnt=0,dfn[N],low[N],bl[N],sz[N];
stack<int>s;bool inst[N];
void tarjan(int u){
    dfn[u]=low[u]=++idx,s.push(u),inst[u]=1;
    for(int i=head[u],v;i;i=e[i].nxt)
        if(!dfn[v=e[i].v]) tarjan(v),low[u]=Min(low[u],low[v]);
        else if(inst[v]&&dfn[v]<low[u]) low[u]=dfn[v];
    if(low[u]==dfn[u]){
        ++Bcnt;int v;
        do{
            v=s.top(),s.pop(),inst[v]=0;
            bl[v]=Bcnt,++sz[Bcnt];
        }while(u!=v);
    }
}

queue<int>q;
int f[N];
void topsort(){
    q.push(0);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=hd[u],v;i;i=E[i].nxt){
            f[v=E[i].v]=Max(f[v],f[u]+E[i].w);
            if(!--in[v]) q.push(v);
        }
    }
    long long ans=0;
    for(int i=1;i<=Bcnt;++i) ans+=(long long)f[i]*sz[i];
    printf("%lld",ans);
}

int main(){
    freopen("in.txt","r",stdin);
    rd(n),rd(m);
    for(int i=1,opt,x,y;i<=m;++i){
        rd(opt),rd(x),rd(y);
        if(opt==1) add(x,y,0),add(y,x,0);
        else if(opt==2) add(x,y,1);
        else if(opt==3) add(y,x,0);
        else if(opt==4) add(y,x,1);
        else add(x,y,0);
    }
    for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);
    for(int u=1;u<=n;++u)
        for(int i=head[u],v;i;i=e[i].nxt){
            if(bl[u]==bl[v=e[i].v]&&e[i].w) return puts("-1"),0;
            else if(bl[u]!=bl[v]) Add(bl[u],bl[v],e[i].w),++in[bl[v]];
        }
    for(int i=1;i<=Bcnt;++i) Add(0,i,1),++in[i];
    topsort();
    return 0;
}

[IOI2008]Island

我太难了 打的时候脑子里一团浆糊

acwing上面要T 洛谷上面可以过

#include<bits/stdc++.h>
using namespace std;
#define Max(x,y) ((x)>(y)?(x):(y))
#define Min(x,y) ((x)<(y)?(x):(y))
#define ll long long
const int N=1e6+10,M=100000 +10,inf=0x3f3f3f3f;
int n,in[N];
template <class t>void rd(t &x){
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

int head[N],tot=0;
struct edge{int v,w,nxt;}e[N<<1];
void add(int u,int v,int w){
    e[++tot]=(edge){v,w,head[u]},head[u]=tot,++in[v];
}

int bl[N],Bcnt=0;
queue<int>q;
void bfs(int u){
    q.push(u);
    while(!q.empty()){
        u=q.front(),q.pop();
        for(int i=head[u],v;i;i=e[i].nxt)
            if(!bl[v=e[i].v]) bl[v]=Bcnt,q.push(v);
    }
}

ll dis[N],f[N],a[N],b[N];
void topsort(){
    for(int i=1;i<=n;++i) if(in[i]==1) q.push(i);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u],v;i;i=e[i].nxt)
        if(in[v=e[i].v]>1){
            dis[bl[u]]=Max(dis[bl[u]],f[u]+f[v]+e[i].w),
            f[v]=Max(f[v],f[u]+e[i].w);
            if((--in[v])==1) q.push(v);
        }
    }
}

int Q[N];
bool vis[N];
void dp(int x,int t)
{
    int m=0,i,l=0,r,y=x;
    do{
        a[++m]=f[y];in[y]=1;
        for(i=head[y];i;i=e[i].nxt)
            if(in[e[i].v]>1)
                {y=e[i].v;b[m+1]=b[m]+e[i].w; break;}
    }while(i);
    if(m==2)
    {
        for(i=head[y];i;i=e[i].nxt)
            if(e[i].v==x) {l=max(l,e[i].w);}
        dis[t]=max(dis[t],f[x]+f[y]+l);
        return;
    }
    for(i=head[y];i;i=e[i].nxt)
        if(e[i].v==x) {b[m+1]=b[m]+e[i].w; break;}
    for(i=1;i<m;i++) {a[m+i]=a[i]; b[m+i]=b[m+1]+b[i];}
    Q[l=r=1]=1;
    for(i=2;i<2*m;i++)
    {
        while(l<=r&&i-Q[l]>=m) l++;
        dis[t]=max(dis[t],a[i]+a[Q[l]]+b[i]-b[Q[l]]);
        while(l<=r&&a[Q[r]]+b[i]-b[Q[r]]<=a[i]) r--;
        Q[++r]=i;
    }
}


int main(){
    freopen("in.txt","r",stdin);
    rd(n);ll ans=0;
    for(int i=1,v,w;i<=n;++i) rd(v),rd(w),add(i,v,w),add(v,i,w);
    for(int i=1;i<=n;++i)
    if(!bl[i]) ++Bcnt,bfs(i);
    topsort();
    for(int i=1;i<=n;++i)
    if(in[i]>1&&!vis[bl[i]]) vis[bl[i]]=1,dp(i,bl[i]),ans+=dis[bl[i]];
    printf("%lld",ans);
    return 0;
}

猜你喜欢

转载自www.cnblogs.com/lxyyyy/p/11552927.html
za
今日推荐
面壁智能发布 Eurux-8x22B 开源大模型 —— 堪称「理科状元」
开源日报 | 谷歌扶持鸿蒙上位；开源Rabbit R1；Docker加持的安卓手机；微软的焦虑和野心；海尔电器把开放平台关了
中国码农的“35岁魔咒”
蘭雅 CorelDRAW 插件 2024.5.1 国际劳动节版，免费下载
Arc Browser for Windows 1.0 正式 GA
90后程序员开发视频搬运软件、不到一年获利超 700 万，结局很刑！
周排行
OOP第二次作业
java web 乱码问题
android 禁止scrollview 因控件变化自动滚动到底的方法
mysql服务解压版的安装(5.7)
centos7 nginx+tomcat配置https 安装免费SSL Let’s Encrypt
使用Mosquitto遗嘱机制实现感知客户端上下线功能的方法
面向对象之------多态与多态性
开发Teams Tabs应用程序
C# 希尔排序
第2章 Jupyter Notebooks
每日归档
更多
2024-05-06(40)
2024-05-05(0)
2024-05-04(7)
2024-05-03(19)
2024-05-02(0)
2024-05-01(4)
2024-04-30(1)
2024-04-29(40)
2024-04-28(0)
2024-04-27(56)

代码天地 © 2018 codetd.com

境外服务器   最新电视剧2022