BZOJ 2976: [Poi2002]出圈游戏 HDU 5668 CRT

2976: [Poi2002]出圈游戏

题目连接：

http://www.lydsy.com/JudgeOnline/problem.php?id=2976

Description

Input

中第一行有一个正整数n, 2 <= n <= 20，第二行有n 个整数其中第i个整数表示编号为i 的小朋友第i个出圈。

Output

求最小的K，如果不存在，则输出一个单词“NO”

Sample Input

4

1 4 2 3

Sample Output

5

Hint

题意

题解：

转化一下，其实就是模线性方程组，这个玩意儿用CRT怼一波就好了

然而BZOJ这道题不是输出NO，而是输出NIE

BC这道题却题面写错了，我比赛后才反应过来这道题的题意

sad……

代码是BZOJ的那道题代码，HDU的需要加上test，和 CA是SB这句话

代码

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <cstring>
using namespace std;
const int M = 30;

long long Egcd (long long a, long long b, long long &x, long long &y)
{
    if (b == 0)
    {
        x = 1, y = 0;
        return a;
    }
    long long d, tp;
    d = Egcd (b, a%b, x, y);
    tp = x;
    x = y;
    y = tp - a/b*y;
    return d;
}

long long CRT2 (long long b[], long long n[], int num)
{
    int i;
    bool flag = false;
    long long n1 = n[0], n2, b1 = b[0], b2, bb, d, t, k, x, y;
    for (i = 1; i < num; i++)
    {
        n2 = n[i], b2 = b[i];
        bb = b2 - b1;
        d = Egcd (n1, n2, x, y);
        if (bb % d)     //模线性解k1时发现无解
        {
            flag = true;
            break;
        }
        k = bb / d * x;    //相当于求上面所说的k1【模线性方程】
        t = n2 / d;
        if (t < 0) t = -t;
        k = (k % t + t) % t;    //相当于求上面的K`
        b1 = b1 + n1*k;
        n1 = n1 / d * n2;
    }
    if (flag)
        return -1;           //无解
    return b1;    //形成的解：b1, b1+n1, b1+2n1,..., b1+xni...
}
int vis[30];
long long b[M], n[M];
long long t, num, i, cc = 1;
pair<int,int>P[M];
int main()
{
    scanf ("%d",&num);
    for (i = 0; i < num; i++)
        n[i]=num-i;
    for (i = 0; i < num; i++)
    {
        scanf("%d",&P[i].first);
        P[i].second = i;
    }
    sort(P,P+num);
    int las = num-1;
    for (i = 0; i < num; i++)
    {
        int x;x=P[i].second;
        int ttt = 0;
        while(1)
        {
            las++;
            las%=num;
            if(vis[las])continue;
            if(las==x)
            {
                b[i]=ttt;
                vis[las]=1;
                break;
            }
            ttt++;
        }
    }
    long long p = CRT2(b,n,num);
    if(p==-1)printf("NIE\n");
    else printf("%lld\n",p+1);
    return 0;
}