BZOJ 1124 [POI2008] 枪战Maf 乱搞。。

先考虑特殊情况

1、只有独立的一个点,即自杀:必死,最多最少都为1


2、一个独立的强联通分量,里面点数为num:最多num-1,从任意一个点倒着杀。最少为(num+1)/2,从任意一个点正向杀,隔一个杀一个。


3、非独立的强联通分量:入度为零的点是一定不死的,最多为总点数-入度为零点数。

      最困难的是找这种情况下的最少数,因为入度为零的点是一定不死,所以把它们都放入一个队列,

      每次找队首元素的目标,如果这个目标没有死,就把它标记为死去,同时把这个已死去点的目标的入度减1,如果目标入度减为零,将它入队。

      最后剩下的一定是许多环,再按照2求。


#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
queue<int> q;
struct bian
{
  int qi,zhong,next;
};
bian c[2000005];
int n,a[1000005],ji=0,duo=0,shao=0,jishu=0,jishu1=0,jishu2=0,jishu3=0,jishu4=0;
bool in[1000005],flag[1000005],si[1000005],biao[1000005];
int ru[1000005];
int sum,ji1=0;
int head[1000005],jilu[1000005],dfn[1000005],num[1000005],low[1000005],belong[1000005],zhan[1000005];
void add(int x,int y)
{
  c[++jishu].qi=x;
  c[jishu].zhong=y;
  c[jishu].next=head[x];
  head[x]=jishu;
}
void tarjan(int x)
{
  low[x]=++jishu2;
  dfn[x]=jishu2;
  zhan[++jishu3]=x;
  in[x]=1;
  for(int i=head[x];i;i=c[i].next)
  {
    if(dfn[c[i].zhong]==-1)
    {
      tarjan(c[i].zhong);
      low[x]=min(low[x],low[c[i].zhong]);
    }
    else  if(in[c[i].zhong])
      low[x]=min(low[x],dfn[c[i].zhong]);
  }
  if(low[x]==dfn[x])
  {
    int y;
    ji++;
    jishu4=0;
    while(1)
    {
      y=zhan[jishu3--];
      in[y]=0;
      jishu4++;
      belong[y]=ji;
      if(y==x)
        break;
    }
    num[ji]=jishu4;
  }
}
int main()
{
  freopen("maf.in","r",stdin);
  freopen("maf.out","w",stdout);
  memset(dfn,-1,sizeof(dfn));
  scanf("%d",&n);
  sum=n;
  for(int i=1;i<=n;++i)
  {
    scanf("%d",&a[i]);
    add(i,a[i]);
    ru[a[i]]++;
  }
  for(int i=1;i<=n;++i)
    if(dfn[i]==-1)
      tarjan(i);
  for(int i=1;i<=jishu;++i)
    if(c[i].qi==c[i].zhong)
    {
      duo++;
      shao++;
      sum--;
      flag[belong[c[i].qi]]=1;
      si[c[i].qi]=1;
    }
  for(int i=1;i<=jishu;++i)
    if(belong[c[i].qi]!=belong[c[i].zhong])
    {
      flag[belong[c[i].qi]]=1;
      flag[belong[c[i].zhong]]=1;
    }
  for(int i=1;i<=ji;++i)
    if(!flag[i])
    {
      duo+=num[i]-1;
      shao+=(num[i]+1)/2;
      sum-=num[i];
    }
  for(int i=1;i<=n;++i)
    if(ru[i]==0)
      ji1++;
  duo+=sum-ji1;
  for(int i=1;i<=n;++i)
    if(ru[i]==0)
      q.push(i);
  while(!q.empty())
  {
    int k=q.front();
    q.pop();
    if(!si[a[k]]&&!si[k])
    {
      si[a[k]]=1;
      biao[belong[a[k]]]=1;
      shao++;
      int e=a[k];
      ru[a[e]]--;
      if(ru[a[e]]==0&&!si[a[e]])
        q.push(a[e]);
    }
  }
  for(int i=1;i<=ji;++i)
    if(flag[i]&&num[i]>1&&!biao[i])
      shao+=(num[i]+1)/2;
  cout<<shao<<" "<<duo;
  return 0;
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值