牛客练习赛25

前言

---

我好菜啊。本来是奔着T恤去的

A

---

要求
ans=∑ni=1∑j|i1 a n s = ∑ i = 1 n ∑ j | i 1
注意到我们并不需要什么奇迹银壳，只需要交换枚举主体就可以√n做了

#include 
#include 
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

typedef long long LL;

void solve(LL n) {
    LL ans=0;
    for (LL i=1,j;i<=n;i=j+1) {
        j=n/(n/i);
        ans=(ans+(n/i)*(j-i+1));
    }
    printf("%lld\n", ans);
}

int main(void) {
    int T; scanf("%d",&T);
    while (T--) {
        LL n; scanf("%lld",&n);
        solve(n);
    }
    return 0;
}

B

---

线段树裸题

#include 
#include 
#include 
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

const int N=200005;

int s[N<<2],ls[N<<2],rs[N<<2],a[N];

int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}

void modify(int now,int tl,int tr,int x,int v) {
    if (tl==tr) {
        ls[now]=rs[now]=s[now]=1;
        return (void) (a[tl]=v);
    }
    int mid=(tl+tr)>>1;
    if (x<=mid) modify(now<<1,tl,mid,x,v);
    else modify(now<<1|1,mid+1,tr,x,v);
    ls[now]=ls[now<<1];
    if (ls[now]==(mid-tl+1)&&a[mid]1]) ls[now]+=ls[now<<1|1];
    rs[now]=rs[now<<1|1];
    if (rs[now]==(tr-mid)&&a[mid]1]) rs[now]+=rs[now<<1];
    s[now]=std:: max(s[now<<1],s[now<<1|1]);
    if (a[mid]1]) s[now]=std:: max(s[now],rs[now<<1]+ls[now<<1|1]);
}

int main(void) {
    int n=read(),m=read();
    rep(i,1,n) {
        int x=read();
        modify(1,1,n,i,x);
    }
    printf("%d\n", s[1]);
    for (;m--;) {
        int x=read(),y=read();
        modify(1,1,n,x,y);
        printf("%d\n", s[1]);
    }
    return 0;
}

C

---

注意到每一个答案前的系数是一样的，因此我们只需要O(t)求一波系数然后乘上对应位置的a即可

#include 
#include 
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)

typedef long long LL;
const int MOD=1000000007;
const int N=200005;

LL a[N],b[N],sum;
int n,m;

int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}

int main(void) {
    n=read(),m=read();
    rep(i,1,n) {
        a[i]=read();
        sum=(sum+a[i])%MOD;
    }
    LL tmp=1;
    rep(i,1,100000) {
        b[i]=(sum*tmp%MOD-b[i-1]+MOD)%MOD;
        tmp=tmp*(n-1)%MOD;
    }
    for (;m--;) {
        int x=read(),t=read();
        if (!t) printf("%lld\n", a[x]);
        else printf("%lld\n", (b[t]+((t&1)?(-1):(1))*a[x]+MOD)%MOD);
    }
    return 0;
}

D

---

我并没有做出来，不过可以猜想一波是奇怪筛+指针乱扫的做法

E

---

首先图不联通肯定impossible
考虑还有什么时候会impossible，只有存在桥的时候，我们不能够给这条桥边定向使得新图强连通
于是tarjan求桥然后随意dfs

#include 
#include 
#include 
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

const int N=2000005;
const int E=2000005;

struct edge {int x,y,next;} e[E];
int ls[N],prt[N],edCnt=1;

void add_edge(int x,int y) {
    e[++edCnt]=(edge) {x,y,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {y,x,ls[y]}; ls[y]=edCnt;
}

int dfn[N],low[N],vis[N];

int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}

void dfs(int now,int fa) {
    dfn[now]=low[now]=++dfn[0];
    for (int i=ls[now];i;i=e[i].next) {
        if (!dfn[e[i].y]) {
            dfs(e[i].y,now);
            low[now]=std:: min(low[now],low[e[i].y]);
        } else if (e[i].y!=fa) {
            low[now]=std:: min(low[now],dfn[e[i].y]);
        }
    }
}

void solve(int now) {
    vis[now]=1;
    for (int i=ls[now];i;i=e[i].next) {
        prt[i/2]=i&1;
        if (vis[e[i].y]) continue;
        solve(e[i].y);
    }
}

int main(void) {
    freopen("data.in","r",stdin);
    freopen("myp.out","w",stdout);
    int n=read(),m=read();
    rep(i,1,m) {
        int x=read(),y=read();
        add_edge(x,y);
    }
    dfs(1,-1);
    if (dfn[0]!=n) {
        puts("impossible");
        return 0;
    } int cnt=0;
    for (int i=2;i<=edCnt;i+=2) {
        if (low[e[i].x]!=low[e[i].y]) {
            cnt++;
        }
    }
    if (!cnt) {
        solve(1);
        rep(i,1,m) printf("%d", prt[i]);
    } else puts("impossible");
    return 0;
}

F

---

我并没有做出来