[USACO15JAN]草鉴定Grass Cownoisseur 题解

本文解析了一道关于图论的竞赛题目,通过强连通分量缩点的方法将原图转化为无环图,并利用SPFA算法预处理出1到每个点及每个点到1的最大距离。最终通过枚举每条边的逆行来找出能经过最多的不同点的数量。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目大意:有n个点,有一只牛从1号点出发,最后要走回1号点,可以逆行一次,要使得路径上经过的不同的点最多,问最多能经过多少个不同的点。

首先可以发现,假如它走到了一个环里面,那么它一定是可以走完环内所有点的,并且可以从环内任意一个点出发继续向下走,那么首先可以先缩点,那么缩完点之后的图是一定没有环的,就是说你从一号点出发之后就一定走不回来了,那么这个时候就要用到逆行了,我们只要枚举缩完点后的图中的边,尝试着逆行这条边,看看最多能经过多少个点,记录最大值就好了,但是如果每一次都跑一遍spfa肯定超时,所以我们可以预处理出1到每个点的最大距离以及每个点到1的最大距离,那这样子每一次逆行所能经过的点数都可以O(1)求出。

代码如下:

#include <cstdio>
#include <cstring>

int dfn[100010],low[100010];
int belong[100010],cnt=0;
int n,m,len=0,lenn=0;
struct node{int x,y,next;};
node e[100010],ee[100010];
int first[100010];
void buildroad(int x,int y)//这是没缩点时的
{
    len++;
    e[len].x=x;
    e[len].y=y;
    e[len].next=first[x];
    first[x]=len;
}
void buildotherroad(int x,int y)//这是缩点后的
{
    lenn++;
    ee[lenn].x=x;
    ee[lenn].y=y;
    ee[lenn].next=first[x];
    first[x]=lenn;
}
int tot=0;
int zhan[100010],t=0;
bool v[100010];
int size[100010];
int minn(int x,int y){return x<y?x:y;}
void dfs(int x)//强连通缩点
{
    zhan[++t]=x;
    dfn[x]=low[x]=++tot;
    v[x]=true;
    for(int i=first[x];i;i=e[i].next)
    {
        int y=e[i].y;
        if(!dfn[y])
        {
            dfs(y);
            low[x]=minn(low[x],low[y]);
        }
        else if(v[y])low[x]=minn(low[x],dfn[y]);
    }
    if(dfn[x]==low[x])
    {
        int xx;cnt++;
        do
        {
            size[cnt]++;
            xx=zhan[t--];
            v[xx]=false;
            belong[xx]=cnt;//belong表示每个点属于哪个强连通分量
        }while(xx!=x);
    }
}
int q[100010];
int dis[2][100010];
void spfa(int id)
{
    int st=1,ed=2;
    q[st]=belong[1];
    memset(v,false,sizeof(v));
    v[q[st]]=true;
    dis[id][belong[1]]=0;//我这里初始化成0,也就是我少算了size[1],但是我会在下面加上
    					//如果这里改成size[1]的话,下面就要减掉size[1],想想为什么 
    while(st!=ed)
    {
        int x=q[st];
        for(int i=first[x];i;i=ee[i].next)
        {
            int y=ee[i].y;
            if(dis[id][y]<dis[id][x]+size[y])//注意这里要加size[y],size表示这个点是多少个点缩成的
            {
                dis[id][y]=dis[id][x]+size[y];
                if(!v[y])
                {
                	v[y]=true;
                	q[ed++]=y;
                }
            }
        }
        st++;
        v[x]=false;
    }
}
void go(int x)//预处理出1到每个点的最大距离以及每个点到1的最大距离,1到每个点的最大距离直接跑一次spfa,每个点到1的距离就把所有边反过来再跑一边就好了
{
    lenn=0;
    memset(first,0,sizeof(first));
    for(int i=1;i<=len;i++)
    if(belong[e[i].x]!=belong[e[i].y])
    {
        if(!x)buildotherroad(belong[e[i].x],belong[e[i].y]);
        else buildotherroad(belong[e[i].y],belong[e[i].x]);
    }
    spfa(x);
}
int main()
{
    scanf("%d %d",&n,&m);
    memset(first,0,sizeof(first));
    memset(dis,-1,sizeof(dis));
    memset(size,0,sizeof(size));
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d %d",&x,&y);
        buildroad(x,y);
    }
    memset(dfn,0,sizeof(dfn));
    memset(v,false,sizeof(v));
    for(int i=1;i<=n;i++)
    if(dfn[i]==0)dfs(i);
    go(1);go(0);//1是求每个点到1,0是1到每个点,注意go(1)和go(0)不能互换,因为go(0)跑的时候用的才是正常的图,而这幅图下面还用得到
    int ans=0;
    for(int i=1;i<=lenn;i++)
    {
        int x=ee[i].x,y=ee[i].y;//x,y原本表示这条边的起点和终点,但我们要把枚举的边反过来,所以在这里y是起点x是终点 
        if(dis[0][y]!=-1&&dis[1][x]!=-1&&dis[0][y]+dis[1][x]+size[belong[1]]>ans)//假如1能够到达y并且x能够到达1,就用他们路径上的总点数尝试更新ans 
        ans=dis[0][y]+dis[1][x]+size[belong[1]];
    }
    printf("%d",ans);
}

### USACO P2035 iCow 题目解析 #### 问题描述 Farmer John 购买了一台新的 MP3 播放器 iCow,其中存储了 N (1 ≤ N ≤ 1,000) 首歌曲,每首歌都有一个初始权值 Ri (1 ≤ Ri ≤ 10,000)[^4]。播放顺序由 FJ 设计的独特算法决定: - 下一首播放的是当前所有未播放过的歌曲中权值最高的那一首;若有多个最高权值,则选择编号最小的一首。 - 当某首歌曲播放完成后,其权值会均匀分配给其余 N − 1 首歌曲,并将其自身的权值设为零。 - 如果该歌曲的权值不能被 N − 1 整除,则剩余部分将以 1 单位的形式依次给予排名靠前但尚未获得额外权重的歌曲。 任务是求出按照上述规则最先播放的 T (1 ≤ T ≤ 1000) 首歌曲的具体情况。 #### 解决方案思路 为了模拟这个过程并找到最开始播放的 T 首歌曲,可以采用优先队列(最大堆)来管理待播列表及其对应的权重。每次取出具有最大权重的元素作为即将播放的对象,在更新其他成员的新权重之后重新加入到队列当中继续循环直至达到所需次数为止。 具体步骤如下: - 初始化数据结构:创建一个包含所有歌曲 ID 和它们各自起始分数的最大堆; - 输出此目标的信息; - 更新剩余项目的得分并将已处理项移回至集合内等待下次轮转; 下面给出完整的 C++ 实现代码示例: ```cpp #include <iostream> #include <queue> using namespace std; struct Song { int id; long score; }; bool operator<(const Song& a, const Song& b){ return !(a.score > b.score || (a.score == b.score && a.id < b.id)); } int main(){ priority_queue<Song> pq; int n,t,r; cin>>n>>t; for(int i=1;i<=n;++i){ cin >> r; pq.push({i,r}); } while(t--){ auto top_song=pq.top(); cout<<top_song.id<<"\n"; vector<int> remainders; pq.pop(); if(top_song.score%(n-1)!=0){ for(int j=0;j<top_song.score%(n-1);++j) remainders.push_back(j); } while(!pq.empty()){ Song current = pq.top(); pq.pop(); current.score += top_song.score/(n-1); if (!remainders.empty()) { current.score++; remainders.erase(remainders.begin()); } pq.push(current); } // Reinsert the played song with zero points back into queue. pq.push({top_song.id, 0}); } } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值