#######################################################################
#   Copyright (C) 2015 All rights reserved.
#   
#   文件名称:p.pl
#   创 建 者:
#   创建日期:2015年05月11日
#   描    述:
#
#   备    注:
#
#########################################################################

use strict;
use warnings;
use Time::HiRes;
# please add your code here!
my $tmStarted=Time::HiRes::time;
if (2 != scalar @ARGV)
{
    &PrintUsage();
    exit 1;
}
=pod
    Subroutine(s);
=cut
sub PrintUsage
{
    print STDERR "perl program.pl [IN] ref_file [IN] input_file [OUT] output_file\n";
}
my $line = "";
my $linecount = 0;
my %hash = ();
open FIN, "<$ARGV[0]" or die "can not open inputfile:$!\n";
open FOUT, ">$ARGV[1]" or die "can not create outputfile :$!\n";
while(defined($line=<FIN>))
{
    chomp $line;
    my @vec = split /\t/, $line;
    my $size = scalar @vec;
    if ($size < 3)
    {
        print STDERR "$line\n";
    }
    if (not defined($hash{$vec[0]}))
    {
        my @vectemp = ();
        push @vectemp, $vec[2];
        $hash{$vec[0]} = \@vectemp;
    }
    else
    {
        push @{$hash{$vec[0]}}, $vec[2];
    }

}
close FIN;
foreach my $key (keys%hash)
{
    my @vec2 = @{$hash{$key}};
    my $size = scalar @vec2;
    my $part = "";
    for (my $i = 0; $i <$size; $i++)
    {
        $part.=$vec2[$i];
        if ($i < $size - 1)
        {
            $part.=",";
        }
    }
    print FOUT "$key\t$part\n";
}
close FOUT;
print STDERR "$0 has finished,congratulations!\n";
print STDERR "Time elapsed:".(Time::HiRes::time-$tmStarted)."\n";

 

posted on 2015-05-11 18:56  finallyly  阅读(587)  评论(0编辑  收藏  举报