#include <stdio.h>
#include <omp.h>

#define N 1000000

int main() {
    double sum = 0.0;
    double a[N], b[N];

    // Initialize vectors a and b

    #pragma omp parallel for reduction(+:sum)
    for (int i = 0; i < N; ++i) {
        sum += a[i] * b[i];
    }

    printf("Dot product: %f\n", sum);

    return 0;
}