@description@

给定一个长度为 N 的序列 a,问有多少排列 p,满足对于每一个 i,都有 \(a_i = p_i\) 或 \(a_i = p_{p_i}\) 成立。

原题传送门。

@solution@

为了更直观地理解问题,不妨建个图:连有向边 \(i -> p_i\)。

对于任意序列 a,这样建出的图是一个基环树森林;对于排列 p,这样建出的图是若干个环。

当 \(a_i = p_i\) 或 \(a_i = p_{p_i}\) 成立时,相当于我们把排列 p 对应的图的某个边 \(i -> p_i\) 替换成 \(i -> p_{p_i}\),最终得到序列 a 对应的图。

如果将某一个奇环的边全部替换,将得到另一个大小相同的奇环。

如果将某一个偶环的边全部替换,将得到两个大小为原先一半的环。

如果将某一个环的边不完全替换(不包括不替换的情况),一定得到一个基环树。

那么我们可以将 a 序列中相同大小的环放在一起计数,把 a 序列中每一个基环树进行计数。

环的情况一定是全部替换/全部不替换,上面的讨论已经包含了所有情况。

考虑基环树。首先你自己手玩一下,发现基环树如果要合法,一定是一个环 + 某些环上的点延伸出去恰好一条链。

其实比较容易理解。替换后每个点的入度最多为 2,且入度为 2 的一定是环上的点。

理论上延伸出去一条链的环上的点会贡献 2 倍方案数:可能延伸出去的链是被替换的;可能环上的入边是被替换的。但是有不合法的情况。

不妨记延伸出去链长为 x,顺着入边走找到的最近的有链延伸出去的点距离为 y。

当 x < y 时,延伸出去的链才可能是被替换的;当 x <= y 时,环上的入边才可能是被替换的。

这个自己手玩一下就能发现了。

@accepted code@

#include <cstdio>
#include <algorithm>
using namespace std; #define ILLEGAL puts("0"), exit(0) const int MAXN = 100000;
const int MOD = int(1E9) + 7;
const int INV2 = (MOD + 1) / 2; int pow_mod(int b, int p) {
int ret = 1;
for(int i=p;i;i>>=1,b=1LL*b*b%MOD)
if( i & 1 ) ret = 1LL*ret*b%MOD;
return ret;
} int fct[MAXN + 5], ifct[MAXN + 5], f[MAXN + 5];
int comb(int n, int m) {
return 1LL*fct[n]*ifct[m]%MOD*ifct[n-m]%MOD;
}
void init() {
fct[0] = 1;
for(int i=1;i<=MAXN;i++)
fct[i] = 1LL*fct[i-1]*i%MOD;
ifct[MAXN] = pow_mod(fct[MAXN], MOD - 2);
for(int i=MAXN-1;i>=0;i--)
ifct[i] = 1LL*ifct[i+1]*(i+1)%MOD;
f[0] = 1;
for(int i=2;i<=MAXN;i+=2)
f[i] = 1LL*comb(i, 2)*f[i-2]%MOD*pow_mod(i/2, MOD-2)%MOD;
} int a[MAXN + 5], N;
bool tag[MAXN + 5], vis[MAXN + 5];
int cnt[MAXN + 5], b[MAXN + 5], ind[MAXN + 5]; int main() {
init(), scanf("%d", &N);
for(int i=1;i<=N;i++)
scanf("%d", &a[i]), ind[a[i]]++;
for(int i=1;i<=N;i++)
if( ind[i] >= 3 ) ILLEGAL;
for(int i=1;i<=N;i++) {
if( ind[i] == 2 ) {
if( tag[i] ) continue;
int p = i;
while( !vis[p] )
vis[p] = true, p = a[p];
if( p != i ) ILLEGAL;
p = i;
do {
tag[p] = true, p = a[p];
}while( p != i );
}
}
for(int i=1;i<=N;i++) {
if( ind[i] == 0 ) {
int t = 0, p = i;
while( !tag[p] )
vis[p] = true, t++, p = a[p];
b[p] = t;
}
}
int ans = 1;
for(int i=1;i<=N;i++) {
if( b[i] ) {
int t = 0, p = i;
do {
t++, p = a[p];
}while( !b[p] );
if( t < b[p] ) ILLEGAL;
else if( t > b[p] ) ans = 2LL*ans%MOD;
}
}
for(int i=1;i<=N;i++) {
if( !vis[i] ) {
int t = 0, p = i;
do {
vis[p] = true, t++, p = a[p];
}while( !vis[p] );
cnt[t]++;
}
}
for(int i=1;i<=N;i++) {
int del = 0;
for(int j=0;j<=cnt[i];j+=2) {
int tmp = 1;
tmp = 1LL*tmp*comb(cnt[i], j)%MOD*f[j]%MOD*pow_mod(i, j/2)%MOD;
if( (i & 1) && i != 1 ) {
tmp = 1LL*tmp*pow_mod(2, cnt[i] - j)%MOD;
}
del = (del + tmp) % MOD;
}
ans = 1LL*ans*del%MOD;
}
printf("%d\n", ans);
}

@details@

其实细节比较多。不过这道题最奇妙的是数形结合利用建图简化思维。

@atcoder - AGC008E@ Next or Nextnext的更多相关文章

  1. [AGC008E] Next or Nextnext [环套树森林+结论讨论]

    题面 传送门 思路 p到a 首先,本题中如果对于所有的$i$,连边$<i,p_i>$,那么可以得到一批环 那么这个题另外一点就是,可以变成连边$<i,p_{p_i}>$ 我们分 ...

  2. AGC008E Next or Nextnext(组合计数,神奇思路)

    神仙题. 排列计数,一种常见的做法是 \(i\) 向 \(p_i\) 连边. 然而这里这个就逼迫我们只能从 \(i\) 向 \(a_i\) 连边. 不过没关系,考虑从 \(i\) 向 \(p_i\) ...

  3. AtCoder刷题记录

    构造题都是神仙题 /kk ARC066C Addition and Subtraction Hard 首先要发现两个性质: 加号右边不会有括号:显然,有括号也可以被删去,答案不变. \(op_i\)和 ...

  4. NOIp2018模拟赛四十

    今天太晚了...题解到时候补吧(flag立好) 成绩:100+0+0=100 感觉A题本质暴力贪心?C题一道水题我居然没做...亏爆 A:[agc011e]increasing numbers B:[ ...

  5. Atcoder Grand Contest 008 E - Next or Nextnext(乱搞+找性质)

    Atcoder 题面传送门 & 洛谷题面传送门 震惊,我竟然能独立切掉 AGC E 难度的思维题! hb:nb tea 一道 感觉此题就是找性质,找性质,再找性质( 首先看到排列有关的问题,我 ...

  6. AGC008E:Next or Nextnext

    传送门 考虑转化成图论问题,\(i\) 向 \(p_i\) 连边,那么合法方案一定是形成了若干个简单环或自环 考虑一个环内的情况: 如果 \(a_i=p_i\),那么 \(i\) 向 \(a_i\) ...

  7. AtCoder Grand Contest 008

    AtCoder Grand Contest 008 A - Simple Calculator 翻译 有一个计算器,上面有一个显示按钮和两个其他的按钮.初始时,计算器上显示的数字是\(x\),现在想把 ...

  8. AtCoder Regular Contest 061

    AtCoder Regular Contest 061 C.Many Formulas 题意 给长度不超过\(10\)且由\(0\)到\(9\)数字组成的串S. 可以在两数字间放\(+\)号. 求所有 ...

  9. AtCoder Grand Contest 001 C Shorten Diameter 树的直径知识

    链接:http://agc001.contest.atcoder.jp/tasks/agc001_c 题解(官方): We use the following well-known fact abou ...

随机推荐

  1. 测试工程中引入Masonry记录

    测试工程中需要引入Masonry,在进行添加新库时发现了几个问题,记录如下,方便有相同问题的朋友查找解决:   1,podfile中添加 pod ‘Masonry’ 后,pod install --v ...

  2. .NET Core HttpClientFactory+Consul实现服务发现

    前言 上篇文章.NET Core HttpClient+Consul实现服务发现提到过,HttpClient存在套接字延迟释放的问题,高并发情况导致端口号被耗尽引起服务器拒绝服务的问题.好在微软意识到 ...

  3. {dede:channelartlist} 改变偶数的class

    {dede:channelartlist} <div {dede:global.itemindex runphp='yes'} if((@me %2) == 0){ @me = 'class=& ...

  4. WordPress安全 - 隐藏保护wp-login.php后台登陆入口

    我们在基本的设置账户用户名和密码安全基础上,最好把这个登录入口限制访问或者隐藏,之前也有看到一些教程说安装插件,比如安装Stealth Login Page插件可以设置登录页面后的参数,与我要设置的非 ...

  5. windows10安全及性能优化

    一.关闭一些服务. Google 更新服务 (gupdate) Google 更新服务 (gupdatem) HomeGroupListener HomeGroupProvider Xbox Live ...

  6. Spring MVC介绍和第一个例子

    1.Spring mvc概述 spring mvc是spring提供给web应用框架设计,实际上MVC框架是一个设计理念.它不仅存在java世界中而且广泛在于各类语言和开发中,包括web的前端应用.对 ...

  7. vue使用IntersectionObserver实现无限下拉信息流

    基于 IntersectionObserver 异步监听方法,实现无线信息流下拉加载, https://developer.mozilla.org/zh-CN/docs/Web/API/Interse ...

  8. Android系统签名简介

    apk的签名,简单说开发者可以通过签名 对应用进行标识和更新.包名在一个设备上是唯一的,这样可以避免被相同包名应用随意覆盖安装.这是一个非常重要的安全功能.系统中的签名文件,也是对系统中应用进行签名, ...

  9. unicode、encode、decode

    1.encode与decode:unicode经过encode -> utf-8,反过来为decode. 爬虫读取网页内容和pandas读取csv时,会把读取到的文字内容转成unicode,当我 ...

  10. Chisel3 - Tutorial - ByteSelector

    https://mp.weixin.qq.com/s/RQg2ca1rwfVHx_QG-IOV-w   字节选择器.   参考链接: https://github.com/ucb-bar/chisel ...