hdu 4944 FSF’s game 数学(公因子)+递推公式
来源:互联网 发布:做网站必备软件 编辑:程序博客网 时间:2024/04/26 22:51
题意:
给定一个整数n,求∑fun(i,j)(1<=i<=j<=n)。其中fun(i,j)=∑i*j/gcd(i/k,j/k)(k为i和j的公因子)。
题解:
令答案ans=dp[n],那么dp[n]=dp[n-1]+∑fun(i,n)(1<=i<=n),dp[1]=1。所以我们只要求出所有的∑fun(i,n)(1<=i<=n),就能求得多有的dp[i]。也就能在O(1)时间内输出结果了。
∑fun(i,n)(1<=i<=n)怎么求?我们先讨论fun(i,n)=n*(i/c1+i/c2+...),其中cj为n和i的最大公约数d的所有因子。
那么我们可以枚举所有n的因子j,对于因子j,存在这个因子且小于等于n的数字中有j,2*j,3*j,...,n。我们可以得到s[j]=(1+2+3+...n/j)*n=(n/j+1)*(n/j)/2*n。那么val[n]=∑fun(i,n)(1<=i<=n)=∑s[j](j为n的因子)。
之后用dp[n]=dp[n-1]+val[n]递推公式就能求得结果了。
#include <cstdio>#include <cmath>#include <cstdlib>#include <cstring>#include <vector>#include <map>#include <iostream>#include <algorithm>using namespace std;#define LL __int64const int maxn=5e5+20;LL mod;LL val[maxn],dp[maxn];void fun(int n){ for(LL i=1;i<n;i++)//枚举因子 { for(LL j=i;j<n;j+=i)//枚举含有因子的数 { val[j]+=(j/i+1)*(j/i)/2; } }}void init(){ mod=1; mod=(mod<<32); memset(val,0,sizeof(val)); fun(maxn); dp[1]=1; for(LL i=2;i<maxn;i++)//递推公式 { dp[i]=dp[i-1]+val[i]*i; dp[i]=dp[i]%mod; }}int main(){ init(); int n,tt=0,T; scanf("%d",&T); while(T--) { scanf("%d",&n); printf("Case #%d: %I64d\n",++tt,dp[n]); } return 0;}
阅读全文
0 0
- hdu 4944 FSF’s game 数学(公因子)+递推公式
- hdu 4944 FSF’s game 数学(公因子)+递推公式
- HDU 4944 FSF’s game 解题报告(递推)
- HDU 4944 FSF’s game(数论+递推)
- HDU4944 FSF’s game【因子和】【递推】
- HDU4944-FSF’s game(递推)
- HDU - 4944 FSF’s game
- hdu 4944 FSF’s game
- hdoj 4944 FSF’s game 【数学思维】
- hdu 4944 FSF’s game(数论)
- HDU 4944 FSF’s game 规律题
- 多校hdu FSF’s game 4944
- HDU 4944 FSF’s game(计数游戏)
- HDOJ 4944 FSF’s game
- hdu4944 FSF’s game
- hdu FSF’s game 欧拉定理,素数筛法
- SOJ-2857(数学递推公式)
- HDU 1297 Children’s Queue 递推 错排公式
- python显示图片
- 学习记录(8)
- POJ 1088 滑雪——DP
- python 公有属性和私有属性
- Codeforces Round #428 (Div. 2) D
- hdu 4944 FSF’s game 数学(公因子)+递推公式
- test ATL OCX by vs2015
- 面试运维
- iptables: No config file. [WARNING]!
- nodejs-socket网络编程
- 关于在Aptana Studio中编写js效果在浏览器中的刷新问题,光标变黑色矩形问题
- 利用Frida从TeamViewer内存中提取密码
- python 访问属性和属性装饰器
- python运算符