【ContestHunter】【弱省胡策】【Round6】

KMP/DP+树链剖分+线段树/暴力

---

　　今天考的真是……大起大落……

String

　　QwQ题意理解又出错了……（还是说一开始理解了，后来自己又忘了为什么是这样了？）

　　反正最后的结果就是……我当成：后面每行只需要和第一行check一下就可以了。

　　因为那个图真的很像在搞串的匹配啊……一格一格往过移，看能不能匹配

　　然后不就是每个s[i]对A取一下模，然后KMP？看有多少个位置能匹配咯。。。

　　（其实这是A=B时的做法）

　　将原序列增长一倍（复制一遍），对A取模做一遍，然后再对B取模做一遍，看有哪些位置是合法的，统计一下……

　　然而其实正解可以证明出来当$A ot = B$的时候，可以当A=B做……也就是说，上面那种做法……是可以AC的-_-b

　　P.S.正解是KMP来找循环节= = n如果可以整除n-next[n]，那么循环节长度就是n-next[n]……可匹配长度就是……那么多……啊反正差不多啦

//Round6 A
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
    int r=1,v=0; char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-1;
    for(; isdigit(ch);ch=getchar()) v=v*10-'0'+ch;
    return r*v;
}
const int N=2000010;
/*******************template********************/


int a[N],b[N],c[N],A,B,n,ans,nxt[N];
bool yes[N];
void init(){
    n=getint(); A=getint(); B=getint();
    F(i,1,n) a[i]=getint();
    memset(yes,0,sizeof yes);
    ans=0;
}
void KMP(int *s,int n){
    int m=n*2;
    nxt[1]=0; int j=0;
    F(i,2,n+1){
        while(j && s[i]!=s[j+1]) j=nxt[j];
        if (s[i]==s[j+1]) j++;
        nxt[i]=j;
    }
    j=0;
    F(i,2,m-1){
        while(j && b[i]!=s[j+1]) j=nxt[j];
        if (b[i]==s[j+1]) j++;
        if (j==n) yes[i]=1,j=nxt[j];
    }
}
void work(){
    F(i,1,n) c[i]=b[i]=b[i+n]=a[i]%A;
    KMP(c,n);
    if (A!=B){
        F(i,1,n) c[i]=b[i]=b[i+n]=a[i]%B;
        KMP(c,n);
    }
    F(i,n,n*2-1) ans+=yes[i];
    printf("%d
",ans?ans+1:0);
}
int main(){
    int T=getint();
    while(T--){
        init();
        work();
    }
    return 0;
}

Tree

　　……30分的做法其实就是直接TreeDP，40~50的做法是线段树维护最大连续子段和，那么满分做法？两者结合……Orzzzzzz

　　所以就是树链剖分+DP（并没有？其实是贪心？）啊……

　　Orz zyf，这题代码感觉细节很多啊……

　　树链剖分一下，权值沿轻链按TreeDP的思路将权值加到fa上面，建出整个线段树……

　　然后每次修改的时候，先将这个点x到top[x]这一条重链上的ans更新一下（这里重点是maxl，即从左端点开始的最大连续子段和），然后！更新fa[top[x]]的DP值！因为它在线段树上的权值是自身权值+轻链的maxl，怎么搞？区间查询一下top[x]所在的重链的maxl……然后更新。。。【这里我并不知道如何记重链的结尾……因为dfs序记start和end是子树的，所以并不能那样搞。。。zyf的做法是：往出连轻链的时候多++一下dfs_clock，这样就可以将不同的重链断开，（当然是将这些没有编号的空节点的权值置为-INF）

　　查询的时候……直接出解呗。。。

//Round6 B
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
    int r=1,v=0; char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-1;
    for(; isdigit(ch);ch=getchar()) v=v*10-'0'+ch;
    return r*v;
}
const int N=500010,INF=1000000000;
/*******************template********************/

int head[N],to[N],nxt[N],cnt;
void add(int x,int y){
    to[++cnt]=y; nxt[cnt]=head[x]; head[x]=cnt;
    to[++cnt]=x; nxt[cnt]=head[y]; head[y]=cnt;
}

int n,m,f[N],id[N],v[N],ans;
int pos[N],dfs_clock,fa[N],son[N],size[N],top[N];
struct node{
   int max,maxl,maxr;
   LL sum;
}t[N];
bool vis[N];

void dfs(int x){
    size[x]=1; son[x]=0;
    int mx=0;
    for(int i=head[x];i;i=nxt[i])
        if (to[i]!=fa[x]){
            fa[to[i]]=x;
            dfs(to[i]);
            size[x]+=size[to[i]];
            if (size[to[i]]>mx) son[x]=to[i],mx=size[to[i]];
        }
}
void connect(int x,int f){
    pos[x]=++m; top[x]=f; vis[x]=1;
    id[m]=x;
    if (son[x]) connect(son[x],f);
    for(int i=head[x];i;i=nxt[i])
        if (!vis[to[i]])
            ++m,connect(to[i],to[i]);
}

#define L (o<<1)
#define R (o<<1|1)
#define mid (l+r>>1)
node maintain(node a,node b){
    node ans;
    ans.maxl=max((LL)a.maxl,a.sum+b.maxl);
    ans.maxr=max((LL)b.maxr,b.sum+a.maxr);
    ans.max=max(a.maxr+b.maxl,max(a.max,b.max));
    ans.sum=a.sum+b.sum;
    return ans;
}
void update(int o,int l,int r,int p,int v){
    if (l==r){
        t[o].sum+=v;
        t[o].maxl=t[o].maxr=t[o].max=max(t[o].sum,(LL)0);
    }else{
        if (p<=mid) update(L,l,mid,p,v);
        else update(R,mid+1,r,p,v);
        t[o]=maintain(t[L],t[R]);
    }
}
node query(int o,int l,int r,int ql,int qr){
    if (l==ql && r==qr) return t[o];
    else if (qr<=mid) return query(L,l,mid,ql,qr);
    else if (ql>mid) return query(R,mid+1,r,ql,qr);
    else return maintain(query(L,l,mid,ql,mid),query(R,mid+1,r,mid+1,qr));
}
void change(int x,int y1,int y2){
    while(x){
        int y=top[x],z=query(1,1,m,pos[y],m).maxl;
        update(1,1,m,pos[x],y2-y1);
        y1=z; y2=query(1,1,m,pos[y],m).maxl;
        x=fa[top[x]];
    }
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("sub.in","r",stdin);
    freopen("sub.out","w",stdout);
#endif
    n=getint(); int T=getint();
    F(i,1,n) v[i]=getint();
    F(i,2,n){
        int x=getint(),y=getint();
        add(x,y);
    }
    dfs(1); connect(1,1);
    F(i,1,m) if (!id[i]) update(1,1,m,i,-INF);
    F(i,1,n) change(i,0,v[i]);
    while(T--){
        int cmd=getint(),x,y;
        if (cmd==1){
            x=getint(), y=getint();
            change(x,v[x],y);
            v[x]=y;
        }else printf("%d
",t[1].max);
    }
    return 0;
}

Transport

　　……我就无语了……

　　一开始我想的是二分！然后模拟过程进行判断！然而连样例都过不了……QwQ

　　仔细想了想好像确实并不满足单调性……会有小的波动……

　　然后蒟蒻就不会了……实在没办法了，二分的时候边界不是 $L=max{a[i]}，R=sum a[i]$ 嘛= =然后我就直接将二分的过程删掉……改成枚举……然后……TM居然就过了……maya我也是惊呆了，这个枚举的上下界明明很大的好吗……如果全都走一遍肯定要T的啊！然而我就是过了……就是这么神奇……>_>

　　P.S.这个边界 L 还是比较松的，其实应该弄成$max{ frac{sum a[i]}{k}, max(a[i]) }$ ？

//Round6 C
#include<cstdio>
#include<set>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
    int r=1,v=0; char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-1;
    for(; isdigit(ch);ch=getchar()) v=v*10-'0'+ch;
    return r*v;
}
const int N=100010;
/*******************template********************/

int n,k,a[N],sum;
multiset<int>s,s2;
multiset<int>::iterator it;
bool check(int c){
    if (k*c<sum) return 0;
    s2=s;
    int num;
    F(i,1,k){
        num=c;
        do{
            if (s2.empty()) return 1;
            it=s2.lower_bound(num);
            if (it==s2.end()) it--;
            if (*it>num){
                if (it==s2.begin()) break;
                else it--;
            }
            num-=*it;
            s2.erase(it);
        }while(1);
    }
    return s2.empty();
}
void init(){
    n=getint(); k=getint();
    s.clear();
    int L=0,R=0;
    sum=0;
    F(i,1,n){
        a[i]=getint();
        sum+=a[i]; R+=a[i]; L=max(L,a[i]);
        s.insert(a[i]);
    }
    F(i,L,R)
        if (check(i)){
            printf("%d
",i);
            return;
        }
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("C.in","r",stdin);
    freopen("C.out","w",stdout);
#endif
    int T=getint();
    while(T--) init();
    return 0;
}