hdu CA Loves GCD(dp)
一道我想骂人的题,差点把我气炸了。
题意:
求一个数的集合中(非多重集,每个数只出现一次)所有子集的gcd的和。结果MOD10^8+7输出。
输入输出不说了,自己看吧,不想写了。
当时我真把它当作数论题来写了,以为可以推导出什么公式然后化简大量重复的操作的。结果最后也没找到。最后题解说是dp,我同学说是暴力,吐血10升。
然后弄出来dp方程之后还是反复的wa,方程明明没啥问题,愣是卡了2个小时找不出错误,心情烦躁的要命,坑爹的室友还各种看视频打游戏,还不带耳机,我自己只好带着耳机大声放音乐,最后连音乐都听不下去了,恶心的想吐。
后来实在无奈了查了下题解,但是没人用dp写,有个用莫比乌斯反演的orz,还有个用暴力的,不过其实有dp的思想在里面。当然这不重要,重要的是我看见了他MOD加的位置挺有意思的,然后猛然想到我的int爆了!因为需要一个小于10^8的数×一个小于1000的数,这个数有可能爆!我叉!特么这不是故意卡int的意思吗?最后把这个改了终于过了……此时距离比赛结束已经5个小时了,我*!
状态转移方程:
dp[i][a[i]] += 1;
dp[i][j] += dp[i-1][j];
dp[i][gcd[j][a[i]]] += dp[i-1][j];
其中gcd[][]是预处理离线出来的,要不然可能会超时。
状态dp[i][j]表示在前n个数的集合中,gcd为j的集合有多少个。
方程表示三种情况:
- 只有a[i]的集合。
- 不存在a[i],只存在前i-1个数中若干数的集合。
- 存在a[i],且存在前i-1个数中若干数的集合。
时间复杂度为O(n*maxn),其中maxn为a[]数组中的最大值。
具体见代码——
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
#define LL long long const int N = ;
const int Mod = ; int a[N];
LL dp[N][N];
int gcd[N][N];
int t, n; int Gcd(int x, int y)
{
if(x < y)
{
int t = x;
x = y;
y = t;
}
while(y != )
{
int t = y;
y = x%y;
x = t;
}
return x;
} void Table()
{
for(int i = ; i < ; i++)
{
for(int j = ; j <= i; j++)
{
gcd[i][j] = gcd[j][i] = Gcd(i, j);
}
}
} int main()
{
//freopen("test.in", "r", stdin);
Table();
scanf("%d", &t);
for(int tm = ; tm <= t; tm++)
{
scanf("%d", &n);
int maxn = ;
for(int i = ; i < n; i++)
{
scanf("%d", &a[i]);
maxn = maxn > a[i] ? maxn : a[i];
}
memset(dp, , sizeof(dp));
dp[][a[]] = ;
for(int i = ; i < n; i++)
{
dp[i][a[i]] += ; //转移方程1
for(int j = ; j <= maxn; j++)
{
dp[i][j] += dp[i-][j]; //转移方程2
dp[i][gcd[j][a[i]]] += dp[i-][j]; //转移方程3
dp[i][j] %= Mod;
dp[i][gcd[j][a[i]]] %= Mod;
}
}
int ans = ;
for(int i = ; i <= maxn; i++)
{
ans += (dp[n-][i]*i)%Mod; //这里小心dp如果是int可能会爆
ans %= Mod;
} printf("%d\n", ans);
}
return ;
}
自己确实挺弱的,还需要努力,但是今天确实非常烦!所有认为这些没什么好烦的,都是因为他没有身临其境的感觉。
hdu CA Loves GCD(dp)的更多相关文章
- HDU 5656 CA Loves GCD dp
CA Loves GCD 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5656 Description CA is a fine comrade w ...
- hdu-5656 CA Loves GCD(dp+数论)
题目链接: CA Loves GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Ot ...
- HDU 5656 CA Loves GCD (数论DP)
CA Loves GCD 题目链接: http://acm.hust.edu.cn/vjudge/contest/123316#problem/B Description CA is a fine c ...
- HDU 5656 ——CA Loves GCD——————【dp】
CA Loves GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)To ...
- hdu 5656 CA Loves GCD(n个任选k个的最大公约数和)
CA Loves GCD Accepts: 64 Submissions: 535 Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 2 ...
- HDU 5656 CA Loves GCD 01背包+gcd
题目链接: hdu:http://acm.hdu.edu.cn/showproblem.php?pid=5656 bc:http://bestcoder.hdu.edu.cn/contests/con ...
- CA Loves GCD (BC#78 1002) (hdu 5656)
CA Loves GCD Accepts: 135 Submissions: 586 Time Limit: 6000/3000 MS (Java/Others) Memory Limit: ...
- 数学(GCD,计数原理)HDU 5656 CA Loves GCD
CA Loves GCD Accepts: 135 Submissions: 586 Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 2621 ...
- hdu 5656 CA Loves GCD
CA Loves GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)To ...
随机推荐
- sql数据库连接字符串在APP.config配置文件内的两种写法
第一种方法 写法: <configuration> <appSettings> <add key="connectionstring" value=& ...
- Spring Mvc 一个接口多个继承; (八)
在 spring 注解实现里,一个接口一般是不能多继承的! 除非在 bean 配置文件里有 针对这个 实现类的配置: <beans:bean id="icService" c ...
- input同名
form表单传值,当,input的 name相同时候,可以同时传值. 当你提交时,你获取的同name数据会被后加载的那个input替代.容易造成数据混乱,但是,可以在某些特殊情况使用.比如,selec ...
- Spring Boot实战系列-----------邮件发送
快速导航 添加Maven依赖 配置文件增加邮箱相关配置 Service.Test项目代码构建 五种邮件发送类型讲解 文本邮件 html邮件 附件邮件 html内嵌图片邮件 模板邮件 问题汇总 添加ma ...
- 第13月第25天 ios11 uitableview reloaddata contentsize
1. [tableView reloadData]; dispatch_after(dispatch_time(DISPATCH_TIME_NOW, (int64_t)(0.1 * NSEC_PER_ ...
- Dream_Spark-----Spark 定制版:005~贯通Spark Streaming流计算框架的运行源码
Spark 定制版:005~贯通Spark Streaming流计算框架的运行源码 本讲内容: a. 在线动态计算分类最热门商品案例回顾与演示 b. 基于案例贯通Spark Streaming的运 ...
- 【API】网络编程模型、多线程
1.网络通信编程 1)网络通信模型基础知识 TCP Server: WSAStartup() socket() bind() linsten() accept() send/recv() closes ...
- 串口硬流控原理验证RTS与CTS
物理连接(交叉连接) 主机的RTS(输出)信号,连接到从机的CTS(输入)信号. 主机的CTS(输入)信号,连接到从机的RTS(输出)信号. 主机发送过程: 主机查询主机的CTS脚信号,此信号连接到从 ...
- vue 兼容IE报错解决方案
IE 页面空白 报错信息 此时页面一片空白 报错原因 Babel 默认只转换新的 JavaScript 语法(如箭头函数),而不转换新的 API ,比如 Iterator.Generator.Set. ...
- eclipse中 EAR Libraries 是什么?
eclipse中 EAR Libraries 是 开发EJB工程所需的库包. 由于新建web工程时,eclipse并不能智能化的判断是否该项目以后会用到ejb, 所以为了全面考虑 就已经帮用户导入了E ...