植物杂交

题目:
作物杂交是作物栽培中重要的一步。

已知有 N 种作物 (编号 1 至 N)，第 i 种作物从播种到成熟的时间为 Ti。

作物之间两两可以进行杂交，杂交时间取两种中时间较长的一方。

如作物 A 种植时间为 5 天，作物 B 种植时间为 7 天，则 AB 杂交花费的时间为 7 天。

作物杂交会产生固定的作物，新产生的作物仍然属于 N 种作物中的一种。

初始时，拥有其中 M 种作物的种子 (数量无限，可以支持多次杂交)。

同时可以进行多个杂交过程。

求问对于给定的目标种子，最少需要多少天能够得到。

如存在 4 种作物 ABCD，各自的成熟时间为 5 天、7 天、3 天、8 天。

初始拥有 AB 两种作物的种子，目标种子为 D，已知杂交情况为 A×B→C，A×C→D。

则最短的杂交过程为：

第 1 天到第 7 天 (作物 B 的时间)，A×B→C。

第 8 天到第 12 天 (作物 A 的时间)，A×C→D。

花费 12 天得到作物 D 的种子。

输入格式
输入的第 1 行包含 4 个整数 N,M,K,T，N 表示作物种类总数 (编号 1 至 N)，M 表示初始拥有的作物种子类型数量，K 表示可以杂交的方案数，T 表示目标种子的编号。

第 2 行包含 N 个整数，其中第 i 个整数表示第 i 种作物的种植时间 Ti。

第 3 行包含 M 个整数，分别表示已拥有的种子类型 Kj，Kj 两两不同。

第 4 至 K+3 行，每行包含 3 个整数 A,B,C，表示第 A 类作物和第 B 类作物杂交可以获得第 C 类作物的种子。

输出格式
输出一个整数，表示得到目标种子的最短杂交时间。

样例解释
1≤N≤2000,
2≤M≤N,
1≤K≤105,
1≤T≤N,
1≤Ti≤100,
1≤Kj≤M,
保证目标种子一定可以通过杂交得到。
不保证作物 A 和 B 杂交只能生成作物 C（也就是说，A×B→C 和 A×B→D 可能同时在输入中出现）
不保证作物 C 只能由作物 A 和 B 杂交生成（也就是说，A×B→D 和 A×C→D 可能同时在输入中出现）。
不保证同一杂交公式不在输入中重复出现。

建边： 若 AXB->C
建边： A->C,B->C,边权为max(A,B)
而且在建边的时候 w[i],设置为PII
既要储存边的大小，也要储存是和哪个点杂交的边

#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
typedef pair<int,int> PII;
const int N = 2010,M = 2e5+10;
bool st[N];
int e[M],ne[M],h[N],idx;
PII w[M];
int a[N];
bool sq[N];
int d[N];
int n,m,k,t;
vector<int> g[N][N];
void add(int a,int b,int c,int d)
{
    e[idx] = b,ne[idx] = h[a],w[idx] = {c,d},h[a] = idx++;
}

int spfa()
{
    queue<int> q;
    memset(d,0x3f,sizeof d);
    for(int i= 1;i<=n;i++)
    if(sq[i])
    {
        q.push(i);
        st[i] = true;
        d[i] = 0;
    }
    while(q.size())
    {
        int ver =q.front();
        q.pop();
        st[ver]  = false;
        for(int i  = h[ver];~i;i=ne[i])
        {
            int j = e[i];
          int distance = w[i].first,x = w[i].second;
            if(sq[x]&&sq[ver] &&d[j]>max(d[ver],d[x])+distance)
            {
                d[j] = max(d[ver],d[x])+distance;
                sq[j] = true;
               if(!st[j])
               {
                   q.push(j);
                   st[j] = true;
               }
            }
            
        }
    }
    
    return d[t];
    
}
int main()
{
    
    memset(h,-1,sizeof h);
    scanf("%d%d%d%d",&n,&m,&k,&t);
    for(int i=  1;i<=n;i++) scanf("%d",&a[i]);
    for(int i = 1;i<=m;i++)
    {
        int x;
        scanf("%d",&x);
        sq[x]= true;
    }
    for(int i = 0;i<k;i++)
    {
        int x,b,c;
        scanf("%d%d%d",&x,&b,&c);
        int y = max(a[x],a[b]);
        add(x,c,y,b),add(b,c,y,x);
    }
    printf("%d\n",spfa());
    return 0;
}