Solution -「CF 840C」On the Bench
\(\mathcal{Description}\)
Link.
给定数列 \(\{a_n\}\),求排列 \(\{p_n\}\) 的个数,使得 \((\forall i\in[1,n))(a_{p_i}a_{p_{i+1}}\not=k^2)\),其中 \(k\in\mathbb N\)。
\(\mathcal{Solution}\)
首先消掉每个数的平方因子,那么限制条件转化为 \(a_{p_i}\not=a_{p_{i+1}}\),我们可以把相等的数放在一个桶里。设桶的大小 \(s_1,s_2,\cdots,s_m\)。
若直接计数,难免需要考虑用过的 \(a\) 的信息。引入容斥,令 \(f_i\) 表示至少有 \(i\) 对数相邻数相等。则答案为:
\]
再用 DP 求 \(f\),令 \(g(i,j)\) 表示把前 \(i\) 个桶分为 \(j\) 块,保证块内元素相同的方案数。转移考虑第 \(i\) 个桶的贡献:
\]
比较显然嘛,其中 \(\binom{s_m-1}{k-1}\) 是隔板法划分第 \(i\) 个桶,\(s_m!\) 表示元素有序,\(k!\) 表示块无序。
最后,找到 \(f\) 和 \(g\) 的关系:、
\]
可以发现 \(g(m,k)\) 实质上就是“至多有 \(k-1\) 个数不相邻”的方案数。
暴力求这两个东西就好,复杂度 \(\mathcal O(n^3)\)。
\(\mathcal{Code}\)
#include <map>
#include <cstdio>
const int MAXN = 300, MOD = 1e9 + 7;
int n, fac[MAXN + 5], ifac[MAXN + 5], f[MAXN + 5], g[MAXN + 5][MAXN + 5];
std::map<int, int> num;
inline void addeq ( int& a, const int b ) { if ( ( a += b ) >= MOD ) a -= MOD; }
inline int qkpow ( int a, int b, const int p = MOD ) {
int ret = 1;
for ( ; b; a = 1ll * a * a % p, b >>= 1 ) ret = 1ll * ret * ( b & 1 ? a : 1 ) % p;
return ret;
}
inline void init ( const int n ) {
fac[0] = 1;
for ( int i = 1; i <= n; ++ i ) fac[i] = 1ll * i * fac[i - 1] % MOD;
ifac[n] = qkpow ( fac[n], MOD - 2 );
for ( int i = n - 1; ~ i; -- i ) ifac[i] = ( i + 1ll ) * ifac[i + 1] % MOD;
}
inline int C ( const int n, const int m ) {
return n < m ? 0 : 1ll * fac[n] * ifac[m] % MOD * ifac[n - m] % MOD;
}
int main () {
scanf ( "%d", &n ), init ( n );
for ( int i = 1, a; i <= n; ++ i ) {
scanf ( "%d", &a );
for ( int j = 2; j * j <= a; ++ j ) for ( ; ! ( a % ( j * j ) ); a /= j * j );
++ num[a];
}
int indx = 0, las = 0;
g[0][0] = 1;
for ( auto p: num ) {
++ indx, las += p.second;
for ( int j = 1; j <= las; ++ j ) {
int& cur = g[indx][j];
for ( int k = 1; k <= j && k <= p.second; ++ k ) {
addeq ( cur, 1ll * g[indx - 1][j - k]
* C ( p.second - 1, k - 1 ) % MOD * fac[p.second] % MOD * ifac[k] % MOD );
}
}
}
for ( int i = 1; i <= n; ++ i ) f[n - i] = 1ll * g[indx][i] * fac[i] % MOD;
int ans = 0;
for ( int i = 0; i < n; ++ i ) addeq ( ans, ( i & 1 ? MOD - 1ll : 1ll ) * f[i] % MOD );
printf ( "%d\n", ans );
return 0;
}
Solution -「CF 840C」On the Bench的更多相关文章
- Solution -「CF 1342E」Placing Rooks
\(\mathcal{Description}\) Link. 在一个 \(n\times n\) 的国际象棋棋盘上摆 \(n\) 个车,求满足: 所有格子都可以被攻击到. 恰好存在 \(k\ ...
- Solution -「CF 1622F」Quadratic Set
\(\mathscr{Description}\) Link. 求 \(S\subseteq\{1,2,\dots,n\}\),使得 \(\prod_{i\in S}i\) 是完全平方数,并最 ...
- Solution -「CF 923F」Public Service
\(\mathscr{Description}\) Link. 给定两棵含 \(n\) 个结点的树 \(T_1=(V_1,E_1),T_2=(V_2,E_2)\),求一个双射 \(\varph ...
- Solution -「CF 923E」Perpetual Subtraction
\(\mathcal{Description}\) Link. 有一个整数 \(x\in[0,n]\),初始时以 \(p_i\) 的概率取值 \(i\).进行 \(m\) 轮变换,每次均匀随机 ...
- Solution -「CF 1586F」Defender of Childhood Dreams
\(\mathcal{Description}\) Link. 定义有向图 \(G=(V,E)\),\(|V|=n\),\(\lang u,v\rang \in E \Leftrightarr ...
- Solution -「CF 1237E」Balanced Binary Search Trees
\(\mathcal{Description}\) Link. 定义棵点权为 \(1\sim n\) 的二叉搜索树 \(T\) 是 好树,当且仅当: 除去最深的所有叶子后,\(T\) 是满的: ...
- Solution -「CF 623E」Transforming Sequence
题目 题意简述 link. 有一个 \(n\) 个元素的集合,你需要进行 \(m\) 次操作.每次操作选择集合的一个非空子集,要求该集合不是已选集合的并的子集.求操作的方案数,对 \(10^9 ...
- Solution -「CF 1023F」Mobile Phone Network
\(\mathcal{Description}\) Link. 有一个 \(n\) 个结点的图,并给定 \(m_1\) 条无向带权黑边,\(m_2\) 条无向无权白边.你需要为每条白边指定边权 ...
- Solution -「CF 599E」Sandy and Nuts
\(\mathcal{Description}\) Link. 指定一棵大小为 \(n\),以 \(1\) 为根的有根树的 \(m\) 对邻接关系与 \(q\) 组 \(\text{LCA}\ ...
随机推荐
- C# - 逆变的具体应用场景
前言 早期在学习泛型的协变与逆变时,网上的文章讲解.例子算是能看懂,但关于逆变的具体应用场景这方面的知识,我并没有深刻的认识. 本文将在具体的场景下,从泛型接口设计的角度出发,逐步探讨逆变的作用,以及 ...
- 【Java常用类】BigDecimal
BigDecimal 一般的Float类和Double类可以用来做科学计算或工程计算,但在商业计算中, 要求数字精度比较高,故用到java.math.BigDecimal类. BigDecimal类支 ...
- Javascript实现让小图片一直跟着鼠标移动
Javascript实现让小图片一直跟着鼠标移动实例 注意:图片可能加载不出来,注意更换 <!doctype html> <html> <head> <met ...
- Solon Web 开发,四、请求上下文
Solon Web 开发 一.开始 二.开发知识准备 三.打包与运行 四.请求上下文 五.数据访问.事务与缓存应用 六.过滤器.处理.拦截器 七.视图模板与Mvc注解 八.校验.及定制与扩展 九.跨域 ...
- How To Remove Systemd Service
Method systemctl stop [servicename] systemctl disable [servicename] rm /your/service/locations/[serv ...
- IoC容器-Bean管理(bean作用域)
IoC操作Bean管理(bean作用域) 1,在Spring里面,设置创建bean实例是单实例还是多实例 2,在Spring里面,默认情况下,bean是单实例对象 3,如何设置单实例还是多实例 (1) ...
- 总是记不住但又总是要用的css
有没有经常遇到一些样式每次写都要用百度呢?我收集了一些我平时经常要用到的但又总是记不住的样式.有错误的地方欢迎指正.转载请注明出处. 一.设置input 的placeholder的字体样式 input ...
- vue开源项目有点全
目录 UI组件 开发框架 实用库 服务端 辅助工具 应用实例 Demo示例 UI组件 element ★31142 - 饿了么出品的Vue2的web UI工具套件 Vux ★14104- 基于Vue和 ...
- ApacheCN 计算机视觉译文集 20211110 更新
OpenCV3 和 Qt5 计算机视觉 零.前言 一.OpenCV 和 Qt 简介 二.创建我们的第一个 Qt 和 OpenCV 项目 三.创建一个全面的 Qt + OpenCV 项目 四.Mat和Q ...
- vue组件中使用<transition></transition>标签过渡动画
直接上代码: <!DOCTYPE html> <html lang="en"> <head> <meta charset="UT ...