[agc008E]Next or Nextnext-[dp+思考题]
Description
Solution
然后我谈下个人理解。由于我们的两个条件只要任意满足,则在p的图中i有两种连边法:i->p[i],i->p[p[i]]。
我们考虑在a的图中i->a[i]。可得我们要把p图塞到a图里。
具体分析看题解吧,题解图画的很清晰呀。然后。。就各种dp+乱搞了。
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
const int mod=1e9+;
int n,x;
int a[];
int in[],col[],cir[];
int len[],cnt[];
ll dp[],ans=;
int main()
{
scanf("%d",&n);
for (int i=;i<=n;i++) scanf("%d",&a[i]);
for (int i=;i<=n;i++)
{
in[a[i]]++;
x=i;while (!col[x]) {col[x]=i;x=a[x];}
if (col[x]!=i) continue;
while (!cir[x]) {cir[x]=i;x=a[x];}
}
for (int i=;i<=n;i++)
if ((cir[i]&&in[i]>)||(!cir[i]&&in[i]>)) return printf(""),; int js;
for (int i=;i<=n;i++)
{
if (in[i]) continue;
x=i;js=;
while (!cir[x]) {js++;x=a[x];}
len[x]=js;
}
int _len,fir,st,id;
for (int i=;i<=n;i++)
{
if (!cir[i]) continue;
x=i;
_len=fir=st=id=;
for(;cir[x];x=a[x])
{
id++;cir[x]=;
if(len[x])
{
if (!fir){fir=st=id;_len=len[x];}
else
{
ans=ans*((len[x]<id-st)+(len[x]<=id-st))%mod;
if (!ans) return printf(""),;st=id;
}
}
}
if (fir)
{
ans=ans*((_len<id+fir-st)+(_len<=id+fir-st))%mod;
if (!ans)return printf(""),;
}
else cnt[id]++;
}
for (int i=;i<=n;i++)
{
if (!cnt[i]) continue;
dp[]=;
if (i>&&i%)
for (int j=;j<=cnt[i];j++)
{
dp[j]=dp[j-]*%mod;
if (j>) dp[j]=(dp[j]+dp[j-]*(j-)*i%mod)%mod;
} else
for (int j=;j<=cnt[i];j++)
{
dp[j]=dp[j-];
if (j>) dp[j]=(dp[j]+dp[j-]*(j-)*i%mod)%mod;
}
ans=ans*dp[cnt[i]]%mod;
if (!ans) return printf(""),;
}
cout<<ans;
}
[agc008E]Next or Nextnext-[dp+思考题]的更多相关文章
- [AGC008E] Next or Nextnext [环套树森林+结论讨论]
题面 传送门 思路 p到a 首先,本题中如果对于所有的$i$,连边$<i,p_i>$,那么可以得到一批环 那么这个题另外一点就是,可以变成连边$<i,p_{p_i}>$ 我们分 ...
- AGC008E Next or Nextnext(组合计数,神奇思路)
神仙题. 排列计数,一种常见的做法是 \(i\) 向 \(p_i\) 连边. 然而这里这个就逼迫我们只能从 \(i\) 向 \(a_i\) 连边. 不过没关系,考虑从 \(i\) 向 \(p_i\) ...
- @atcoder - AGC008E@ Next or Nextnext
目录 @description@ @solution@ @accepted code@ @details@ @description@ 给定一个长度为 N 的序列 a,问有多少排列 p,满足对于每一个 ...
- AtCoder刷题记录
构造题都是神仙题 /kk ARC066C Addition and Subtraction Hard 首先要发现两个性质: 加号右边不会有括号:显然,有括号也可以被删去,答案不变. \(op_i\)和 ...
- NOIp2018模拟赛四十
今天太晚了...题解到时候补吧(flag立好) 成绩:100+0+0=100 感觉A题本质暴力贪心?C题一道水题我居然没做...亏爆 A:[agc011e]increasing numbers B:[ ...
- AGC008E:Next or Nextnext
传送门 考虑转化成图论问题,\(i\) 向 \(p_i\) 连边,那么合法方案一定是形成了若干个简单环或自环 考虑一个环内的情况: 如果 \(a_i=p_i\),那么 \(i\) 向 \(a_i\) ...
- [agc011F]Train Service Planning-[线段树优化dp+神秘思考题]
Description 传送门 Solution 请围观lhx大佬的博客(大佬写的太好了我都没有写的动力了em) Code #include<iostream> #include<c ...
- hdu4939 Stupid Tower Defense (DP)
2014多校7 第二水的题 4939 Stupid Tower Defense Time Limit: 12000/6000 MS (Java/Others) Memory Limit: 131 ...
- POJ2677 Tour(DP+双调欧几里得旅行商问题)
Tour Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 3929 Accepted: 1761 Description ...
随机推荐
- Aheadof Time Compilation(AOT) vs (JIT)Just In Time compilation approach
像java这样的高级语言,往往先做好部分编译,在运行时,再使用JIT将前面编译的中间件输出编译为机器语言,放到机器上运行.这可能会影响到运行的性能. 再比如,像angular的web 应用,angul ...
- Django之modelform简介
在django中内置了form类和model类,当页面中的form值和model字段值完全一样时,此时可以通过model生成一个完全一样的form,Django中的modelForm就因此而生. 目标 ...
- "docker build" requires exactly 1 argument(s).
Docker 是怎么样的东西,这里就不说了,这里说说dockerfile创建容器时遇到的问题. 首先我想达到的目的很简单,就是用dockerfile去创建容器,步骤如下: 创建并编辑dockerfil ...
- 转:Eclipse+webservice开发实例
原文地址:http://blog.csdn.net/xw13106209/article/details/7049614 1.参考文献: 1.利用Java编写简单的WebService实例 http ...
- Linode VPS主机套餐方案降低处理方法且不影响数据
使用Linode VPS主机产品经历过512MB内存升级至1GB内存,再升级至2GB内存,以及目前推出1GB内存方案月付10美元.比如我们在使用Linode 2GB内存方案的时候占用资源不是太多,其实 ...
- 解决Oracle11g密码180天过期,账号锁住的问题
1.查看用户的proifle是哪个,一般是default: sql>SELECT username,PROFILE FROM dba_users; 2.查看指定概要文件(如default)的密码 ...
- HTML-head头部浅析
HTML结构 在sublime或HBuildr新建HTML文件,输入html:5,按下tab键后,自动生成的代码大致如下: <!DOCTYPE html> <html lang=&q ...
- springmvc常用的组件,注解,跳转
路径映射 XXXHandlerMapping 随开发配置越来越多 注解到java代码中来简化xml配置 请求到哪个Controller 控制器bean Controller 随着开发配置越来越多 注解 ...
- Mycat问题总结
Mycat问题总结 一丶自增主键设置 Mycat提供了几种设置自增主键的方式 本地文件方式 数据库方式 服务器时间戳方式 分布式ZK-ID生成器 第一种和第二种只适合单点设置,对于集群不适用.第四种方 ...
- 跟我一起阅读Java源代码之HashMap(二)
上一节中实现的SimpleHashMap,没有解决冲突的问题,这一节我们继续深入 由于table的大小是有限的,而key的集合范围是无限大的,所以寄希望于hashcode散落,肯定会出现多个key散落 ...