zoj 3644 记忆化搜索

题目：给出一个有向图，从1到n，每个结点有个权值，每走一步，分值为结点权值的LCM，而且每一步的LCM都要有变化，问到达N的时候分值恰好为K的路径有多少条

记忆化搜索，虽然做过很多了，但是一直比较慢，这次总结出几点

1.注意确定终点状态

2.状态的初始化

3.不可能状态的排除

代码是参考cxlove写的，kuangbin博客要是刷完了，下一个就刷他啦

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 #include<map>
 using namespace std;
 #define MOD 1000000007
 #define for0n for(i=0;i<n;i++)
 #define for1n for(i=1;i<=n;i++)
 #define for0m for(i=0;i<m;i++)
 #define for1m for(i=1;i<=m;i++)
 #define cl(a) memset(a,0,sizeof(a))
 #define ts printf("*****\n");
 const int MAXN=;
 const int MAXM=;
 int n,m,tt,k;
 int dp[MAXN][MAXN];
 struct Edge
 {
     int to,next;
 }edge[MAXM];
 int head[MAXN],a[MAXN];
 int tot;
 void addedge(int u,int v)
 {
     edge[tot].to=v;
     edge[tot].next=head[u];
     head[u]=tot++;
 }
 void init()
 {
     tot=;
     memset(head,-,sizeof(head));
 }
 map<int,int> fac;
 int gcd(int a,int b)
 {
     return b==?a:gcd(b,a%b);
 }
 int lcm(int a,int b)
 {
     return a/gcd(a,b)*b;
 }
 int dfs(int u,int val)
 {
     if(dp[u][fac[val]]!=-) return dp[u][fac[val]];
     if(u==n)    return val==k;
     dp[u][fac[val]]=;
     for(int i=head[u];i!=-;i=edge[i].next)
     {
         int v=edge[i].to;
         if(k%a[v])  continue;
         int temp=lcm(val,a[v]);
         if(temp>k||temp==val) continue;
         dp[u][fac[val]]=(dp[u][fac[val]]+dfs(v,temp))%MOD;
     }
     return dp[u][fac[val]];
 }
 int main()
 {
     int i,j;
     #ifndef ONLINE_JUDGE
     freopen("1.in","r",stdin);
     #endif
     while(scanf("%d%d%d",&n,&m,&k)!=EOF)
     {
         int u,v;
         int q=;
         init();
         for0m
         {
             scanf("%d%d",&u,&v);
             addedge(u,v);
         }
         for(i=;i*i<=k;i++)
         {
             if(k%i==)
             {
                 fac[i]=q++;
                 if(k/i!=i)  fac[k/i]=q++;
             }
         }
         for1n
         {
             scanf("%d",&a[i]);
         }
         if(fac.find(a[])==fac.end())   printf("0\n");
         else
         {
             memset(dp,-,sizeof(dp));
             printf("%d\n",dfs(,a[]));
         }
     }
 }