Codeforces Round #785 (Div. 2) D Lost Arithmetic Progression

题意：给定两个等差数列B,C 其中C是A和B的所有的公共部分 求这样的等差数列A的个数

若无限个A 输出-1

细节题

思路：

首先判断C的项是不是全部在B中 若不是 答案为0 直接输出

然后A必须有所有C的项 所以A的公差必然是C的公差的因子

记 Di 为 i 的公差

即 Dc % Da == 0

想到 A和B不能有别的相同的项 且都是等差数列

那么 C的两项之间 ,A，B的每一项没有交集

所以lcm（Da，Db）== Dc

还有无限个A的情况

C的后一项/前一项如果在B中没有出现过 那么只要答案非0（前面判过0了） 都是无限个

因为A数组可以往左/右方向无限扩展

至于有效贡献的计算 就是在C首项左边，末项右边分别可以添加多少个A的项

从（0~Dc/Da-1） 所以就是ans+=(Dc/Da)^2

特殊情况

公差为1时只有在Db==Dc时 可能有贡献

Da==Dc时 除了-1的情况 只提供+1的贡献

const int mod = 1e9 + 7;
const int INF = 0x3f3f3f3f;
const int N = 1000010;
int n, m, k, ans, a, b, c, T, cnt, tot, top, num, sum, root, mas, Q, x, y, t;
char ch;

int u[N];
void dd(int x) {
    m = 0;
    int i = 2;
    for (; i * i < x; i++) {
        if (x % i == 0) {
            u[++m] = i;
            u[++m] = x / i;
        }
    }
    if (x == i * i) u[++m] = i;
}
int lcm(int x, int y) {
    return x / __gcd(x, y) * y;
}

signed main() {
    cf {
        ans = 0;
        a = read();
        b = read();
        c = read();  //首项/公差/项数
        int mx1 = a + (c - 1) * b;
        x = read();
        y = read();
        t = read();
        int mx2 = x + (t - 1) * y;

        //先判断b是不是全部在c中
        if (y % b != 0 || mx2 > mx1 || x < a || b > y || (x - a) % b != 0 || (mx2 - a) % b != 0) {
            cout << 0 << "\n";
            continue;
        }
        //再判断无限的情况
        if (x - y < a || mx2 + y > mx1) {
            cout << -1 << "\n";
            continue;
        }
        dd(y);
        rep(i, 1, m) {
            int now = lcm(b, u[i]);
            if (now != y) { continue; }
            ans += y * y / u[i] / u[i];
            ans %= mod;
        }
        ans++;
        if (y == b) { ans += y * y; }
        ans %= mod;
        cout << ans << "\n";
    }
    return 0;
}