100211D Police Cities

传送门

分析

看到这个题我们的第一反应自然是Tarjan缩点，在这之后我们可以发现实际只要在缩点之后所有出度或入度为0的点布置警察局就可以达到要求，我们用dpij表示考虑前i个出度或入度为0的点共布置了j个警察局，s[i]表示这个点原先由几个点构成，tot表示出度或入度为0的点的总数，所以不难得到转移方程

所以不难得到最终答案

注意此题需要使用高精度。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
#define sp cout<<"---------------------------------------------------"<<endl
struct mint {
      int _[],__;
};
mint operator + (mint _x,mint _y){
      int _k,_g=,_i;
      mint _z;
      for(_i=_x.__+;_i<=;_i++)_x._[_i]=;
      for(_i=_y.__+;_i<=;_i++)_y._[_i]=;
      if(_x.__>_y.__)_k=_x.__;
        else _k=_y.__;
      for(_i=;_i<=_k;_i++){
          _z._[_i]=(_x._[_i]+_y._[_i]+_g)%;
          _g=(_x._[_i]+_y._[_i]+_g)/;
      }
      if(_g>){
          _z._[++_k]=_g;
      }
      _z.__=_k;
      return _z;
}
mint operator * (mint _x,mint _y){
      int _k,_g=,_i,_j;
      mint _z;
      for(_i=_x.__+;_i<=;_i++)_x._[_i]=;
      for(_i=_y.__+;_i<=;_i++)_y._[_i]=;
      _k=_x.__+_y.__-;
      for(_i=;_i<=;_i++)
        _z._[_i]=;
      for(_i=;_i<=_x.__;_i++)
        for(_j=;_j<=_y.__;_j++)
          _z._[_i+_j-]+=_x._[_i]*_y._[_j];
      for(_i=;_i<=_k;_i++){
          int _a=_z._[_i]+_g;
          _z._[_i]=_a%;
          _g=_a/;
      }
      while(_g){
        _z._[++_k]=_g%;
        _g/=;
      }
      while(_k>&&_z._[_k]==)_k--;
      _z.__=_k;
      return _z;
}
mint read(){
      mint _x;
      string _s;
      int _L,_i;
      for(_i=;_i<=;_i++)
        _x._[_i]=;
      cin>>_s;
      _L=_s.length();
      for(_i=;_i<=_L;_i++)
        _x._[_i]=_s[_L-_i]-'';
      _x.__=_L;
      return _x;
}
void pr(mint _x){
      int _i;
      for(_i=_x.__;_i>;_i--)
        printf("%d",_x._[_i]);
      puts("");
}
void test(){
      mint _a,_b;
      _a=read(),_b=read();
      _a=_a+_b;
      pr(_a);
}
int n,m,K,sum,s[],dfn[],low[],ist[],cnt,belong[];
int tot,o_d[],i_d[],wh[];
mint dp[][],c[][];
stack<int>a;
vector<int>o_v[];
vector<int>v[];
inline void tarjan(int x){
      dfn[x]=low[x]=++cnt;
      a.push(x);
      ist[x]=;
      for(int i=;i<o_v[x].size();i++)
        if(!dfn[o_v[x][i]]){
            tarjan(o_v[x][i]);
            low[x]=min(low[x],low[o_v[x][i]]);
        }else if(ist[o_v[x][i]]){
            low[x]=min(low[x],dfn[o_v[x][i]]);
        }
      if(low[x]==dfn[x]){
          sum++;
        while(){
          int u=a.top();
          a.pop();
          ist[u]=;
          s[sum]++;
          belong[u]=sum;
          if(u==x)break;
        }
      }
}
int main(){
      freopen("police.in","r",stdin);
      freopen("police.out","w",stdout);
      int i,j,k;
      //test();
      for(i=;i<=;i++)
        for(j=;j<=;j++){
          c[i][j]._[]=;
          c[i][j].__=;
        }
      for(i=;i<=;i++){
        c[i][]._[]=;
        c[i][].__=;
        c[i][i]._[]=;
        c[i][i].__=;
      }
      for(i=;i<=;i++)
        for(j=;j<i;j++)
          c[i][j]=c[i-][j-]+c[i-][j];
      scanf("%d%d%d",&n,&m,&K);
      for(i=;i<=m;i++){
          int x,y;
          scanf("%d%d",&x,&y);
          o_v[x].push_back(y);
      }
      for(i=;i<=n;i++)
        if(!dfn[i])tarjan(i);
      for(i=;i<=n;i++)
        for(j=;j<o_v[i].size();j++)
          if(belong[i]!=belong[o_v[i][j]]){
              v[belong[i]].push_back(belong[o_v[i][j]]);
              o_d[belong[i]]++;
              i_d[belong[o_v[i][j]]]++;
          }
      int ant=;
      for(i=;i<=sum;i++)
        if(!o_d[i]||!i_d[i]){
          tot++;
          wh[tot]=i;
        }else ant+=s[i];
      for(i=;i<=;i++)
        for(j=;j<=;j++){
          dp[i][j].__=;
          dp[i][j]._[]=;
        }
      dp[][]._[]=;
      dp[][].__=;
      for(i=;i<=tot;i++)
        for(j=;j<=K;j++)
          for(k=;k<=s[wh[i]];k++)
            dp[i][j]=dp[i][j]+(dp[i-][j-k]*c[s[wh[i]]][k]);
      mint ans;
      ans.__=;
      ans._[]=;
      for(i=tot;i<=K;i++)
        ans=ans+(dp[tot][i]*c[ant][K-i]);
      pr(ans);
      return ;
}