POJ 2800 Joseph's Problem
2015-02-02 15:57
399 查看
给n 和 k 求: ∑1<=i<=n(k
mod i).
p = k/i
k mod i = k - p * i
k mod ( i + 1 ) = k - p * ( i + 1 ) = k mod i - p
k mod ( i + 2 ) = k - p * ( i + 2 ) = k mod i - 2 * p
对于连续的 i ,很多p都是一样的 . 相差的部分是一个等差数列 ,
i 的 范围是 从 i 到 min(k/p,n) 如果 p == 0 则 一直延续到最后
[Submit] [Go Back] [Status]
[Discuss]
/* ***********************************************
Author :CKboss
Created Time :2015年02月02日 星期一 10时53分55秒
File Name :POJ2800.cpp
************************************************ */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long int LL;
LL n,k;
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
while(cin>>n>>k)
{
LL ret=0;
for(int i=1;i<=n;i++)
{
LL p = k/i;
if(p==0)
{
ret += (n-i+1)*k;
break;
}
else
{
LL q = k/p;
if(q>n) q = n;
LL md = k%i;
LL num = q-i+1;
LL temp = md*num-(num-1)*num/2*p;
ret += temp;
i=q;
}
}
cout<<ret<<endl;
}
return 0;
}
mod i).
p = k/i
k mod i = k - p * i
k mod ( i + 1 ) = k - p * ( i + 1 ) = k mod i - p
k mod ( i + 2 ) = k - p * ( i + 2 ) = k mod i - 2 * p
对于连续的 i ,很多p都是一样的 . 相差的部分是一个等差数列 ,
i 的 范围是 从 i 到 min(k/p,n) 如果 p == 0 则 一直延续到最后
Joseph's Problem
Joseph likes taking part in programming contests. His favorite problem is, of course, Joseph's problem. It is stated as follows. There are n persons numbered from 0 to n - 1 standing in a circle. The person numberk, counting from the person number 0, is executed. After that the person number k of the remaining persons is executed, counting from the person after the last executed one. The process continues until only one person is left. This person is a survivor. The problem is, given n and k detect the survivor's number in the original circle. Of course, all of you know the way to solve this problem. The solution is very short, all you need is one cycle: r := 0; for i from 1 to n do r := (r + k) mod i; return r; Here "x mod y" is the remainder of the division of x by y, But Joseph is not very smart. He learned the algorithm, but did not learn the reasoning behind it. Thus he has forgotten the details of the algorithm and remembers the solution just approximately. He told his friend Andrew about the problem, but claimed that the solution can be found using the following algorithm: r := 0; for i from 1 to n do r := r + (k mod i); return r; Of course, Andrew pointed out that Joseph was wrong. But calculating the function Joseph described is also very interesting. Given n and k, find ∑1<=i<=n(k mod i). Input The input file contains n and k (1<= n, k <= 109). Output Output the sum requested. Sample Input 5 3 Sample Output 7 Source Northeastern Europe 2005 |
[Discuss]
/* ***********************************************
Author :CKboss
Created Time :2015年02月02日 星期一 10时53分55秒
File Name :POJ2800.cpp
************************************************ */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long int LL;
LL n,k;
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
while(cin>>n>>k)
{
LL ret=0;
for(int i=1;i<=n;i++)
{
LL p = k/i;
if(p==0)
{
ret += (n-i+1)*k;
break;
}
else
{
LL q = k/p;
if(q>n) q = n;
LL md = k%i;
LL num = q-i+1;
LL temp = md*num-(num-1)*num/2*p;
ret += temp;
i=q;
}
}
cout<<ret<<endl;
}
return 0;
}
相关文章推荐
- Poj 2800 Joseph's Problem - 找规律
- poj 2800 Joseph’s Problem(数论)
- POJ 2800 Joseph’s Problem 数论找规律
- POJ 2800 Joseph's Problem 笔记
- (one day one problem)poj 2800 Joseph's Problem (数学)
- UVa 1363 POJ 2800 Joseph's Problem
- POJ 2800 Joseph's Problem(数论)
- POJ 2800 : Joseph\'s Problem (须仔细分析)
- POJ1012 Joseph
- UVa1363 - Joseph's Problem(数论)
- poj 1681 Painter's Problem(开关问题)
- POJ 2480 Longge's problem(积性函数之欧拉函数与gcd)
- POJ 2800 垂直直方图 解题报告
- 【POJ】1012 Joseph
- poj 1681 Painter's Problem(高斯消元法)
- Joseph - POJ 1012 打表
- [POJ 1012] Joseph(约瑟夫)
- POJ 1012 Joseph题解与心得
- Poj 1012 Joseph
- poj-1012 Joseph