hduoj 5223 GCD 构造检查
有一个数组 A1⋯AN ,数组中每一个元素是在区间 [1,109] 之间。现在有Q个问题,每一个问题是GCD( ALi, ALi+1, ALi+2, ..., ARi ),即给定区间的所有数的最大公约数,现在知道了每个问题和问题的答案。要构造出来一个数组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
*/