Intel Code Challenge Elimination Round (Div.1 + Div.2, combined)F. Cyclic Cipher

题意： 
给n个序列，一个序列内没有重复元素，不同序列内可能有重复元素，每个序列长度<=40，每秒钟每个序列会循环左移一次，要求算出 10100 秒内，[1,m]的答案。 
一个数x在某时刻的答案为，把此时所有序列的第一个元素顺序排成一个新序列，这个新序列中最长的连续x的长度。 
思路：考虑对m个答案分开求，设当前要求解的数为x，那么我们只用关心所有x的位置，可能的答案一定是连续出现x的序列，假设i，j是答案中的两个序列，长度为li，lj，对应x的位置为pi，pj，显然要让pi和pj同时出现在第一个元素，则有 pi+li∗k==pj+lj∗k ，然后有 pi−pj==lj∗k−li∗k ，两边同时对gcd(li,lj)取模，则有 (pi−pj)=0(modgcd(li,lj))

，然后可以得到一个结论，两个长度相同的序列，其pi和pj一定相同。因为序列长度<=40，所以长度相同的序列肯定是很多的，一个序列出现的数是很少的，我们就可以考虑two pointer了，过程中要保证这一段出现的x是连续的

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long ll;
const int N = 1e5+5;
int n, m;
struct node{
    int r, len, c;
    node(int _a = 0, int _b = 0, int _c = 0) {
        r = _a, len = _b, c = _c;
    }
};
vector<node>l[N];
int pos[50];
int cnt[50];
int gcd[50][50];
int _gcd(int a, int b){ return b? _gcd(b, a%b) : a; }
// (xi+a*li) == (xj+b*lj)
// xi-xj == b*lj-a*li
// (xi-xj)%gcd(lj,li) == 0
bool add(node a) {
    if(cnt[a.len] && pos[a.len] != a.c) return 0;
    for(int i = 1; i <= 40; ++i) {
        if(cnt[i] && (pos[i]-a.c)%gcd[i][a.len] != 0) return 0;
    }
    pos[a.len] = a.c;
    cnt[a.len]++;
    return 1;
}
void del(node a) {
    cnt[a.len]--;
}
int solve(vector<node>& x) {
    int res = 0;
    memset(pos, 0, sizeof(pos));
    memset(cnt, 0, sizeof(cnt));
    int l = 0, r = 0;
    while(l < x.size()) {
        while(r < x.size()) {
            if(x[r].r-x[l].r != r-l) break;
            if(!add(x[r])) break;
            ++r;
        }
        res = max(res, r-l);
        if(x[r].r - x[l].r != r-l) {
            l = r;
            memset(pos, 0, sizeof(pos));
            memset(cnt, 0, sizeof(cnt));
        }
        else del(x[l++]);
    }
    return res;
}
int main(int argc, const char * argv[]) {
    for(int i = 1; i < 50; ++i) {
        for(int j = 1; j < 50; ++j) {
            gcd[i][j] = _gcd(i, j);
        }
    }
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; ++i) {
        int k;
        scanf("%d", &k);
        for(int x, j = 1; j <= k; ++j) {
            scanf("%d", &x);
            l[x].emplace_back(node(i, k, j));
        }
    }
    for(int i = 1; i <= m; ++i) {
        printf("%d\n", solve(l[i]));
    }
    return 0;
}