Use openmp reduction to parallelize the sample code below

#include

#include

#include

#include

int main (int argc, char *argv[])

{

int i, n;

float a[1000000], b[1000000], sum;

/* Some initializations */ n = 1000000;

for (i=0; i < n; i++)

a[i] = b[i] = i * 1.0; sum = 0.0;

// Insert a parallel code here

for (i=0; i < n; i++)

sum = sum + (a[i] * b[i]*b[i]*b[i]);

printf(" Sum = %f ",sum);

}