ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Luogu2421 [NOI2002]荒岛野人(拓欧)题解

2020-11-16 18:03:30  阅读:235  来源: 互联网

标签:return Luogu2421 题解 野人 int 枚举 Peo NOI2002 拓欧


算法

暴力枚举+拓欧检查

思路

两个野人\(i\)与\(j\)会相遇,当且仅当方程\(c_i + xp_i \equiv c_j + xp_j(mod M)\)有正整数解且该解\(\leq min(l_i,l_j)\)。

又发现题目中保证 \(M \leq 1e6\),于是就可以愉快的暴力枚举了。

而上面的方程可以化为\(x(p_i - p_j) \equiv c_j - c_i(mod M)\),即求不定方程\((p_i - p_j)x + My = c_j-c_i\)的正整数解,可以使用拓展欧几里得定理。

复杂度:\(O(Mn^2logC_i)\)。

Tips

  • \(M\)是不单调的,所以不能二分!
  • \(M\)不能从\(1\)开始枚举,而要从\(max(c_i)\)开始,毕竟每个野人都得有地方住
  • 要是不会拓欧的,可以康康这篇博客

参考代码

#include <cstdio>
#include <algorithm>

using namespace std;

int n,ans;
struct Peo{
    int c,p,l;
}p[20];

bool cmp(Peo x, Peo y){return x.p > y.p;}

int ex_gcd(int a, int b, int &x, int &y){
    if(b == 0){
        x = 1, y = 0;
        return a;
    }
    int tmp = ex_gcd(b, a % b, x, y);
    int ret = x - (a / b) * y;
    x = y, y = ret;
    return tmp;
}

int main(){
    int Max = 0;
    scanf("%d", &n);
    for(int i = 1; i <= n; ++ i)
        scanf("%d%d%d", &p[i].c, &p[i].p, &p[i].l), Max = max(Max, p[i].c);
    sort(p + 1, p + 1 + n, cmp);
    for(ans = Max; ans <= ((int)1e6 + 10); ++ ans){
        int flag = 1;
        for(int i = 1; i < n; ++ i){
            for(int j = i + 1; j <= n; ++ j){
                int a = p[i].p - p[j].p, b = ans, c = p[j].c - p[i].c, x = 0, y = 0;
                int val = ex_gcd(a, b, x, y);
                if(c % val != 0) continue;
                int minn = (c / val * x % (b / val) + b / val) % (b / val);
                if(minn <= min(p[i].l, p[j].l)){
                    flag = 0;
                    break;
                } 
            }
            if(!flag) break;
        }
        if(flag) return printf("%d\n", ans), 0;
    }
    return 0;
}

标签:return,Luogu2421,题解,野人,int,枚举,Peo,NOI2002,拓欧
来源: https://www.cnblogs.com/whenc/p/13986309.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有