csp2021-04-4 校门外的树

先预处理因子，用1e5个vector保存每个数字所有因子(不包括数字本身)
然后i和j双重循环
dp[i]表示从第0个障碍到第i个障碍的方案数量
障碍i固定，j从i-1到1 遍历两个障碍间长度的所有因子
对于每一个i，用一个set来保存已经使用过的约数
如果第j个障碍到第i个障碍间存在一个约数符合条件，那么在j前面的障碍到i之间，这个约数就不再符合条件

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <vector>
#include <map>
#include <queue>
#include <algorithm>
#include <math.h>
#include <cstdio>
#include <set>
#define MIN(x,y) (((x)<(y))?(x):(y))
#define MAX(x,y) (((x)>(y))?(x):(y))
typedef long long ll;
using namespace std;
const ll mod = 1e9+7;
const ll MAXN =100000;
ll ans=0;
int n;
ll b[1005];
ll dp[1005];
vector <vector <ll> > f(MAXN);
vector <ll> ::iterator it;
set <ll> st;
void fact()
{
	for(ll i=1; i<=MAXN/2; i++)
		{
			for(ll j=i*2; j<=MAXN; j+=i)
				f[j].push_back(i);
		}
}
ll solve(ll x,ll y)
{
	ll left=b[x];
	ll right=b[y];
	ll len=right-left;
	ll res=0;
	st.insert(len);
	for(it=f[len].begin(); it!=f[len].end(); it++)
		{
			if(st.find(*it)!=st.end())
				continue;
			res++;
			st.insert(*it);
		}

	return res;
}
int main()
{
	fact();
	memset(dp,0,sizeof(dp));
	scanf("%d",&n);
	for(int i=0; i<n; i++)
		scanf("%lld",b+i);
	dp[0]=1;
	for(int i=1; i<n; i++)
		{
			st.clear();
			for(int j=i-1; j>=0; j--)
				{
					dp[i]+=(dp[j]*solve(j,i))%mod;
					dp[i]%=mod;
				}
		}
	printf("%lld\n",dp[n-1]);
}