1 #include<stdio.h>
 2 #include<algorithm>
 3 #include<iostream>
 4 #include<string.h>
 5 #include<queue>
 6 #include<set>
 7 #include<math.h>
 8 #include<map>
 9 using namespace std;
10 typedef struct node
11 {
12         int m[4][4];
13         node()
14         {
15                 memset(m,0,sizeof(m));
16         }
17 } maxtr;
18 void Init(maxtr *ans,int x,int y);
19 maxtr E();
20 maxtr quick_m(maxtr ak,int m);
21 const int mod = 10007;
22 int main(void)
23 {
24         int n,x,y;
25         while(scanf("%d %d %d",&n,&x,&y)!=EOF)
26         {
27                 int f1 = 2;
28                 int a1 = 1;
29                 int a0 = 1;
30                 int xx = 1;
31                 maxtr ask ;
32                 Init(&ask,x,y);
33                 maxtr tp = quick_m(ask,n-1);
34                 printf("%d\n",(tp.m[0][0]*2+tp.m[0][1]*a1+tp.m[0][2]*a0+tp.m[0][3]*xx)%mod);
35         }
36         return 0;
37 }
38 void Init(maxtr *ans,int x,int y)
39 {      memset(ans->m,0,sizeof(ans->m));
40         x%=mod;y%=mod;
41         ans->m[0][0] = 1;
42         ans->m[0][1] = x*x%mod;
43         ans->m[0][2] = y*y%mod;
44         ans->m[0][3] =2*x*y%mod;
45         ans->m[1][1] = x*x%mod;
46         ans->m[1][2] = y*y%mod;
47         ans->m[1][3] = 2*x*y%mod;
48         ans->m[2][1] = 1;
49         ans->m[3][1] = x%mod;
50         ans->m[3][3] = y%mod;
51 }
52 maxtr E()
53 {
54         maxtr ak;
55         int i,j;
56         for(i = 0; i < 4; i++)
57         {
58                 for(j = 0; j < 4; j++)
59                 {
60                         if(i == j)
61                                 ak.m[i][j]  = 1;
62                 }
63         }
64         return ak;
65 }
66 maxtr quick_m(maxtr ak,int m)
67 {
68         int i,j;
69         maxtr ac = E();
70         while(m)
71         {
72                 if(m&1)
73                 {
74                         maxtr a;
75                         for(i = 0; i < 4; i++)
76                                 for(j = 0; j < 4; j++)
77                                         for(int s= 0; s < 4; s++)
78                                                 a.m[i][j] = (a.m[i][j] + ak.m[i][s]*ac.m[s][j]%mod)%mod;
79                         ac = a;
80                 }
81                 maxtr b;
82                 for(i = 0; i < 4; i++)
83                         for(j = 0; j < 4; j++)
84                                 for(int s = 0; s < 4; s++)
85                                         b.m[i][j] = (b.m[i][j] + ak.m[i][s]*ak.m[s][j]%mod)%mod;
86                 ak = b;
87                 m>>=1;
88         }
89         return ac;
90 }