苟为蒟蒻又何妨

Gcd 题目背景 SOURCE：NOIP2015-SHY-2 题目描述 给出n个正整数，放入数组 a 里。 问有多少组方案，使得我从 n 个数里取出一个子集，这个子集的 gcd 不为 1 ，然后我再从剩下的数中取出一个数，把他放进刚刚取出的子集里，使得 gcd 为 1 。 输出方案数 mod （10^9 + 7）。 输入格式 第一行一个数 n 。 第二行 n 个数，表示 a 数组...

【问题描述】 根据哥德巴赫猜想（每个不小于 6 的偶数都可以表示为两个奇素数之和），定义 哥德巴赫矩阵 A 如下：对于正整数对(i,j)，若 i+j 为偶数且 i,j 均为奇素数，则 Ai,j = 1， 否则 Ai,j = 0。现在有若干询问(x1,y1,x2,y2)，你需要回答下列式子的值 【输入】 第一行一个整数 m 接下来 m 行，每行四个整数 x1 y1 x2 y2，表示一...

传送门 一道挺有趣的。 第一眼以为每个数都用miller_rabin判一次，但感觉会被卡时间啊。 继续分析发现可以晒出sqrt(r)中的所有素数，然后用类似埃式筛法的方法晒出[l,r][l,r][l,r]之间的素数。 然后就没了。 代码： #include<iostream> #include<cstdio> #include<cstring> ...

传送门 要将所有置换变成一个轮换，显然轮换的周期是所有置换长度的最小公倍数。 于是我们只需要求长度不超过n，且长度最小公倍数为t的不同置换数。 而我们知道，lcm只跟所有素数的最高位有关。 因此lcm=∏iprimeipi∏iprimeipi\prod _i {prime_i}^{p_i} 。 于是我们可以定义状态f[i][j]表示前i个素数凑出的和为j的方案数。 这个可以用类似于背包...