HDU 5833 Zhu and 772002
HDU 5833 Zhu and 772002
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Description |
题目描述 |
Zhu and 772002 are both good at math. One day, Zhu wants to test the ability of 772002, so he asks 772002 to solve a math problem. But 772002 has a appointment with his girl friend. So 772002 gives this problem to you. There are n numbers a1,a2,...,an. The value of the prime factors of each number does not exceed 2000, you can choose at least one number and multiply them, then you can get a number b . How many different ways of choices can make b is a perfect square number. The answer maybe too large, so you should output the answer modulo by 1000000007. |
Zhu和772002都擅长数学。 某天,Zhu想试试772002的能耐,就给772002出了道数学题。 但772002与女票有约在先,果断甩锅予你。 有n个数a1,a2,...,an。每个数的质因数不超过2000,你可以选择至少一个数并将所选的数相乘,得到b。 有多少种选择方式可以使b为完全平方数。结果可能非常大,输出时模1000000007。 |
Input |
输入 |
First line is a positive integer T, represents there are T test cases. For each test case: First line includes a number n(1≤n≤300),next line there are n numbers a1,a2,...,an,(1≤ai≤1018). |
第一行是一个整数T(T≤100),表示有T个测试用例的数量。 对于每个测试用例: 第一行有一个整数n(1≤n≤300),下一行有n个数a1,a2,...,an,(1≤ai≤1018)。 |
Output |
输出 |
For the i-th test case, first output Case #i: in a single line. Then output the answer of i-th test case modulo by 1000000007. |
对于第i个测试用例,先输出一行Case #i:。 然后输出第i个测试用例的答案模1000000007后的结果。 |
Sample Input - 输入样例 |
Sample Output - 输出样例 |
2 |
Case #1: |
【题解】
分解质因数 + 高斯消元
对于容易一个输入的数进行分解质因数,保存到矩阵的各个行中。
每个数都能看作若干个质数相乘,若每个质数的指数均为偶数,这个数即是完全平方数。因此用10表示奇偶,每个数的合并/化简就能用异或来操作。
然后用异或高斯消元,得到矩阵的秩r,接着得到(n - r)个可自由组合的完全平方数。最后组合数求和,去掉什么都不选的情况,结果为2n - 1。
(啪!迷之巴掌)
好吧,其实我一开始根本不知道高斯消元是什么鬼。刚刚开始的理解就是分解质因数(形成多项式),压入矩阵方便求基底(当成向量来看),然后默默地发现高斯消元就是矩阵化简吧……
【代码 C++】
#include <cstdio>
#include <cstring>
#define pMX 2005
#define mod 1000000007
int prim[] = { }, iP, mtx[][];
void getPrim() {
int i, j;
bool mark[pMX];
memset(mark, , sizeof(mark));
for (i = ; i < pMX; i += ) {
if (mark[i]) continue;
prim[++iP] = i;
for (j = i << ; j < pMX; j += i) mark[j] = ;
}
}
void tPrim(int y, __int64 a) {
int i, j;
for (i = ; i <= iP && a >= prim[i]; ++i) {
if (a / prim[i] * prim[i] == a) {
for (j = ; a / prim[i] * prim[i] == a; a /= prim[i]) ++j;
mtx[y][i] = j & ;
}
}
}
int mTS(int n){
bool inUS[]; memset(inUS, , sizeof(inUS));
int size = , i, j, k, ik;
for (j = ; j <= iP; ++j){
for (i = ; i < n; ++i){
if (inUS[i] == && mtx[i][j] == ) break;
}
if (i == n) continue;
inUS[i] = ; ++size;
for (k = ; k < n; ++k){
if (inUS[k] || mtx[k][j] == ) continue;
for (ik = j; ik <= iP; ++ik) mtx[k][ik] ^= mtx[i][ik];
}
}
return n - size;
}
__int64 qMod(__int64 a, int n){
__int64 opt = ;
while (n){
if (n & ) opt = (opt*a) % mod;
n >>= ;
a = (a*a) % mod;
}
return opt;
}
int main() {
getPrim();
int t, iT, i, n;
__int64 ai;
scanf("%d", &t);
for (iT = ; iT <= t; ++iT) {
printf("Case #%d:\n", iT);
memset(mtx, , sizeof(mtx));
scanf("%d", &n);
for (i = ; i < n; ++i){
scanf("%I64d", &ai);
tPrim(i, ai);
}
printf("%I64d\n", qMod(, mTS(n)) - );
}
return ;
}
HDU 5833 Zhu and 772002的更多相关文章
- hdu 5833 Zhu and 772002 ccpc网络赛 高斯消元法
传送门:hdu 5833 Zhu and 772002 题意:给n个数,每个数的素数因子不大于2000,让你从其中选则大于等于1个数相乘之后的结果为完全平方数 思路: 小于等于2000的素数一共也只有 ...
- HDU 5833 Zhu and 772002 (高斯消元)
Zhu and 772002 题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5833 Description Zhu and 772002 are b ...
- hdu 5833 Zhu and 772002 高斯消元
Zhu and 772002 Problem Description Zhu and 772002 are both good at math. One day, Zhu wants to test ...
- HDU 5833 Zhu and 772002(高斯消元)
题意:给n个数,从n个数中抽取x(x>=1)个数,这x个数相乘为完全平方数,求一共有多少种取法,结果模1000000007. 思路:每个数可以拆成素数相乘的形式,例如: x1 2=2^1 * 3 ...
- HDU 5833 Zhu and 772002 (数论+高斯消元)
题目链接 题意:给定n个数,这n个数的素因子值不超过2000,从中取任意个数使其乘积为完全平方数,问有多少种取法. 题解:开始用素筛枚举写了半天TLE了,后来队友说高斯消元才想起来,果断用模板.赛后又 ...
- hdu 5833 Zhu and 772002 异或方程组高斯消元
ccpc网赛卡住的一道题 蓝书上的原题 但是当时没看过蓝书 今天又找出来看看 其实也不是特别懂 但比以前是了解了一点了 主要还是要想到构造异或方程组 异或方程组的消元只需要xor就好搞了 数学真的是硬 ...
- HDU - 5833: Zhu and 772002 (高斯消元-自由元)
pro:给定N个数Xi(Xi<1e18),保证每个数的素因子小于2e3:问有多少种方案,选处一些数,使得数的乘积是完全平方数.求答案%1e9+7: N<300; sol:小于2e3的素数只 ...
- HDU 5833 Zhu and 772002 ——线性基
[题目分析] 这题貌似在UVA上做过,高精度高斯消元. 练习赛T2,然后突然脑洞出来一个用Bitset的方法. 发现代码只需要30多行就A掉了 Bitset大法好 [代码] #include < ...
- 【HDU 5833】Zhu and 772002(异或方程组高斯消元)
300个最大质因数小于2000的数,选若干个它们的乘积为完全平方数有多少种方案. 合法方案的每个数的质因数的个数的奇偶值异或起来为0. 比如12=2^2*3,对应的奇偶值为01(2的个数是偶数为0,3 ...
随机推荐
- mysql相关总结
mysql设置初始密码和更改密码(ZIP文件解压安装): http://blog.csdn.net/stypace/article/details/38232393
- Linux修改系统以及pip更新源
Linux修改系统以及pip更新源 时间:2015-08-01来源:csdn 作者:henulwj 修改系统更新源 你是否跟我一样在刚看时接触Linux时被系统更新源问题搞得晕头转向,不同的Linux ...
- 在Win8下无法打开 hlp 帮助文件的问题
需要安装Win8针对该问题的补丁程序,并且修改注册表,详细的解决方案: http://support.microsoft.com/kb/917607/zh-cn#fixit4me
- python sklearn环境配置
os:win10 python2.7 主要参照 1.现下载pip.exe,因为很多安装文件都变成whl格式了,这里要注意下载对应python版本的,要用管理员权限,可以参照https://pypi ...
- Python代码项目目录规范v1.0
程序目录规范:bin # 存放可执行程序 xxxx.py # 程序主程序(入口文件)config # 存放配置信息 settings.py # 全局配置文件(可能暂时未应用)db # 存放数据文件 c ...
- linux命令:du 命令
Linux du命令也是查看使用空间的,但是与df命令不同的是Linux du命令是对文件和目录磁盘使用的空间的查看,还是和df命令有一些区别的. 1.命令格式: du [选项][文件] 2.命令功能 ...
- setw和setfill控制输出间隔
在C++中,setw(int n)用来控制输出间隔.例如:cout<<'s'<<setw(8)<<'a'<<endl;则在屏幕显示s a //s与a之间 ...
- 每日一九度之 题目1033:继续xxx定律
时间限制:1 秒 内存限制:32 兆 特殊判题:否 提交:5502 解决:1351 题目描述: 当n为3时,我们在验证xxx定律的过程中会得到一个序列,3,5,8,4,2,1,将3称为关键数, ...
- YTU 2296: KMP模式匹配 二(串)
2296: KMP模式匹配 二(串) 时间限制: 1 Sec 内存限制: 128 MB 提交: 29 解决: 17 题目描述 输入一个主串和一个子串,用KMP进行匹配,问进行几趟匹配才成功,若没成 ...
- I love sneakers!(分组背包HDU3033)
I love sneakers! Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) ...