套题 codeforces 359

A题：Free Ice Cream

注意要使用LL,避免爆int

#include <bits/stdc++.h>
#define scan(x,y) scanf("%d%d",&x,&y)
using namespace std;
typedef long long LL;
const int Max=1e5+;
int main()
{
    int n,has,x;
    while(~scan(n,has))
    {
        LL have=has;
        char op;LL dis=;
        for(int i=;i<n;i++)
        {
            cin>>op>>x;
            if(op=='+')
            {
                have+=x;
            }
            else
            {
                if(have>=x) have-=x;
                else  dis++;
            }
        }
        printf("%I64d %I64d\n",have,dis);
    }
}

B题：Little Robber Girl's Zoo
不要想得太复杂，直接暴力就好

#include <bits/stdc++.h>
#define scan(x,y) scanf("%d%d",&x,&y)
using namespace std;
typedef long long LL;
const int Max=1e3+;
int a[Max][Max],L[Max],R[Max],best[Max],l,r;
int g=,n;
bool check()
{
    int i=,cnt=,cunt=;
    memset(L,,sizeof(L));memset(R,,sizeof(R));
    while(i<n)
    {
        if(best[i]>best[i+])
        {
            cnt++;
            swap(best[i],best[i+]);
            if(cnt==) {L[cunt]=i;R[cunt++]=i+;}
            else
            {
                if(i-==R[cunt-]) R[cunt-]=i+;
                else {L[cunt]=i;R[cunt++]=i+;}
            }
            i+=;
        }
        else
        {
            i++;
        }
    }
    if(cnt==) return ;
    for(int i=;i<cunt;i++) printf("%d %d\n",L[i],R[i]);
    return ;
}
int main()
{
    while(~scanf("%d",&n))
    {
        g=;
        for(int i=;i<=n;i++) scanf("%d",&best[i]);
        memset(a,,sizeof(a));
        while(!check());
    }
    return ;
}

C题：Robbers' watch
观察可知，由于时间表是7进制的,所以整个可行的序列长不会超过7位，因此可以套用数位dp的模式进行dfs搜索

#include <bits/stdc++.h>
using namespace std;
int ans,n,m;
int vis[];
int bit[][];
int pos[];
stack<int>st1,st;
void dfs2(int len,int flag,int id)
{
    if(len==)
    {
        ans++;
        return;
    }
    int end=flag?bit[id][len]:;
    for(int i=;i<=end;i++)
    {
        if(!vis[i])
        {
            vis[i]=;
            dfs2(len-,i==end&&flag,id);
            vis[i]=;
        }
    }
}
void dfs(int len,int flag,int id)
{
    if(len==)
    {
        if(id==) dfs2(pos[],,);
        return;
    }
    int end=flag?bit[id][len]:;
    for(int i=;i<=end;i++)
    {
        if(!vis[i])
        {
            vis[i]=;
            dfs(len-,i==end&&flag,id);
            vis[i]=;
        }
    }
}
int cacu(int x,int y)
{
    pos[]=;
    bit[][]=bit[][]=;
    if(x<) bit[][++pos[]]=x;
    else
    {
        while(x)
        {
            bit[][++pos[]]=x%;
            x/=;
        }
    }
    pos[]=;
    if(y<) bit[][++pos[]]=y;
    else
    {
        while(y)
        {
            bit[][++pos[]]=y%;
            y/=;
        }
    }
    if(pos[]+pos[]>) return ;
    else dfs(pos[],,);
    return ;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        ans=;
        n--;m--;
        printf("%d\n",cacu(n,m)==?:ans);
    }
    return ;
}

D题：Kay and Snowflake
1.重心的另一个定义：以这个点为根，那么所有的子树（不算整个树自身）
的大小都不超过整个树大小的一半。（这个对解这道题很有用）
题解思路：
对于以u为根节点的子树，重心肯定在它的重儿子之中，但是有一个范围
可以想到重心的下区间是重儿子的重心，上区间就是u自己
可以对于这个范围内的节点，利用性质1逐个检测。

#include <bits/stdc++.h>
using namespace std;
const int Max=3e5+;
const int E=Max*;
int son[Max],maxson[Max],fa[Max],zx[Max];
int head[Max],nex[E],pnt[E],edge;
void Init()
{
    edge=;
    memset(son,,sizeof(son));
    memset(maxson,,sizeof(maxson));
    memset(fa,-,sizeof(fa));
    memset(head,-,sizeof(head));
}
void Addedge(int u,int v)
{
    pnt[edge]=v;
    nex[edge]=head[u];
    head[u]=edge++;
}
void dfs1(int u)
{
    son[u]=;
    maxson[u]=;
    int v;
    for(int e=head[u];e!=-;e=nex[e])
    {
        v=pnt[e];
        dfs1(v);
        son[u]+=son[v];
        maxson[u]=max(maxson[u],son[v]);
    }
}
bool is_ok(int u,int v)
{
    if(son[u]>=maxson[v]*&&son[u]>=(son[u]-son[v])*) return ;
    //利用定义4，检测是否是正确的重心，存在u自身的情况
    //所以不能son[u]>=son[v]*2;
    return ;
}
void dfs2(int u)
{
    if(son[u]==)
    {
        zx[u]=u;
        return;
    }
    int v,z=u;
    for(int e=head[u];e!=-;e=nex[e])
    {
        v=pnt[e];
        dfs2(v);
        if(maxson[u]==son[v])
        {
            z=zx[v];
        }
    }
    while(!is_ok(u,z)&&fa[z]!=-)
    {
        z=fa[z];
    }
    zx[u]=z;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
       Init();
       int u,v;
       for(int i=;i<=n;i++)
       {
           scanf("%d",&u);
           fa[i]=u;
           Addedge(u,i);
       }
       dfs1();
       dfs2();
       while(m--)
       {
           scanf("%d",&u);
           printf("%d\n",zx[u]);
       }
    }
    return ;
}

E题：Optimal Point
题解链接：http://codeforces.com/blog/entry/45558

#include <bits/stdc++.h>
using namespace std;
#define MP(x,y) make_pair(x,y);
typedef long long LL;
const int Max=1e5+;
const LL inf=*(LL)1e9*(LL)1e9+;
//(1e18不能直接写，会爆int)
const LL LM=(LL()<<)*(LL()<<)*(LL()<<)-;
const LL LMM=-(LL()<<)*(LL()<<)*(LL()<<)+;
LL max(LL x,LL y){return x>=y?x:y;}
LL min(LL x,LL y){return x<=y?x:y;}
int n;
struct coord{
     LL x,y,z;
     coord(){};
     coord(LL xx,LL yy,LL zz):x(xx),y(yy),z(zz){}
}node[Max];
struct equa{
     pair<LL,LL> s;
     pair<LL,LL> a;
     pair<LL,LL> b;
     pair<LL,LL> c;
     equa operator + (const equa e)
     {
         equa ans;
         ans.s=MP(max(s.first,e.s.first),min(s.second,e.s.second));
         ans.a=MP(max(a.first,e.a.first),min(a.second,e.a.second));
         ans.b=MP(max(b.first,e.b.first),min(b.second,e.b.second));
         ans.c=MP(max(c.first,e.c.first),min(c.second,e.c.second));
         return ans;
     }
};
void Acc()
{
    std::iostream::sync_with_stdio(false);
    cin.tie(NULL);cout.tie(NULL);
}
LL Getmid(LL x)
{
    return (x-(x&))/;
}

struct coord Getsolu(const equa &x)
{
    if(x.s.first>x.s.second||
       x.a.first>x.a.second||
       x.b.first>x.b.second||
       x.c.first>x.c.second) return coord(LM,LM,LM);
    if(x.a.first+x.b.first+x.c.first>x.s.second||
       x.a.second+x.b.second+x.c.second<x.s.first)
       return coord(LM,LM,LM);

     coord res;
     res.x=x.a.first;
     res.y=x.b.first;
     res.z=x.c.first;
     LL delta=max(LL(),x.s.first-res.x-res.y-res.z);
     res.x+=min(delta,x.a.second-x.a.first);
     delta-=min(delta,x.a.second-x.a.first);
     res.y+=min(delta,x.b.second-x.b.first);
     delta-=min(delta,x.b.second-x.b.first);
     res.z+=min(delta,x.c.second-x.c.first);
     delta-=min(delta,x.c.second-x.c.first);
     assert(delta==);
     return res;
}
struct coord Cacu(LL Maxlen)
{
    //设Maxlen为最大距离
    equa eq;
    eq.s=eq.a=eq.b=eq.c=MP(LMM,LM);
    for(int i=;i<=n;i++)
    {
        equa ne;
        //     -Maxlen<=A+B+C<=Maxlen
        //     -Maxlen<=A<=Maxlen
        //     -Maxlen<=B<=Maxlen
        //     -Maxlen<=C<=Maxlen
        ne.s=MP(node[i].x+node[i].y+node[i].z-Maxlen,
                node[i].x+node[i].y+node[i].z+Maxlen);
        ne.a=MP(-node[i].x+node[i].y+node[i].z-Maxlen,
                -node[i].x+node[i].y+node[i].z+Maxlen);
        ne.b=MP(node[i].x-node[i].y+node[i].z-Maxlen,
                node[i].x-node[i].y+node[i].z+Maxlen);
        ne.c=MP(node[i].x+node[i].y-node[i].z-Maxlen,
                node[i].x+node[i].y-node[i].z+Maxlen);
        eq=eq+ne;  //取最大距离
    }
    for(LL r=;r<=;r++)
    {
        //对二取余，可能余1或者余0
        equa tr=eq;
        tr.s.first=Getmid(tr.s.first-(*r-));
        tr.a.first=Getmid(tr.a.first-(r-));
        tr.b.first=Getmid(tr.b.first-(r-));
        tr.c.first=Getmid(tr.c.first-(r-));

        tr.s.second=Getmid(tr.s.second-(*r));
        tr.a.second=Getmid(tr.a.second-(r));
        tr.b.second=Getmid(tr.b.second-(r));
        tr.c.second=Getmid(tr.c.second-(r));

        coord solu=Getsolu(tr);
        if(solu.x!=LM)
        {
            coord ans;
            //(A,B奇偶性相同才会有整数解)
            ans.x=r+solu.y+solu.z;
            ans.y=r+solu.x+solu.z;
            ans.z=r+solu.x+solu.y;
            return ans;
        }
    }
    return coord(LM,LM,LM);  //有超过Maxlen的距离或者无整数解
}
int main()
{
    Acc();
    int T;
    for(cin>>T;T;T--)
    {
        cin>>n;
        for(int i=;i<=n;i++)
        {
            cin>>node[i].x>>node[i].y>>node[i].z;
        }
        LL l=,r=LM,answ=;
        while(l<=r)
        {
            LL mid=l+(r-l)/;
            if(Cacu(mid).x!=LM) r=mid-,answ=mid;
            else l=mid+;
        }
        coord ans=Cacu(answ);
        cout<<ans.x<<" "<<ans.y<<" "<<ans.z<<endl;
    }
    return ;
}