poj1062——昂贵的聘礼（dijkstra变形）

Description

年轻的探险家来到了一个印第安部落里。在那里他和酋长的女儿相爱了，于是便向酋长去求亲。酋长要他用10000个金币作为聘礼才答应把女儿嫁给他。探险家拿不出这么多金币，便请求酋长降低要求。酋长说：”嗯，如果你能够替我弄到大祭司的皮袄，我可以只要8000金币。如果你能够弄来他的水晶球，那么只要5000金币就行了。”探险家就跑到大祭司那里，向他要求皮袄或水晶球，大祭司要他用金币来换，或者替他弄来其他的东西，他可以降低价格。探险家于是又跑到其他地方，其他人也提出了类似的要求，或者直接用金币换，或者找到其他东西就可以降低价格。不过探险家没必要用多样东西去换一样东西，因为不会得到更低的价格。探险家现在很需要你的帮忙，让他用最少的金币娶到自己的心上人。另外他要告诉你的是，在这个部落里，等级观念十分森严。地位差距超过一定限制的两个人之间不会进行任何形式的直接接触，包括交易。他是一个外来人，所以可以不受这些限制。但是如果他和某个地位较低的人进行了交易，地位较高的的人不会再和他交易，他们认为这样等于是间接接触，反过来也一样。因此你需要在考虑所有的情况以后给他提供一个最好的方案。
为了方便起见，我们把所有的物品从1开始进行编号，酋长的允诺也看作一个物品，并且编号总是1。每个物品都有对应的价格P，主人的地位等级L，以及一系列的替代品Ti和该替代品所对应的”优惠”Vi。如果两人地位等级差距超过了M，就不能”间接交易”。你必须根据这些数据来计算出探险家最少需要多少金币才能娶到酋长的女儿。
Input

输入第一行是两个整数M，N（1 <= N <= 100），依次表示地位等级差距限制和物品的总数。接下来按照编号从小到大依次给出了N个物品的描述。每个物品的描述开头是三个非负整数P、L、X（X < N），依次表示该物品的价格、主人的地位等级和替代品总数。接下来X行每行包括两个整数T和V，分别表示替代品的编号和”优惠价格”。
Output

输出最少需要的金币数。
Sample Input

1 4
10000 3 2
2 8000
3 5000
1000 2 1
4 200
3000 2 1
4 200
50 2 0
Sample Output

5250

把每种物品看成节点，需要附加的金钱看成权值，很明显构造的是一个有向图，接下来就是求到点1的最短路径，因为不确定从谁那里开始买最省钱，所以再设一个虚拟起点0，这个点直接与其他所有点连接，权值为这个点所代表的物品的价格。
另外还有个等级的限制，枚举所有点的等级，不能与他交易的点都设为已访问就行

#include <iostream>
#include <cstring>
#include <string>
#include <vector>
#include <queue>
#include <cstdio>
#include <set>
#include <cmath>
#include <algorithm>
#define INF 0x3f3f3f3f
#define MAXN 10000010
#define Mod 10001
using namespace std;
int p[MAXN],a[MAXN],sq[MAXN],cnt;  //p保存素数，a判断当前数是否为素数，sq保存质数的平方数
int vis[MAXN]; //vis保存奇个数的质数
void Primer()
{
    int i,j, k;
    for(i=2;i<MAXN;++i)   //i从2开始遍历
    {
        if(a[i]==0)    //a[i]==0说明该i是素数
            p[cnt++]=i;   //cnt作为地址标志
        for(j=0;j<i&&(k=i*p[j])<MAXN;++j)   //j从0开始遍历，k作为有质数p[j]的合数
        {
            a[k]=1;    //合数设1
            if(i%p[j]==0)
                break;
        }
    }
}
int sum(int x)
{
    return x*(x-1)/2; //组合公式的化简
}
int main()
{
    cnt=0;
    Primer();
    int t,n,top=0,x;
    for(int i=0;i<cnt;++i)
    {
        sq[top++]=p[i]*p[i];
    }
    cin>>t;
    while(t--)
    {
        memset(vis,0,sizeof(vis));
        cin>>n;
        for(int i=0;i<n;++i)
        {
            cin>>x;
            for(int j=0;sq[j]<=x;++j)
            {
                while(x%sq[j]==0)
                    x/=sq[j];
            }
            vis[x]++;
        }
        int ans=0;
        for(int i=0;i<MAXN;++i)
        {
            ans+=sum(vis[i]);
        }
        cout<<ans<<endl;
    }
    return 0;
}