Cod sursa(job #1195260)

Utilizator apopeid15Apopei Daniel apopeid15 Data 6 iunie 2014 19:55:00
Problema Atac Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.14 kb
#include <fstream>
#include <cmath>
#include <algorithm>
#include <vector>
#include <utility>
 
#define inf 300001
#define vintp vector<pair<int,int> >::iterator
 
using namespace std;
 
ifstream fin ("atac.in");
ofstream fout("atac.out");
 
vector <pair<int,int> >T[32001];
int rmq[16][32001],t[16][32001],logn,L[32001];
int x,y,a,b,c,d,n,m,p,z;
 
void dfs (int x, int father, int lvl)
{
    L[x] = lvl;
 
    for (vintp it = T[x].begin (); it != T[x].end (); ++it)
    {
        if (it->first != father)
        {
           t[0][it->first] = x;
           rmq[0][it->first] = it->second;
           dfs (it->first,x,lvl+1);
        }
    }
}
 
void RMQ ()
{
    for (int i=1; i<=logn; ++i)
    {
        for (int j = 1; j<= n; ++j)
        {
            t[i][j] = t[i-1][t[i-1][j]];
            rmq[i][j] = min (rmq[i-1][j],rmq[i-1][t[i-1][j]]);
        }
    }
}
 
int query (int x, int y)
{
    int val = inf;
 
    if (L[x] != L[y])
    {
        if (L[x] < L[y])
          swap (x,y);
 
        int op = L[x] - L[y];
 
        for (int i=0; i<=logn; ++i)
        {
            if ((op>>i)&1)
            {
                val = min (val,rmq[i][x]);
                x = t[i][x];
            }
        }
    }
 
    if (x == y)
    {
      if (val == inf)
        return 0;
        return val;
    }
 
 
    for (int i = logn; i>=0; --i)
    {
        if (t[i][x] != t[i][y])
        {
            val = min (val,min(rmq[i][x],rmq[i][y]));
            x = t[i][x];
            y = t[i][y];
        }
    }
 
 
    return min (val,min(rmq[0][x],rmq[0][y]));
}
 
int main()
{
    fin>>n>>m>>p;
 
    for (int i=2; i<=n; ++i)
    {
        fin>>x>>y;
 
        T[i].push_back (make_pair(x,y));
        T[x].push_back (make_pair(i,y));
    }
 
    rmq[0][1] = inf;
    dfs (1,0,1);
    logn = log2(n);
    RMQ ();
 
    fin>>x>>y>>a>>b>>c>>d;
 
    for (int i=1; i<=m; ++i)
    {
        z = query (x,y);
 
        if (i >= m-p+1)
        {
            fout<<z<<"\n";
        }
 
        x = (x*a + y*b) % n + 1;
        y = (y*c + z*d) % n + 1;
    }
 
    return 0;
}