hduoj 5223 GCD 构造检查

hduoj 5223 GCD 构造检查

有一个数组$A_1 \cdots A_N$，数组中每一个元素是在区间$[1,10^9]$之间。现在有Q个问题，每一个问题是GCD($A_{L_i},~A_{L_i + 1},~A_{L_i + 2},~...,~A_{R_i}$)，即给定区间的所有数的最大公约数，现在知道了每个问题和问题的答案。要构造出来一个数组A，满足给定的Q个问题，并且要求数组A所有元素的和最小，如果和最小的有多个，那么输出任意一个就可以。

首先，数组A每个元素最小是1，对于每个问题的区间，我们来构造和最小的序列使得最大公约数是给定的数，显然一种特殊的情况是区间中的数都是最大公约数，那么和一定最小；但是为了满足其他区间，就必须将区间中的数置为现在的数和当前问题答案的最小公倍数，全部构造完成后，再来检查一下每个问题的答案，是否满足要求，若满足，输出这个序列，否则，就输出“Stupid BrotherK!”。

代码：

/*************************************************************************
    > File Name: 5223.cpp
    > Author: gwq
    > Mail: gwq5210@qq.com 
    > Created Time: 2015年05月02日 星期六 22时19分57秒
 ************************************************************************/

#include <cmath>
#include <ctime>
#include <cctype>
#include <climits>
#include <cstdio>
#include <cstdlib>
#include <cstring>

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <string>
#include <vector>
#include <sstream>
#include <iostream>
#include <algorithm>

#define INF (INT_MAX / 10)
#define clr(arr, val) memset(arr, val, sizeof(arr))
#define pb push_back
#define sz(a) ((int)(a).size())

using namespace std;
typedef set<int> si;
typedef vector<int> vi;
typedef map<int, int> mii;
typedef long long ll;

const double esp = 1e-5;

#define N 1010

ll n, q;

ll num[N], l[N], r[N], s[N];

ll gcd(ll a, ll b)
{
    if (b == 0) {
        return a;
    } else {
        return gcd(b, a % b);
    }
}

int main(int argc, char *argv[])
{
    int t;
    scanf("%d", &t);
    while (t--) {
        cin >> n >> q;
        for (int i = 0; i < N; ++i) {
            num[i] = 1;
        }
        for (int i = 0; i < q; ++i) {
            cin >> l[i] >> r[i] >> s[i];
            for (int j = l[i]; j <= r[i]; ++j) {
                num[j] = (num[j] * s[i]) / gcd(num[j], s[i]);
            }
        }
        bool ok = true;
        for(int i = 0; i < q; i++){
            ll cur = num[l[i]];
            for(int j = l[i] + 1; j <= r[i]; j++){
                cur=gcd(cur,num[j]);
            }
            if(cur != s[i]){
                ok = false;
                break;
            }
        }
        if(ok){
            for(int i = 1; i <= n; i++){
                cout << num[i];
                if(i != n) {
                    printf(" ");
                }
            }
            printf("\n");
        }else{
            cout<<"Stupid BrotherK!"<<endl;
        }
    }
    return 0;
}

/*
2
3 2
1 2 12
2 3 6
3 2
1 2 5
2 3 6
*/