Codeforces Round #785 (Div. 2) D. Lost Arithmetic Progression （类似质因分解）

传送门：

link

题意：

• 给你两个等差数列的：首项，公差，项数。
• 数列$B$和数列$C$给出.要求有多少满足条件的数列$A$
• 数列$C$ 是 数列$A$ 和 数列$B$的公共项。

思路：

感觉难度在 $[1800,2100]$的样子。

这里为了方便说明：
令$B$的首项，公差，最后一项分别为 $b,d2,bn$
令$C$的首项，公差，最后一项分别为 $c,d3,cn$

先把边界特判了：

1. 由于数列$C$是公共项，所以$C$的范围肯定在$B$里面。
2. 在1的基础上，我们再看 $c(首项)$是否在$B$内，还要看$d3$是否是$d2$的倍数。
3. 之后就是无穷解的情况，假如 $c-d3<b$，那么我们可以在左边无限填充，另一边一样判断。

看完边界后，之后枚举所有可能的$d1$
因为 $lcm(d1,d2)=d3$（题目的条件，公共项），所以我们利用这个条件看枚举的$d1$是否合法。

这里可以反证：
假如$lcm(d1,d2)<d3$，那么与题目公共项条件相违背。

对于怎么找因子，可以在这篇里看oi-wiki

#include <iostream>
#include <cstdio>
#include <vector>
#include <map>
#include <unordered_map>
#include <set>
#include <algorithm>
#include <queue>
#include <stack>
#include <string>
#include <cmath>
#include <cstring>
#define For(i,x,y) for(int i = (x); i <= (y); i ++ )
#define fori(i,x,y) for(int i = (x); i < (y); i ++ )
#define sz(a) (int)a.size()
#define ALL(a) a.begin(), a.end()
#define mst(x,a) memset(x,a,sizeof(x))
#define pb push_back
#define eb emplace_back
#define mp make_pair
#define fi first
#define se second
#define db double
#define endl '\n' 
#define debug(a) cout << #a << ": " << a << endl
using namespace std;
typedef long long LL;
typedef long long ll;
typedef unsigned long long ULL;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int inf = 0x3f3f3f3f;
const int mod = 1e9+7;
typedef pair<int,int>pa;
typedef pair<ll,ll>pai;
typedef pair<db,db> pdd;
const db eps = 1e-6;
const db pi = acos(-1.0);

template<typename T1, typename T2> void ckmin(T1 &a, T2 b) { if (a > b) a = b; }
template<typename T1, typename T2> void ckmax(T1 &a, T2 b) { if (a < b) a = b; }
int read() {
    int x = 0, f = 0; char ch = getchar();
    while (!isdigit(ch)) f |= ch == '-', ch = getchar();
    while (isdigit(ch)) x = 10 * x + ch - '0', ch = getchar();
    return f ? -x : x;
}
template<typename T> void print(T x) {
    if (x < 0) putchar('-'), x = -x;
    if (x >= 10) print(x / 10);
    putchar(x % 10 + '0');
}
template<typename T> void print(T x, char let) {
    print(x), putchar(let);
}
template<class T> bool uin(T &a, T b) { return a > b ? (a = b, true) : false; }
template<class T> bool uax(T &a, T b) { return a < b ? (a = b, true) : false; }

const int maxn = 100000 + 5;

ll gcd(ll a, ll b) {
    if(b == 0) return a;
    return gcd(b,a%b);
}

ll lcm(ll a, ll b){
    ll d = gcd(a,b);
    return a/d*b;
}


void sol(){
    ll b, c, d2, d3, n2, n3;
    cin>>b>>d2>>n2;
    cin>>c>>d3>>n3;
    //考虑边界
    ll cn = c + (n3 - 1)*d3, bn = b + (n2 - 1)*d2;
    if(c < b ||  cn > bn ) {
        cout<<0<<endl;
        return ;
    }
    ll rem = (c - b);
    if(rem % d2 || d3 % d2) {
        cout<<0<<endl;
        return ;
    }

    if(c - d3 < b || cn + d3 > bn) {
        cout<<-1<<endl;
        return ;
    }

    ll ans = 0;
    for(ll i = 1; i * i <= d3; i ++ ) {
        if(d3 % i == 0 ){
            if(lcm(i,d2) == d3) {
                ans = ans + (d3/i)*(d3/i)%mod;
                ans %= mod;
            }
            if(d3/i != i && lcm(d3/i,d2) == d3) {
                ans = ans + (d3/(d3/i))*(d3/(d3/i))%mod;
                ans %= mod;
            }
        }
    }
    cout<<ans<<endl;
}



int main() {
    ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int tt; cin>>tt;
    while(tt--)
        sol();
    return 0;
}