J. Product of GCDs(莫比乌斯反演)(2021牛客暑期多校训练营2)

Product of GCDs

∏d=1nd∑[gcd⁡(s1d,s2d,…,skd)=1]∏d=1nd∑i=1ndμ(i)Cf[id]k\prod_{d = 1} ^{n} d ^{\sum[\gcd(\frac{s_1}{d}, \frac{s_2}{d}, \dots, \frac{s_k}{d}) = 1]}\\ \prod_{d = 1} ^{n} d ^{\sum\limits_{i = 1} ^{\frac{n}{d}} \mu(i) C_{f[id]} ^{k}}\\ d=1nd[gcd(ds1,ds2,,dsk)=1]d=1ndi=1dnμ(i)Cf[id]k
筛出质数,然后得到phi(mod)phi(mod)phi(mod),即可欧拉降幂写(提供一个思路整体复杂度nlog⁡2nn \log ^ 2 nnlog2n了)。

#pragma GCC optimize(3,"Ofast","inline")
#include <bits/stdc++.h>using namespace std;const int M = 1e7 + 10, N = 8e4 + 10, maxn = 80000;int prime[M], mu[M], cnt;int f[N], n, k;long long C[N][32], mod, phi;bool st[M];inline long long read() {long long x = 0;char c = getchar();while (c < '0' || c > '9') {c = getchar();}while (c >= '0' && c <= '9') {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return x;
}void Init() {mu[1] = 1;for (int i = 2; i < M; i++) {if (!st[i]) {prime[++cnt] = i;mu[i] = -1;}for (int j = 1; j <= cnt && 1ll * i * prime[j] < M; j++) {st[i * prime[j]] = 1;if (i % prime[j] == 0) {break;}mu[i * prime[j]] = -mu[i];}}
}long long Phi(long long n) {long long ans = n;for (int i = 1; i <= cnt && 1ll * prime[i] * prime[i] <= n; i++) {if (n % prime[i] == 0) {ans = ans / prime[i] * (prime[i] - 1);while (n % prime[i] == 0) {n /= prime[i];}}}if (n > 1) {ans = ans / n * (n - 1);}return ans;
}long long mul(long long x, long long y) { return (__int128)x * y % mod; 
}long long quick_pow(long long a, long long n) {long long ans = 1;while (n) {if (n & 1) {ans = mul(ans, a);}a = mul(a, a);n >>= 1;}return ans;
}int main() {// freopen("data05.in", "r", stdin);// freopen("out.txt", "w", stdout);Init();int T;T = read();while (T--) {n = read(), k = read(), mod = read();for (int i = 1; i < N; i++) {f[i] = 0;}for (int i = 1, x; i <= n; i++) {x = read();f[x]++;}for (int i = 1; i <= maxn; i++) {for (int j = 2 * i; j <= maxn; j += i) {f[i] += f[j];}}phi = Phi(mod);C[0][0] = 1;for (int i = 1; i <= maxn; i++) {C[i][0] = 1;for (int j = 1; j <= 30; j++) {C[i][j] = C[i - 1][j - 1] + C[i - 1][j];if (C[i][j] > phi) {C[i][j] -= phi;}}}long long ans = 1;for (int d = 2; d <= maxn; d++) {long long cur = 0;for (int i = 1; i <= maxn / d; i++) {if (mu[i] == -1) {cur -= C[f[i * d]][k];if (cur < 0) {cur += phi;}}else if (mu[i] == 1) {cur += C[f[i * d]][k];if (cur > phi) {cur -= phi;}}}ans = mul(ans, quick_pow(d, cur));}printf("%lld\n", ans);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/313661.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

微软推出Python免费在线教程视频

开源中国曾报道过最近微软针对 Python 初学者&#xff0c;推出了一套免费的教程视频。这套视频名为 Python for Beginners&#xff0c;该在线教学视频由微软高级项目经理 Christopher Harrison、以及微软 AI Gaming 的商业开发经理 Susan Ibach 共同讲解&#xff0c;在课程中加…

徒手实现Spring的IOC

Mark https://www.jianshu.com/p/9fe5a3c25ab6 https://juejin.im/post/5abe75f351882577b45f2336 https://blog.csdn.net/u012373815/article/details/74937913 https://xilidou.com/2018/01/08/spring-ioc/ IOC原理介绍&#xff1a; http://www.importnew.com/14751.htm…

HDU6956-Pass!(2021杭电多校一)(BSGS)

Pass! f(1)0,f(2)n−1,f(t)(n−2)f(t−1)(t−1)f(t−2)f(1) 0, f(2) n - 1, f(t) (n - 2) \times f(t - 1) (t - 1) \times f(t - 2)f(1)0,f(2)n−1,f(t)(n−2)f(t−1)(t−1)f(t−2)&#xff0c;考虑对通项两边同时加一个xf(t−1)x \times f(t - 1)xf(t−1)。 可以得到f(t…

推荐neter常用优秀开源项目系列之一

.net社区有很多优秀的开源项目&#xff0c;我们今天先推荐6个开源项目&#xff1b;1.MassTransitMassTransit 是一个自由、开源、轻量级的消息总线, 用于使用. NET 框架创建分布式应用程序。MassTransit在现有消息传输上提供了一组广泛的功能, 从而使开发人员能够友好地使用基于…

spring四种依赖注入方式

平常的java开发中&#xff0c;程序员在某个类中需要依赖其它类的方法&#xff0c;通常是new一个依赖类再调用类实例的方法&#xff0c;这种开发存在的问题是new的类实例不好统一管理。 spring提出了依赖注入的思想&#xff0c;即依赖不由程序员控制&#xff0c;而是通过spring…

3085 吃遍赴丝码(分治)

3085 吃遍赴丝码 给定一个长度为nnn的数组&#xff0c;求∑i1n∑jin(j−i1)min{ai,…,aj}max{ai,…,aj}\sum\limits_{i 1} ^{n} \sum\limits_{j i} ^{n} (j - i 1) \times min\{a_i, \dots, a_j\} \times max\{a_i, \dots, a_j\}i1∑n​ji∑n​(j−i1)min{ai​,…,aj​}max…

聊聊高并发下库存加减那些事儿——“异步扣减库存”

聊聊高并发下库存加减那些事儿不定期福利发放哦聊聊高并发下库存加减那些事儿背景一般在日常开发中经常会遇到打折促销&#xff0c;秒杀活动&#xff0c;就如拼多多最近的4999抢券买爱疯11促销活动&#xff0c;毕竟谁的钱也不是大风刮来的&#xff0c;有秒杀有促销必定带来大量…

什么是MVC

什么是MVC&#xff1f; MVC全名是Model View Controller&#xff0c;是模型(model)&#xff0d;视图(view)&#xff0d;控制器(controller)的缩写&#xff0c;一种软件设计典范。 用一种业务逻辑、数据、界面显示分离的方法&#xff0c;将业务逻辑聚集到一个部件里面&#xf…

P3242 [HNOI2015] 接水果(整体二分、扫描线)

P3242 [HNOI2015] 接水果 给定一棵树&#xff0c;定义给定了ppp个盘子&#xff0c;每个盘子是树上u,vu, vu,v两点的路径&#xff0c;且盘子有权值&#xff0c;定义水果&#xff0c;水果也是树上u,vu, vu,v两点间的路径。 有qqq个询问&#xff0c;每次给定u,v,ku, v, ku,v,k&a…

SpringMVC 、Struts的区别

先说结论&#xff1a;目前Spring MVC已经一统江山&#xff0c;Struts正在被抛弃 1、Spring MVC原理 Spring 体系&#xff1a; Spring MVC工作流程图 SpringMVC的工作流程描述 1. 用户向服务器发送请求&#xff0c;请求被Spring前端控制Servelt DispatcherServlet捕获&#x…

P3700 [CQOI2017]小Q的表格(反演、分块)

P3700 [CQOI2017]小Q的表格 给定一个大小为nnn \times nnn的表格&#xff0c;初始时i,ji, ji,j位置上填的是f(i,j)ijf(i, j) i \times jf(i,j)ij&#xff0c;有mmm个操作&#xff0c;每次操作给定a,b,x,ka, b, x, ka,b,x,k&#xff0c;把格子a,ba, ba,b上的值改成xxx&#xf…

基于Asp.Net Core MVC和AdminLTE的响应式管理后台之侧边栏处理

说明&#xff1a;.NET Core版本为:2.2 AdminLTE版本为:2.4.18 Bootstrap版本为:3.4.1 font-awesome版本为:4.7.01、新建项目&#xff1a;AdminLteDemo&#xff0c;添加区域Admin 在Areas/Admin/Views文件夹添加文件并分别添加如下代码&#xff1a;_ViewImports.cshtmladdTagHel…

Spring boot——起步依赖

一、起步依赖 1、是什么 本质上是一个Maven项目对象模型&#xff08;Project Object Model&#xff0c; POM&#xff09;&#xff0c; 定义了对其他库的传递依赖&#xff0c;这些东西加在一起即支持某项功能。 比如&#xff1a; spring-boot-dependencies^ spring-boot-pare…

[翻译] ASP.NET Core 3.0 的新增功能

全文翻译自微软官方文档英文版 Whats new in ASP.NET Core 3.0本文重点介绍了 ASP.NET Core 3.0 中最重要的更改&#xff0c;并提供相关文档的连接。BlazorBlazor 是 ASP.NET Core 中的一个新的框架&#xff0c;用于使用 .NET 构建交互式的客户端 Web UI&#xff1a;使用 C# 而…

Convolution(2021牛客暑期多校训练营4)

Convolution 定义a⊕babgcd⁡(a,b)2a \oplus b \frac{a \times b}{\gcd(a, b) ^ 2}a⊕bgcd(a,b)2ab​&#xff0c;bm∑i1n∑j1naijc[i⊕jm]b_m \sum\limits_{i 1} ^{n} \sum\limits_{j 1} ^{n}a_i \times j ^ c [i \oplus j m]bm​i1∑n​j1∑n​ai​jc[i⊕jm]&#xff0…

Sprig boot自动配置

1、概述 Spring Boot是Spring旗下众多的子项目之一&#xff0c;其理念是约定优于配置&#xff0c;它通过实现了自动配置&#xff08;大多数用户平时习惯设置的配置作为默认配置&#xff09;的功能来为用户快速构建出标准化的应用。Spring Boot的特点可以概述为如下几点&#x…

E. Mocha and Stars(莫比乌斯反演、简单dp)

E. Mocha and Stars ∑a1l1r1∑a2l2r2⋯∑anlnrn[a1a2⋯an≤m][gcd⁡(a1,a2,…,an)1]\sum_{a_1 l_1} ^{r_1} \sum_{a_2 l_2} ^{r_2} \dots \sum_{a_n l_n} ^{r_n} [a_1 a_2 \dots a_n \le m][\gcd(a_1, a_2, \dots, a_n) 1]\\ a1​l1​∑r1​​a2​l2​∑r2​​⋯an​ln…

.NET Core 3.0稳定版发布

一年一度的 .NET 开发者盛会 .NET Conf 2019 已正式开始了&#xff08;举办时间为 9.23-9.25&#xff09;。大会第一天最重磅也是最激动人心的消息莫过于 .NET Core 3.0 稳定版的发布。.NET 项目管理总监 Scott Hunter 在大会宣布了这则消息&#xff0c;并在主题演讲环节围绕 .…

深入理解 Spring Boot Starters 原理(手写Spring boot Start)

一、Spring Boot Starter诞生原因 Spring Boot Starter是在SpringBoot组件中被提出来的一种概念&#xff0c;stackoverflow上面已经有人概括了这个starter是什么东西&#xff0c;想看完整的回答戳这里 Starter POMs are a set of convenient dependency descriptors that you …

FFT字符串匹配(解决通配符问题)

FFT字符串匹配 定义字符串下标从000&#xff0c;开始&#xff0c;有文本串AAA长度为nnn&#xff0c;模式串BBB长度为mmm&#xff0c;我们可以考虑一个函数f(x,y)A(x)−B(y)f(x, y) A(x) - B(y)f(x,y)A(x)−B(y)。 我们设F(x)(x≥m−1)∑i0m−1f(x−m1i,i)F(x)(x \ge m - 1) …