HDU - 4944 FSF’s game
2016-02-07 17:21
429 查看
[align=left]Problem Description[/align]
FSF has programmed a game.
In this game, players need to divide a rectangle into several same squares.
The length and width of rectangles are integer, and of course the side length of squares are integer.
After division, players can get some coins.
If players successfully divide a AxB rectangle(length: A, width: B) into KxK squares(side length: K), they can get A*B/ gcd(A/K,B/K) gold coins.
In a level, you can’t get coins twice with same method.
(For example, You can get 6 coins from 2x2(A=2,B=2) rectangle. When K=1, A*B/gcd(A/K,B/K)=2; When K=2, A*B/gcd(A/K,B/K)=4; 2+4=6; )
There are N*(N+1)/2 levels in this game, and every level is an unique rectangle. (1x1 , 2x1, 2x2, 3x1, ..., Nx(N-1), NxN)
FSF has played this game for a long time, and he finally gets all the coins in the game.
Unfortunately ,he uses an UNSIGNED 32-BIT INTEGER variable to count the number of coins.
This variable may overflow.
We want to know what the variable will be.
(In other words, the number of coins mod 2^32)
[align=left]Input[/align]
There are multiply test cases.
The first line contains an integer T(T<=500000), the number of test cases
Each of the next T lines contain an integer N(N<=500000).
[align=left]Output[/align]
Output a single line for each test case.
For each test case, you should output "Case #C: ". first, where C indicates the case number and counts from 1.
Then output the answer, the value of that UNSIGNED 32-BIT INTEGER variable.
[align=left]Sample Input[/align]
[align=left]Sample Output[/align]
题意:给你个n,让你求在n的范围内。是否能将一个矩形分成若干个同样大小为k的正方形,相应有val值,让你统计在n内的全部可能的分数总值
思路:首先我们来试着求解∑i=1nn∗igcd(nk,ik),那么我们能够确定的是假设能够把n∗m的矩形分成大小为k的正方形的话,那么k一定是gcd(n,
i)的因子。那么对于一项来说由于公式能够变形
n∗i∗kgcd(n,i)
-> n∗(ic1+ic2+...)
{k枚举全部的可能},那么cj是n的因子,那么icj就是因子相应的系数,我们再从全部的i来讲。对于因子我们能够计算出全部可能的数,比方因子cj,我们能够得到cj,
2∗cj,
3∗cj,
4∗cj....n,那么相应的系数就是我们须要的icj,累加起来计算是:
num[cj]=(1+2+...+ncj)=(1+ncj)∗(ncj)2
val[n]=∑i=1nnum[i]
ans[n]=ans[n−1]+val[n]
FSF has programmed a game.
In this game, players need to divide a rectangle into several same squares.
The length and width of rectangles are integer, and of course the side length of squares are integer.
After division, players can get some coins.
If players successfully divide a AxB rectangle(length: A, width: B) into KxK squares(side length: K), they can get A*B/ gcd(A/K,B/K) gold coins.
In a level, you can’t get coins twice with same method.
(For example, You can get 6 coins from 2x2(A=2,B=2) rectangle. When K=1, A*B/gcd(A/K,B/K)=2; When K=2, A*B/gcd(A/K,B/K)=4; 2+4=6; )
There are N*(N+1)/2 levels in this game, and every level is an unique rectangle. (1x1 , 2x1, 2x2, 3x1, ..., Nx(N-1), NxN)
FSF has played this game for a long time, and he finally gets all the coins in the game.
Unfortunately ,he uses an UNSIGNED 32-BIT INTEGER variable to count the number of coins.
This variable may overflow.
We want to know what the variable will be.
(In other words, the number of coins mod 2^32)
[align=left]Input[/align]
There are multiply test cases.
The first line contains an integer T(T<=500000), the number of test cases
Each of the next T lines contain an integer N(N<=500000).
[align=left]Output[/align]
Output a single line for each test case.
For each test case, you should output "Case #C: ". first, where C indicates the case number and counts from 1.
Then output the answer, the value of that UNSIGNED 32-BIT INTEGER variable.
[align=left]Sample Input[/align]
3 1 3 100
[align=left]Sample Output[/align]
Case #1: 1 Case #2: 30 Case #3: 15662489 HintIn the second test case, there are six levels(1x1,1x2,1x3,2x2,2x3,3x3) Here is the details for this game: 1x1: 1(K=1); 1x2: 2(K=1); 1x3: 3(K=1); 2x2: 2(K=1), 4(K=2); 2x3: 6(K=1); 3x3: 3(K=1), 9(K=3); 1+2+3+2+4+6+3+9=30
题意:给你个n,让你求在n的范围内。是否能将一个矩形分成若干个同样大小为k的正方形,相应有val值,让你统计在n内的全部可能的分数总值
思路:首先我们来试着求解∑i=1nn∗igcd(nk,ik),那么我们能够确定的是假设能够把n∗m的矩形分成大小为k的正方形的话,那么k一定是gcd(n,
i)的因子。那么对于一项来说由于公式能够变形
n∗i∗kgcd(n,i)
-> n∗(ic1+ic2+...)
{k枚举全部的可能},那么cj是n的因子,那么icj就是因子相应的系数,我们再从全部的i来讲。对于因子我们能够计算出全部可能的数,比方因子cj,我们能够得到cj,
2∗cj,
3∗cj,
4∗cj....n,那么相应的系数就是我们须要的icj,累加起来计算是:
num[cj]=(1+2+...+ncj)=(1+ncj)∗(ncj)2
val[n]=∑i=1nnum[i]
ans[n]=ans[n−1]+val[n]
#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #define ll __int64 using namespace std; const int maxn = 500005; const ll mod = 1ll<<32; ll num[maxn], dp[maxn]; void cal() { for (ll i = 1; i < maxn; i++) for (ll j = i; j < maxn; j += i) num[j] += (j/i+1) * (j/i) / 2; } void init() { memset(num, 0, sizeof(num)); cal(); dp[1] = 1; for (ll i = 2; i < maxn; i++) { dp[i] = dp[i-1] + num[i]*i; dp[i] = dp[i] % mod; } } int main() { init(); int t, n, cas = 1; scanf("%d", &t); while (t--) { scanf("%d", &n); printf("Case #%d: %I64d\n", cas++, dp ); } return 0; }
相关文章推荐
- C语言左旋转字符串与翻转字符串中单词顺序的方法
- android Scroller类详解
- 使用commons-daemon启动、关闭java程序
- Spring注解之:@RestController
- css3 z-index属性
- bzoj3620 似乎在梦中见过的样子
- 春节
- Linux:getrusage的用法
- ActiveMQ 快速入门教程系列 第二章 发布-订阅者模式实现
- Onenote如何快速实现首行缩进的功能。
- ajax优点
- 利用简洁的C语言代码解决跳台阶问题与约瑟夫环问题
- 用努力的程度拼别人聪明的程度。
- 机房重构——DataGridView数据显示,窗体显示
- Learn what you can do with Spring Boot
- bzoj1051【HAOI2006】受欢迎的牛
- 错排问题
- HDU-1014
- 自己动手写一个 WordPress 插件
- What is a JavaBean exactly?