洛谷八连测 #R7

LIST

题目链接

T1

T1 每次spfa(或者bfs)，时间复杂度O(spfa*Q)；
每次先把每一个查询的特殊点入队，距离为0；
spfa或者bfs即可。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<ctime>
#define LL long long
using namespace std;
const int N=5009;
int head[N],nxt[2*N],to[2*N],tot;
int n,m,q,k,ans;
int dis[N];
bool vis[N];
queue <int> qu;
void add(int x,int y)
{
    to[++tot]=y;
    nxt[tot]=head[x];
    head[x]=tot;
}
void spfa()
{
    while(!qu.empty())
    {
        int x=qu.front();
        qu.pop();
        vis[x]=0;
        for(int i=head[x];i;i=nxt[i])
        {
            if(dis[to[i]]>dis[x]+1&&dis[x]+1<=k)//小小优化
            {
                dis[to[i]]=dis[x]+1;
                if(!vis[to[i]])
                {
                    qu.push(to[i]);
                    vis[to[i]]=1;
                }
            }
        }
    }
}
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    scanf("%d%d%d",&n,&m,&q);
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y);add(y,x);   
    } 
    for(int i=1;i<=q;i++)
    {
        memset(dis,127/3,sizeof(dis));
        memset(vis,0,sizeof(vis));
        while(!qu.empty()) qu.pop();
        scanf("%d%d",&m,&k);
        for(int x,i=1;i<=m;i++)
        {
            scanf("%d",&x);
            dis[x]=0;
            qu.push(x);
            vis[x]=1;
        }
        spfa();
        ans=0;
        for(int i=1;i<=n;i++)
         if(dis[i]<=k) ans++;

        printf("%d\n",ans);
    }
    //cout<<clock();
    return 0;
} 

T2

T2
离散化+树状数组可以过掉60，时间复杂度O(n^2logn)；
100分的做法：
我们在树状数组的数组再加一个记录坐标和的值。
对于一个满足 i < j,a i > a j 的点对 (i,j)，
其会在满足 1 ≤ l ≤ i,j ≤ r ≤ n 的子区间 [l,r] 内出现。
所以这个点对对答案的贡献为 i × (n − j + 1)。
还需要高精度！
时间复杂度O(nlogn)

60分代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<ctime>
#define LL long long
using namespace std;
const int N=1e6+7;
int n,a[N],c[N];
LL ans,sum,b[N];
LL query(int x)
{
    LL s=0;
    for(int i=x+1;i<=n;i+=i&(-i))
    {
        s+=b[i];
    }
    return s;
}
void update(int x)
{
    for(int i=x;i;i-=i&(-i))
     b[i]++;
}
int main()
{
    freopen("b.in","r",stdin);
    freopen("b.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),c[i]=a[i];

    sort(c+1,c+n+1);
    int m=unique(c+1,c+n+1)-c-1; 
    for(int i=1;i<=n;i++) a[i]=lower_bound(c+1,c+m+1,a[i])-c;

    for(int s=1;s<=n;s++)
    {
        memset(b,0,sizeof(b));
        sum=0;
        for(int i=s;i<=n;i++)
        {
            sum+=query(a[i]);
            ans+=sum;
            update(a[i]);
        }
    }   
    printf("%lld\n",ans);
    //cout<<clock();
    return 0;
}

T3

Tarjan缩点+二分+拓扑排序+贪心

Tarjan缩点，然后处理出缩点后的拓扑序，然后按照拓扑排序 倒着 来一个个满足条件即可。倒着的原因是先满足儿子节点的需求，因为儿子不够可以向父亲要。

极易被卡常。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define LL long long
using namespace std;
const int N=1e6+77;
int n,m,f[N],F[N],a[N],cnt,b[N];
int head[N],nxt[2*N],to[2*N],tot;
int dfn[N],low[N],stack[N],color[N],top,dfs_num,color_num;
LL c[N],cc[N],num[N];
bool vis[N];
void add(int x,int y)
{
    to[++tot]=y;
    nxt[tot]=head[x];
    head[x]=tot;
}
void dfs(int x)//拓扑排序 
{
    b[++cnt]=x;
    for(int i=head[x];i;i=nxt[i])
     dfs(to[i]);
}
void DFS(int x)//Tarjan缩点 
{
    dfn[x]=++dfs_num;
    low[x]=dfs_num;
    stack[++top]=x;
    vis[x]=1;

    for(int i=head[x];i;i=nxt[i])
    {
        if(!dfn[to[i]])
        {
            DFS(to[i]);
            low[x]=min(low[x],low[to[i]]);
        }
        else
        if(vis[to[i]]) low[x]=min(low[x],dfn[to[i]]);
    }

    if(low[x]==dfn[x])
    {
        vis[x]=0;
        color[x]=++color_num;

        while(stack[top]!=x)
        {
            vis[stack[top]]=0;
            color[stack[top]]=color_num;
            top--;
        }
        top--;
    }
}
bool check(LL x)
{
    for(int i=1;i<=color_num;i++) cc[i]=c[i];
    LL last=m;

    for(int i=cnt;i>=1;i--)
    {
        int p=b[i],px=F[p];
        LL need=x*num[i];
        /*while(cc[p]<need&&px)
        {
            if(cc[px]>=need-cc[p])
            {
                cc[px]-=need-cc[p];
                cc[p]=need;
            } 
            else
            {
                cc[p]+=cc[px];
                cc[px]=0;
            }
            px=F[px];
        }
        if(cc[p]<need)
        {
            if(last>=need-cc[p])
            {
                last-=need-cc[p];
                cc[p]=need;
            } 
            else return 0;
        }
        if(cc[p]<need) return 0;*/
        if(need<=cc[p]) continue; 
        if(!px) {
            if(last<need-cc[p]) return 0;
            else 
            {
                last-=need-cc[p];
                cc[p]=need;
            }
        }
        if(cc[p]<need)
        {
            cc[px]-=need-cc[p];
            cc[p]=need;
        }
    }
    return 1;
}
int main()
{
    //freopen("c2.in","r",stdin);
    //freopen("c.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) 
    {
        scanf("%d",&f[i]);
        if(f[i]!=-1) add(f[i],i);
    }
    for(int i=1;i<=n;i++) scanf("%d",&a[i]); 

    for(int i=1;i<=n;i++) if(!dfn[i]) DFS(i);

    memset(head,0,sizeof(head));
    memset(to,0,sizeof(to));
    memset(nxt,0,sizeof(nxt));
    tot=0;

    for(int i=1;i<=n;i++)
    {
        if(color[i]!=color[f[i]])
        {
            add(color[f[i]],color[i]);
            F[color[i]]=color[f[i]];
        }//缩点 
        num[color[i]]++;
        c[color[i]]+=a[i];
    }

    for(int i=1;i<=color_num;i++) if(!F[i]) dfs(i);//拓扑排序 


    LL L=0,R=2e9+7,mid,ans=L;
    while(L<=R)
    {
        mid=(L+R)>>1;
        if(check(mid)) ans=mid,L=mid+1;
        else R=mid-1;
    }
    printf("%lld\n",ans);
    return 0;
}