转载请注明出处,谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526 by---cxlove
题目:给出一个有向图,从1到n,每个结点有个权值,每走一步,分值为结点权值的LCM,而且每一步的LCM都要有变化,问到达N的时候分值恰好为K的路径有多少条
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4834
练习的时候,把K质因子分解,最多只有8个,也就是说路径最多8个结点,傻叉了,连LCM都搞错,把LCM(4,5)=LCM(2,5)了,果断跪了
其实不是考虑K的质因子,每次LCM的变化都为K的约数,所以求出K的所有约数,然后map离散一下
K最大为10^6,约数数目最多为2*sqrt(K);
dp[i][j]表示到达第i个结点的时候,当前为状态j的时候的路径数,记忆化搜索
#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<set>
#include<string>
#include<queue>
#define inf 1<<29
#define M 200005
#define N 2005
#define maxn 300005
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mem(a,b) memset(a,b,sizeof(a))
#define LL long long
#define MOD 1000000007
#define lson step<<1
#define rson step<<1|1
using namespace std;
struct Node
{
int v,next;
}edge[M];
int n,m,k;
int start[N],cnt=0,tot=0;
int val[N],dp[N][N];
map<LL,int>fac;
void addedge(int u,int v)
{
edge[cnt].v=v;
edge[cnt].next=start[u];
start[u]=cnt++;
}
LL gcd(LL a,LL b)
{
return b==0?a:gcd(b,a%b);
}
LL lcm(LL a,LL b)
{
return a/gcd(a,b)*b;
}
LL dfs(int u,LL x)
{
if(dp[u][fac[x]]!=-1) return dp[u][fac[x]];
if(u==n)
return x==k;
dp[u][fac[x]]=0;
for(int i=start[u];i!=-1;i=edge[i].next)
{
int v=edge[i].v;
if(k%val[v]) continue;
LL tmp=lcm(x,val[v]);
if(tmp>k||tmp==x) continue;
dp[u][fac[x]]=(dp[u][fac[x]]+dfs(v,tmp))%MOD;
}
return dp[u][fac[x]];
}
int main()
{
while(scanf("%d%d%d",&n,&m,&k)!=EOF)
{
fac.clear();
mem(start,-1);cnt=0;tot=0;
for(int i=1;i*i<=k;i++)
{
if(k%i==0)
{
fac[(LL)i]=++tot;
if(k/i!=i)
fac[(LL)k/i]=++tot;
}
}
while(m--)
{
int u,v;
scanf("%d%d",&u,&v);
addedge(u,v);
}
for(int i=1;i<=n;i++) scanf("%d",&val[i]);
if(fac.find(val[1])==fac.end()) printf("0\n");
else
{
mem(dp,-1);
printf("%lld\n",dfs(1,val[1]));
}
}
return 0;
}