模拟 nbut1225 NEW RDSP MODE I

传送门：点击打开链接

题意：输入n，m，x，刚开始有一个1~n的排列，然后定义了一种操作，是将数组中的偶数位数字选出来，按照顺序放到数组最前面，奇数位按照顺序放到偶数位的后面，进行m次这样的操作，输出之后前x个数字

思路：找到循环节T，利用T去约m，然后再将很小的m拿去模拟，输出前x个

一开始就想到找循环节，，刚开始只想到去用找规律的方法去找到通项公式，但是找了好久就是没找到，虽然感觉理论上肯定是有的T^T

但是找规律的时候发现了很多特点：T一定小于等于n，还有就是最刚开始的时候数字1是在第一个位置，当数字1再次出现在第一个位置的时候，刚好就是一个循环节！

所以，我们只需要模拟1的位置，一直模拟到1出现在第一个位置时，循环节就算出来了。复杂度O(n)，1的位置还是很好模拟的，因为只研究了一个数字而已，还是很好找到递推式子的。找到循环节T之后，我们令m=m%T.这样m就变成<=n的了，然后就可以再次模拟

接下来，我们对前x个数字，分别倒着模拟m次，因为现在的m<=n,所以复杂度O(xn)，倒着模拟的公式也是很好找的

最后看别人的代码才发现，，其实模拟的时候，就是一个快速幂(fuck)

总之，还是有些感慨，有时候不一定模拟就非要找到通项公式，找到办法能在较短的时间内算出通项公式，这样也并不算差~

#include<map>#include<set>#include<cmath>#include<stack>#include<queue>#include<cstdio>#include<string>#include<vector>#include<cstring>#include<iostream>#include<algorithm>#include<functional>#define FIN freopen("input.txt","r",stdin)#define FOUT freopen("output.txt","w+",stdout)using namespace std;typedef long long LL;typedef pair<int, int> PII;const int MX = 1e5 + 5;const int INF = 0x3f3f3f3f;int find_t(int n) {    int p = 1, ret = 0;    while(true) {        ret++;        if(p % 2) p = n / 2 + (p + 1) / 2;        else p = p / 2;        if(p == 1) return ret;    }}int solve(int p, int n, int m) {    for(int i = 1; i <= m; i++) {        if(p * 2 <= n) p = p * 2;        else p = (p - n / 2 - 1) * 2 + 1;    }    return p;}int main() {    int n, m, x;//FIN;    while(~scanf("%d%d%d", &n, &m, &x)) {        int t = find_t(n);        for(int i = 1; i <= x; i++) {            printf("%d%c", solve(i, n, m % t), i == x ? '\n' : ' ');        }    }    return 0;}