题意:给定两个等差数列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;
}