第九届ECNU Coder K.计软联谊

题目链接：http://acm.ecnu.edu.cn/contest/16/problem/K/

题目：

K. 计软联谊

Time limit per test: 7.0 seconds

Time limit all tests: 7.0 seconds

Memory limit: 512 megabytes

Accept / Submit: 19 / 398

在计算机和软件专业的联谊会上，计算机和软件的同学相间着排成一列。现在要计算相邻两个同学的友谊度。

友谊度 friend(a,b) 是这么计算的：令 a, b 两个整数分别是两个同学的属性，两个同学的友谊度取决于 a,b 第 k 大的公约数。如果不存在，就说明这两个同学之间完全没有友谊，友谊度为 −1。

Input

第一行是数据组数 T (1≤T≤60)。

对于每组数据：
第一行：首先是学生的数量 n (1≤n≤105)，约定的常数 k (1≤k≤106)。
第二行：n 个整数，依次表示这些学生的属性值：m1,m2,…,mn (1≤mi≤106)。

Output

对于每组数据输出一行，以 Case x: 开头（x 表示数据编号，从 1 开始），后面是 n−1 个整数，分别是 friend(m1,m2),friend(m2,m3),…,friend(mn−1,mn)，整数和整数之间用空格隔开。

Examples

input

2
3 1
4 6 12
6 2
13 12 12 24 36 30

output

Case 1: 2 6
Case 2: -1 6 6 6 3

Note

请注意输入输出上的优化！

题解：

　　就每两个数gcd，求出最大公约数。然后求第k大的公约数。

　　一开始理解错，以为第2大的公约数是最大公约数除以2，第3大则是除以4。dalao说是错的。

　　第k大的公约数是最大公约数的约数（不一定是除以2，可以除以3什么的）。

　　就先用筛法将1~1e6的每个数的约数筛出来。

　　

 void init() {
     for(int i = ;i<=L;i++){
         for(int j = i;j<=L;j+=i){
             P[j].push_back(i);
         }
     }
 }

筛法

　　详细代码：

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <cmath>
 #include <string>
 #include <vector>
 #include <map>
 #include <set>
 #include <stack>
 #include <queue>
 #include <sstream>
 #include <algorithm>
 using namespace std;
 #define pb push_back
 #define mp make_pair
 #define ms(a, b)  memset((a), (b), sizeof(a))
 #define eps 0.0000001
 typedef long long LL;
 typedef unsigned long long ULL;
 const int inf = 0x3f3f3f3f;
 const LL INF = 0x7fffffff;
 const int maxn = 1e5+;
 const int L = 1e6;
 const int mod = 1e9+;
 int m[maxn];
 vector <int > P[L+];
 int kase = ;
 void init() {
     for(int i = ;i<=L;i++){
         for(int j = i;j<=L;j+=i){
             P[j].push_back(i);
         }
     }
 }
 int gcd(int a, int b){
     return b==?a:gcd(b, a%b);
 }
 void solve() {
     int n, k;
     scanf("%d%d", &n, &k);
     for(int i = ;i<=n;i++)
         scanf("%d", &m[i]);
     printf("Case %d:", kase++);
     for(int i = ;i<n;i++){
         int x = gcd(m[i], m[i+]);
         int len = P[x].size();
         if(len-k>=){
             printf(" %d", P[x][len-k]);
         }
         else
             printf(" -1");
     }
     printf("\n");
 }
 int main() {
 #ifdef LOCAL
     freopen("input.txt", "r", stdin);
 //      freopen("output.txt", "w", stdout);
 #endif // LOCAL
     init();
     int T;
     scanf("%d", &T);
     while(T--){
         solve();
     }
     return ;
 }

你努力的时候，比你厉害的人也在努力。