【codeforces 757E】Bash Plays with Functions
【题目链接】:http://codeforces.com/problemset/problem/757/E
【题意】
给你q个询问;
每个询问包含r和n;
让你输出f[r][n];
这里f[0][n]是n分解成两个数u,v的乘积的个数;
这里u和v互质;
而f[r][n]当r>0时,有个递推式;
【题解】
那个递推式等价于
即n的所有因子x的f[r][x]的和;
这里需要知道;
f[0][n]=2n的不同质因子个数
且容易得到
当a和b互质的时候,f[0][a*b]=f[0][a]*f[0][b];
则f[0]是积性函数;
有个定理就是;
如果是类似上面那个递推
如果且f为积性函数,则g也为积性函数;
则我们在算
f[r][n]的时候可以对n分解质因数;
即f[r][p1^q1*p2^q2…pm^qm];
这样p1^q1..p2^q2..pm^qm都是互质的了;
则可以利用积性函数把它们分解成
f[r][p1^q1]*f[r][p2^q2]..*f[r][pm^qm];
这里
f[0][p^q]的值除了q=0的时候为1之外,其他情况都为2,因为p是质数,可知不同的质因子个数只有1个,也即这个质数p;所以都是2^1;
q为0的时候,只有1,1这种情况,所以只有一对;为1;
(所以这里的p只要为质数,是几都无所谓,f[0]的值都一样)
而
根据
可知
f[r+1][pq]=∑f[r][p0..q]⋅⋅⋅①
因为p0..q都是pq的因子
这样其实就能看出来了,要算f[r][p^q]的时候,根本不需要知道p是多少;
f[r][p^q]只与q和r有关;
根据①式能够很容易写个dp;
具体的
g[0][0] = 1,g[0][1..n] = 2;
g[i][0] = g[i-1][0];
g[i][j]=∑g[i−1][0..j]
对于要算的
f[r][p1q1]∗f[r][p2q2]..∗f[r][pmqm];
直接输出
g[r][q1]∗g[r][q2]∗...∗g[r][qm]
就好;
在对n分解质因数的时候,好像用到了更牛逼的方式.
类似筛法的东西吧.
快速获取某个值的质因子;
(用cin竟然会超时…再也不用了!…还是算了吧,cin还是比较方便的,记住可能会因为这个超时就好…)
【Number Of WA】
5
【完整代码】
#include <bits/stdc++.h>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define ms(x,y) memset(x,y,sizeof x)
#define Open() freopen("F:\\rush.txt","r",stdin)
#define Close() ios::sync_with_stdio(0),cin.tie(0)
typedef pair<int,int> pii;
typedef pair<LL,LL> pll;
const int dx[9] = {0,1,-1,0,0,-1,-1,1,1};
const int dy[9] = {0,0,0,-1,1,-1,1,-1,1};
const double pi = acos(-1.0);
const int N = 1e6;
const int MOD = 1e9+7;
int dp[N+100][21],lp[N+10];
void pre()
{
lp[1] = 1;
rep1(i,2,N)
if (lp[i]==0)
for (int j = i;j<=N;j+=i)
lp[j] = i;
rep1(i,1,20)
dp[0][i] = 2;
rep1(i,0,N)
dp[i][0] = 1;
rep1(i,1,N)
rep1(j,1,20)
dp[i][j] = (dp[i][j-1] + dp[i-1][j])%MOD;
}
int main()
{
//Open();
pre();
int q;
scanf("%d",&q);
while (q--)
{
int r,n;
scanf("%d%d",&r,&n);
LL ans = 1;
while (n>1)
{
int cnt = 0,x = lp[n];
while (n%x==0)
{
cnt++;
n/=x;
}
ans = (ans*dp[r][cnt])%MOD;
}
printf("%lld\n",ans);
}
return 0;
}
【codeforces 757E】Bash Plays with Functions的更多相关文章
- 【codeforces 757B】 Bash's Big Day
time limit per test2 seconds memory limit per test512 megabytes inputstandard input outputstandard o ...
- 【Codeforces 757B】 Bash's big day
[题目链接] 点击打开链接 [算法] 若gcd(s1,s2,s3....sk) > 1, 则说明 : 一定存在一个整数d满足d|s1,d|s2,d|s3....,d|sk 因为我们要使|s|尽可 ...
- Codeforces 757 E Bash Plays with Functions
Discription Bash got tired on his journey to become the greatest Pokemon master. So he decides to ta ...
- 【codeforces 415D】Mashmokh and ACM(普通dp)
[codeforces 415D]Mashmokh and ACM 题意:美丽数列定义:对于数列中的每一个i都满足:arr[i+1]%arr[i]==0 输入n,k(1<=n,k<=200 ...
- 【codeforces 789C】Functions again
[题目链接]:http://codeforces.com/contest/789/problem/C [题意] 看式子. [题解] 考虑最后的答案区间; 如果那个区间是从奇数位置的数字开始的; 那么奇 ...
- 【30.49%】【codeforces 569A】Music
time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard o ...
- 【codeforces 757A】Gotta Catch Em' All!
time limit per test1 second memory limit per test256 megabytes inputstandard input outputstandard ou ...
- 【codeforces 754B】 Ilya and tic-tac-toe game
time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard o ...
- 【codeforces 604D】Moodular Arithmetic
time limit per test1 second memory limit per test256 megabytes inputstandard input outputstandard ou ...
随机推荐
- 使用Android design support library在Eclipse和Android Studio
背景:为了在低版本号下使用Android 5.0的一些新特新 Eclipse篇: 1.将Android Manager中的Android support library升级到最新,我的为22.2.1' ...
- 纪念2014 TI DSP大奖赛
偶然发现TI官网有新闻报道大奖赛,还有沈洁女士给我们颁奖的照片.纪念一下. 第六届TI DSP及嵌入式大奖赛决赛暨颁奖典礼在厦门大学成功举行 Frances Han 2013-2014 TI DSP ...
- UNIX环境高级编程(5):文件I/O(1)
UNIX系统中的大多数文件I/O仅仅须要用到5个函数:open.read.write.lseek以及close.本章说明的函数常常称为"不带缓冲的I/0",术语不带缓冲指的是每一个 ...
- TS3
let [first, ...rest] = [1, 2, 3, 4]; console.log(first); // outputs 1 console.log(rest); // outputs ...
- poj--2631--Roads in the North(树的直径 裸模板)
Roads in the North Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 2389 Accepted: 117 ...
- bfs初学
BFS: ** 当知道初始和目标状态的,用双向BFS: 无权图最好用BFS 不用重复如队** 实现框架: 抄来的(来源:https://www.luogu.org/blog/stephen2333/s ...
- JavaScript学习笔记——对象的创建
对象是JavaScript基本数据类型,在JavaScript中除了Undefined.Null.布尔型(ture.false).字符串和数字之外,其他的都属于对象. 在JavaScript中,一个对 ...
- 制作ubuntu的U盘启动盘
在制作U盘启动盘之前,请各位先格式化你的U盘. 制作U盘启动盘的工具有很多种,我们这里为大家介绍的是用软碟通制作.所有我们需要有这个软件,如果大家没有可以百度“软碟通”,下载安装一个.然后点击打开.在 ...
- hive2.0函数大全
Hive2.0函数大全(中文版) 摘要 Hive内部提供了很多函数给开发者使用,包括数学函数,类型转换函数,条件函数,字符函数,聚合函数,表生成函数等等,这些函数都统称为内置函数. 目录 数学函数 ...
- 欢迎来到Flask的世界
不多说,直接上文档链接:Flask的文档 教程 API 快速上手