#include <bits/stdc++.h>

using namespace std;
struct csucs
{
    int cs, p;
};
vector<vector<int>>vec;
int osv[100000];
bool veg[100000], volt[100000];
int req(int n, int el)
{
    int sum=0;
    for(auto i:vec[el])sum+=req(osv[i], i);
    if(sum >= n)
    {
        volt[el] = 1;
        for(auto i:vec[el])if(!volt[i])
        {
            veg[i]=1;
            volt[i]=1;
        }
        return sum;
    }
    else
    {
        return n;
    }
}
void tis(int n, bool van)
{
    if(van)
    {
        veg[n] = 0;
        for(auto i:vec[n])tis(i, van);
    }
    else
    {
        if(veg[n])van =1;
        for(auto i: vec[n])tis(i, van);
    }
}
int main()
{
    int i, n, a, j;
    cin >> n;
    csucs e;
    vec.resize(n);
    for(i=0; i<n; i++)cin >> osv[i];
    for(i=1; i<n; i++)
    {
        cin >> a;
        vec[a-1].push_back(i);
    }
    cout << req(osv[0], 0) << endl;
    a=0;
    tis(0, 0);
    for(i=0; i<n; i++)if(veg[i])a++;
    cout << a<<endl;
    for(i=0; i<n; i++)if(veg[i])cout << i+1 << ' ';
    return 0;
}