JZOJ 5421. 【NOIP2017提高A组集训10.25】嘟嘟噜

xiaoxiao2021-02-28  8

Description

由于众所周知的原因, 冈部一直欠真由理一串香蕉. 为了封上真由理的嘴, 冈部承诺只要真由理回答出这个问题, 就给她买一车的香蕉: 一开始有n 个人围成一个圈, 从1 开始顺时针报数, 报出m 的人被机关处决. 然后下一个人再从1 开始报数, 直到只剩下一个人. 红莉栖: “这不就是约瑟夫问题吗…” 伦太郎: “助手你给我闭嘴!” 真由理虽然已经晕头转向了, 但听到有一车的香蕉, 两眼便放出了光芒. “那个呢, 真由氏很想要一车子的香蕉呢. 如果可以帮帮我的话, 我可以把一些香蕉分给你哟, 诶嘿嘿. 拜托你啦.”

Input

第一行一个整数T, 表示数据组数. 接下来T 行, 每行两个整数n;m.

Output

对于每组数据, 输出一行一个整数, 表示幸存者的编号.

Sample Input

5 4 6 2 8 2 9 8 8 7 9

Sample Output

3 1 2 4 7

Data Constraint

Solution

经典的 Josephus 问题, 参见 Wikipedia 。网址:Josephus Problem

设函数 f(n,m) 表示 n,m 的答案,那么就有:

时间复杂度 O(m log n)

Code

#include<cstdio> using namespace std; inline int read() { int X=0,w=1; char ch=0; while(ch<'0' || ch>'9') {if(ch=='-') w=-1;ch=getchar();} while(ch>='0' && ch<='9') X=(X<<3)+(X<<1)+ch-'0',ch=getchar(); return X*w; } inline int f(int x,int y) { if(x==1) return 0; if(x<y) return (f(x-1,y)+y)%x; int x1=x-x/y; return (long long)(f(x1,y)-x%y+x1)%x1*y/(y-1); } int main() { int T=read(); while(T--) { int n=read(),m=read(); printf("%d\n",f(n,m)+1); } return 0; }
转载请注明原文地址: https://www.6miu.com/read-2150032.html

最新回复(0)