Problem: Matryoshka

Matryoshka is a famous Russian folk craft doll. Matryoshka can be divided into upper and lower parts, and when opened, there is another smaller doll inside. The nesting structure is such that when the small doll that appears is opened, a smaller doll is contained.

<image>


You found an unusually shaped matryoshka doll on your trip and bought N dolls. The shape of the i-th doll is a rectangular parallelepiped of xi × yi × zi.

After watching Matryoshka for a while, you are about to put away Matryoshka. Before that, I want to reduce the space required by storing some dolls in another. When storing a doll, one other doll can be stored only in the doll that has not yet stored any doll. However, only the dolls that are stored directly are counted, and the doll that contains the doll can be stored in another doll.

The stored doll becomes invisible from the outside. However, the following conditions must be met.

* The doll may rotate, but each side of the rectangular parallelepiped is parallel to any side of the other rectangular parallelepiped.
* After rotation, the length of the doll on the side to be stored is shorter for each of the lengths of the corresponding sides.
* At most one doll can be stored directly in one doll



Since the volume of the closet is limited, we want to minimize the sum of the volumes of the dolls that can be seen from the outside. Your job is to create a program that finds the minimum sum of the volumes of the dolls that are visible from the outside, which can be achieved by repeating the operation of storing the dolls any number of times.

Input

The input consists of multiple datasets. The maximum number of data sets does not exceed 50. Each dataset is represented in the following format.

> N
> x1 y1 z1
>::
>::
> xN yN zN
>

Each dataset consists of N + 1 rows, and the first row of the dataset is given the integer N, which represents the number of dolls. In the i-th line of the following N lines, three integers xi, yi, and zi representing the size of the i-th doll are given, separated by a half-width space. These integers satisfy 1 ≤ N, xi, yi, zi ≤ 100.

> The end of the input is represented by a single zero line.

> ### Output

For each data set, output the minimum value of the sum of the volumes of the dolls that can be seen from the outside in one line.

> ### Sample Input


2
one two Three
4 2 3
3
2 5 2
3 3 4
5 5 5
Five
1 1 1
2 2 2
3 3 3
4 4 4
5 5 5
Five
1 1 1
2 1 1
3 1 1
4 1 1
5 1 1
Ten
3 1 4
1 5 9
2 6 5
3 5 8
9 7 9
3 2 3
8 4 6
2 6 4
3 3 8
3 2 7
0


Output for Sample Input


twenty four
145
125
15
864





Example

Input

2
1 2 3
4 2 3
3
2 5 2
3 3 4
5 5 5
5
1 1 1
2 2 2
3 3 3
4 4 4
5 5 5
5
1 1 1
2 1 1
3 1 1
4 1 1
5 1 1
10
3 1 4
1 5 9
2 6 5
3 5 8
9 7 9
3 2 3
8 4 6
2 6 4
3 3 8
3 2 7
0


Output

24
145
125
15
864

Solution (in CPP):

Solution:
#include<bits/stdc++.h>
using namespace std;
#define ALL(x) x.begin(),x.end()
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}
 
struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
template<typename T1,typename T2>
ostream &operator<<(ostream &os,const pair<T1,T2>&p){
    os<<p.first<<" "<<p.second;
    return os;
}
 
template<typename T>
ostream &operator<<(ostream &os,const vector<T>&v){
    for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");
    return os;
}

template<typename T1,typename T2>
istream &operator>>(istream &is,pair<T1,T2>&p){
    is>>p.first>>p.second;
    return is;
}

template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}

template<typename flow_t, typename cost_t>
struct PrimalDual{
    const cost_t TINF;
 
    struct edge{
        int to;
        flow_t cap;
        cost_t cost;
        int rev;//この辺の逆辺がg[from]の何番目にあるか
        bool isrev;
    };
 
    vector<vector<edge>> graph;
    vector<cost_t> potential,min_cost;
    vector<int> prevv,preve;//点，辺
 
    PrimalDual(int V):graph(V),TINF(numeric_limits<cost_t>::max()){}
 
    void add_edge(int from,int to,flow_t cap,cost_t cost){
        graph[from].push_back((edge){to,cap,cost,(int)graph[to].size(),false});
        graph[to].push_back((edge){from,0,-cost,(int)graph[from].size()-1,true});
    }
 
    cost_t min_cost_flow(int s,int t,flow_t f){
        int V=(int)graph.size();
        cost_t ret=0;
        using Pi=pair<cost_t,int>;
        priority_queue<Pi,vector<Pi>,greater<Pi>> que;
        potential.assign(V,0);
        preve.assign(V,-1);
        prevv.assign(V,-1);
 
        while(f>0){
            min_cost.assign(V,TINF);
            que.emplace(0,s);
            min_cost[s]=0;
            //dijkstraパート
            while(!que.empty()){
                Pi p=que.top();que.pop();
                if(min_cost[p.second]<p.first) continue;
                for(int i=0;i<graph[p.second].size();i++){
                    edge &e=graph[p.second][i];
                    cost_t nextCost=min_cost[p.second]+e.cost+potential[p.second]-potential[e.to];
                    if(e.cap>0 and min_cost[e.to]>nextCost){
                        min_cost[e.to]=nextCost;
                        prevv[e.to]=p.second,preve[e.to]=i;
                        que.emplace(min_cost[e.to],e.to);
                    }
                }
            }
            if(min_cost[t]==TINF) return -1;
            //dijkstraの結果に応じてpotentialを上げ下げ
            for(int v=0;v<V;v++)potential[v]+=min_cost[v];
            flow_t addflow=f;
            for(int v=t;v!=s;v=prevv[v]){
                addflow=min(addflow,graph[prevv[v]][preve[v]].cap);
            }
            f-=addflow;
            ret+=addflow*potential[t];//天才じゃん
            for(int v=t;v!=s;v=prevv[v]){
                edge &e=graph[prevv[v]][preve[v]];
                e.cap-=addflow;
                graph[v][e.rev].cap+=addflow;
            }
        }
        return ret;
    }
 
    void output(){
        for(int i=0;i<graph.size();i++){
            for(auto &e:graph[i]){
                if(e.isrev)continue;
                auto &rev_e=graph[e.to][e.rev];
                cout<<i<<"->"<<e.to<<" (flow: "<<rev_e.cap<<" / "<<rev_e.cap+e.cap<<")"<<endl;
            }
        }
    }
};
 



int n;

struct Box{
    vector<int> p;
    int v;
    Box(int x,int y,int z){
        p.resize(3);
        p[0]=x,p[1]=y,p[2]=z;
        sort(ALL(p));
        v=x*y*z;
    }
    Box(){}
};

// can a cover b?
bool cover(Box a,Box b){
    rep(i,3){
        if(not(a.p[i]>b.p[i])) return false;
    }
    return true;
}

void solve(){
    vector<Box> v(n);
    ll res=0;
    rep(i,n){
        int x,y,z;cin>>x>>y>>z;
        v[i]=Box(x,y,z);
        res+=v[i].v;
    }

    PrimalDual<int,int> flow(2*n+2);
    int src=2*n,sink=2*n+1;
    flow.add_edge(src,sink,n,0);
    rep(i,n){
        flow.add_edge(src,i,1,0);
        flow.add_edge(n+i,sink,1,0);
    }

    rep(i,n){
        rep(j,n){
            if(cover(v[i],v[j])){
                flow.add_edge(i,n+j,1,-v[j].v);
            }
        }
    }
    cout<<res+flow.min_cost_flow(src,sink,n)<<endl;
}

signed main(){
    while(cin>>n,n) solve();
    return 0;
}