P2714-四元组统计【数论,容斥】

正题

题目链接:https://www.luogu.com.cn/problem/P2714


题目大意

给出nnn个数,求有多少个(i,j,k,l)(i,j,k,l)(i,j,k,l)使得gcd(ai,aj,ak,al)=1gcd(a_i,a_j,a_k,a_l)=1gcd(ai,aj,ak,al)=1


解题思路

我们设fif_ifi表示gcdgcdgcd和为iii的方案数。FiF_iFi表示gcdgcdgcd和为iii的倍数的方案数。

先考虑如何统计FiF_iFi,其实就是选出444iii的倍数的方案数,可以统计iii的倍数的个数即可计算。

然后有fi=Fi−∑i∣jfjf_i=F_i-\sum_{i|j}f_jfi=Fiijfj

时间复杂度O(Tnlog⁡n)O(Tn\log n)O(Tnlogn)


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=11000;
ll n,f[N];
int main()
{while(scanf("%lld",&n)!=EOF){memset(f,0,sizeof(f));for(ll i=1;i<=n;i++){ll x;scanf("%lld",&x);f[x]++;}for(ll i=1;i<N;i++)for(ll j=i<<1;j<N;j+=i)f[i]+=f[j];for(ll i=N-1;i>=1;i--){f[i]=f[i]*(f[i]-1)*(f[i]-2)*(f[i]-3)/24;for(ll j=2*i;j<N;j+=i)f[i]-=f[j];}printf("%lld\n",f[1]);}
}

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

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

相关文章

【并查集】Supermarket(poj 1456/luogu-UVA1316)

Supermarket poj 1456 luogu-UVA1316 题目大意&#xff1a; 有一堆物品&#xff0c;每一件物品都有自己的价值和保质期&#xff0c;每天只能卖出一件物品&#xff0c;问最大价值是多少 原题&#xff1a; 题目描述 有一个商店有许多批货&#xff0c;每一批货又有N(0<N…

Jimu : .Net Core 分布式微服务框架介绍

一、前言近些年一直浸淫在 .Net 平台做企业应用开发&#xff0c;用过的 .Net 框架不多&#xff08;具体数量不清&#xff0c;印象深刻的有 Asp.Net MVC&#xff0c;WPF&#xff0c;其他很多都是基于微软开发的框架做些封装而形成新的框架&#xff0c;大都是还没起好名就湮灭在历…

组合数学学习笔记

常见组合计数 n球m盒分配问题 球有别&#xff0c;盒子有别&#xff0c;盒子可空&#xff1a;m^n 每个同学都有m种选择 球无别&#xff0c;盒子有别&#xff0c;盒子不可空&#xff1a;C(n-1,m-1) 隔板法 球无别&#xff0c;盒子有别&#xff0c;盒子可空&#xff1a;C(nm-1,m-1…

P3287-[SCOI2014]方伯伯的玉米田【二维树状数组,dp】

正题 题目链接:https://www.luogu.com.cn/problem/P3287 题目大意 nnn个玉米高度不同&#xff0c;可以选择kkk个区间拔高111个高度&#xff0c;求最长不降子序列长度。 解题思路 显然每次拔高都是拔一个后缀&#xff0c;所以我们设fi,jf_{i,j}fi,j​表示到第iii个玉米&#x…

【背包】买装备

买装备 题目大意&#xff1a; 有n件物品&#xff0c;每件物品有它的物抗&#xff0c;魔抗&#xff0c;价格&#xff0c;现在要在物抗魔抗各不小于一个值的前提下&#xff0c;使价格最小&#xff08;每件物品只能买一件&#xff09; 原题&#xff1a; 题目描述 mxy 沉迷于一…

.NetCore2.1 WebAPI 根据swagger.json自动生成客户端代码

前言上一篇博客中我们可以得知通过Swagger插件可以很方便的提供给接口开发者在线调试&#xff0c;但是实际上Swagger附带的功能还有很多&#xff0c;比如使用NSwag生成客户端调用代码&#xff0c;进一步解放接口开发者。NSwag NSwag是一个发布在GitHub上的开源项目&#xff0c;…

BZOJ1485: [HNOI2009]有趣的数列

题解&#xff1a;考虑按顺序从小到大&#xff0c;加入数字&#xff0c;将加入奇数位置看作入栈&#xff0c;加入偶数位置看作出栈。为什么可以&#xff1f;考虑我们要保证相邻奇数小于偶数&#xff0c;所以必须先填上一个奇数的位置才能填偶数的位置&#xff0c;既时刻保证奇数…

jzoj4223-旅游【并查集】

正题 题目大意 nnn个点mmm条边&#xff0c;qqq次询问走边权小于xxx的能联通的点对数。 解题思路 将边权排序&#xff0c;然后并查集预处理答案即可。 时间复杂度O(mlog⁡m)O(m\log m)O(mlogm) codecodecode #include<cstdio> #include<cstring> #include<alg…

【dfs】买门票

买门票 题目大意&#xff1a; 给出一些字母&#xff0c;求出可组合成的组合&#xff08;要按顺序&#xff0c;和一定的规则&#xff09; 原题&#xff1a; 题目描述 mxy 正要经过新世界的大门。 现在有很多人在门口排队&#xff0c;每个人将会被发到一个有效的通行密码作为…

HDU1812 - Count the Tetris

polya自主ac的第一道&#xff0c;讨论方法&#xff1a;先把奇偶分开(1)顺时针0度&#xff0c;90度&#xff0c;180度&#xff0c;270度 (2)镜像竖线&#xff0c;水平线&#xff0c;两条对角线。分别推出公式计算&#xff0c;实在推不出来&#xff0c;写个模拟暴力找循环节&…

使用.NetCore 控制台演示 熔断 降级(polly)

1、熔断降级的概念&#xff1a; 熔断&#xff1a;我这里有一根长度一米的钢铁&#xff0c;钢铁的熔点1000度&#xff08;假设&#xff09;&#xff0c;现在我想用力把这根钢铁折弯&#xff0c;但是人的力有限达不到折弯的点&#xff0c;然后我使用火给钢铁加热&#xff0c;每隔…

I Got a Matrix!

I Got a Matrix! 题目大意&#xff1a; 给一个矩阵&#xff0c;求出边上的数的和 原题&#xff1a; 题目描述 给定一个 n ∗ m 的矩阵 A&#xff0c;询问位于矩阵边缘的元素之和。所谓矩阵边缘的元素&#xff0c;就是第一行和 最后一行的元素以及第一列和最后一列的元素。…

jzoj6803-NOIP2020.9.26模拟tom【构造】

正题 题目大意 nnn个点的一棵树&#xff0c;给每个点一个权值是1∼a1\sim a1∼a或−1∼−b-1\sim -b−1∼−b。每次选择正负中一个绝对值最小的删去使得无论如何选择都不会将树分成两个联通块。 解题思路 因为可以随意选择&#xff0c;所以aaa和−b-b−b的点一定要连在一起&am…

给正在努力的您几条建议(附开源代码)

前言我是一名在广州的某家互联网公司工作&#xff0c;并有6年工作经验&#xff0c;奔着架构师与微软MVP为目标的老鸟程序员。最近回顾了下多年来走的路&#xff0c;有不少的弯路。今天不说技术&#xff0c;而是总结了一些职业生涯上的建议与大家分享。虽说今天不说技术&#xf…

BZOJ3527: [Zjoi2014]力 [FFT]

化简之后&#xff0c;发现减号左边的式子是一个卷积。右边的式子&#xff0c;把一个函数倒序就是卷积&#xff0c;分别FFT&#xff0c;求解答案。 大佬blog: https://blog.csdn.net/kyleyoung_ymj/article/details/51721495 #include <bits/stdc.h> #define pi acos(-1.0…

jzoj6804-NOIP2020.9.26模拟jerry【dp】

正题 题目大意 nnn个数n−1n-1n−1个或−-−号的一个算式&#xff0c;加上一些括号使得运算结果最大。 解题思路 对于号左右的两个数字&#xff0c;无论如何加括号都不会变结果&#xff0c;所以我们可以合并这些数。然后就只剩下了减号&#xff0c;此时除了第一个就可以随意选…

.Net Core Cors中间件解析

同源策略和资源跨域共享1、同源策略同源策略&#xff0c;它是由Netscape提出的一个著名的安全策略。现在所有支持JavaScript 的浏览器都会使用这个策略。所谓同源是指&#xff0c;域名&#xff0c;协议&#xff0c;端口相同。1.1、目的主要是为了保证用户信息的安全&#xff0c…

BZOJ 3513: [MUTC2013]idiots [FFT]

统计每种长度的木棒数量&#xff0c;先计算出两根棒子能构成的长度&#xff0c;想到卷积。1.拿这个序列卷积自己 2.计算重算的部分&#xff0c;首先是一条边自己和自己的这种情况&#xff0c;另一种是(i,j)和(j,i)这种形式。第一种&#xff0c;可以枚举读入的木棒长度&#xff…

【暴力】I Liked Matrix!

I Liked Matrix! 题目大意&#xff1a; 有一个矩阵&#xff0c;求出某一些子矩阵内最大的数 原题&#xff1a; 题目描述 给定一个 n ∗ m 的矩阵 A&#xff0c;对其进行 q 次询问&#xff1a;以 (x1, y1) 为左上角&#xff0c;(x2, y2) 为右下角的子矩 阵中&#xff0c;所…

jzoj6805-NOIP2020.9.26模拟speike【扫描线】

正题 题目大意 一个平面上&#xff0c;起点是(0,0)(0,0)(0,0)&#xff0c;终点是(0,t)(0,t)(0,t)。有nnn个没有共同面积的矩形障碍物&#xff0c;对于每个障碍物不可以从内部穿过而可以从边上走过。求最短路。 解题思路 因为没有共同面积所以横坐标不会变小。并且只有在边界处…