shor.html 15 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314
  1. <!DOCTYPE HTML>
  2. <html lang="en" class="sidebar-visible">
  3. <head>
  4. <!-- Book generated using mdBook -->
  5. <meta charset="UTF-8">
  6. <title>Shors - QCGPU User Guide</title>
  7. <meta content="text/html; charset=utf-8" http-equiv="Content-Type">
  8. <meta name="description" content="The user guide / prose documentation for QCGPU">
  9. <meta name="viewport" content="width=device-width, initial-scale=1">
  10. <meta name="theme-color" content="#ffffff" />
  11. <base href="../">
  12. <link rel="stylesheet" href="book.css">
  13. <link href="https://fonts.googleapis.com/css?family=Open+Sans:300italic,400italic,600italic,700italic,800italic,400,300,600,700,800" rel="stylesheet" type="text/css">
  14. <link href="https://fonts.googleapis.com/css?family=Source+Code+Pro:500" rel="stylesheet" type="text/css">
  15. <link rel="shortcut icon" href="favicon.png">
  16. <!-- Font Awesome -->
  17. <link rel="stylesheet" href="https://maxcdn.bootstrapcdn.com/font-awesome/4.3.0/css/font-awesome.min.css">
  18. <link rel="stylesheet" href="highlight.css">
  19. <link rel="stylesheet" href="tomorrow-night.css">
  20. <link rel="stylesheet" href="ayu-highlight.css">
  21. <!-- Custom theme -->
  22. <!-- MathJax -->
  23. <script type="text/javascript" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.1/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
  24. <!-- Fetch Clipboard.js from CDN but have a local fallback -->
  25. <script src="https://cdn.jsdelivr.net/clipboard.js/1.6.1/clipboard.min.js"></script>
  26. <script>
  27. if (typeof Clipboard == 'undefined') {
  28. document.write(unescape("%3Cscript src='clipboard.min.js'%3E%3C/script%3E"));
  29. }
  30. </script>
  31. <noscript>
  32. <style type="text/css">
  33. .javascript-only {
  34. display: none;
  35. }
  36. </style>
  37. </noscript>
  38. </head>
  39. <body class="light">
  40. <!-- Work around some values being stored in localStorage wrapped in quotes -->
  41. <script type="text/javascript">
  42. try {
  43. var theme = localStorage.getItem('mdbook-theme');
  44. var sidebar = localStorage.getItem('mdbook-sidebar');
  45. if (theme.startsWith('"') && theme.endsWith('"')) {
  46. localStorage.setItem('mdbook-theme', theme.slice(1, theme.length - 1));
  47. }
  48. if (sidebar.startsWith('"') && sidebar.endsWith('"')) {
  49. localStorage.setItem('mdbook-sidebar', sidebar.slice(1, sidebar.length - 1));
  50. }
  51. } catch (e) { }
  52. </script>
  53. <!-- Set the theme before any content is loaded, prevents flash -->
  54. <script type="text/javascript">
  55. var theme;
  56. try { theme = localStorage.getItem('mdbook-theme'); } catch(e) { }
  57. if (theme === null || theme === undefined) { theme = 'light'; }
  58. document.body.className = theme;
  59. document.querySelector('html').className = theme;
  60. </script>
  61. <!-- Hide / unhide sidebar before it is displayed -->
  62. <script type="text/javascript">
  63. var html = document.querySelector('html');
  64. var sidebar = 'hidden';
  65. if (document.body.clientWidth >= 1080) {
  66. try { sidebar = localStorage.getItem('mdbook-sidebar'); } catch(e) { }
  67. sidebar = sidebar || 'visible';
  68. }
  69. html.classList.remove('sidebar-visible');
  70. html.classList.add("sidebar-" + sidebar);
  71. </script>
  72. <nav id="sidebar" class="sidebar" aria-label="Table of contents">
  73. <ol class="chapter"><li><a href="qcgpu.html"><strong aria-hidden="true">1.</strong> QCGPU</a></li><li><ol class="section"><li><a href="getting-started.html"><strong aria-hidden="true">1.1.</strong> Getting Started</a></li></ol></li><li><a href="user-guide/user-guide.html"><strong aria-hidden="true">2.</strong> User Guide</a></li><li><ol class="section"><li><a href="user-guide/registers.html"><strong aria-hidden="true">2.1.</strong> Quantum Registers</a></li><li><a href="user-guide/gates.html"><strong aria-hidden="true">2.2.</strong> Quantum Gates</a></li><li><a href="user-guide/operations.html"><strong aria-hidden="true">2.3.</strong> Quantum Operations</a></li><li><a href="user-guide/examples.html"><strong aria-hidden="true">2.4.</strong> Examples</a></li><li><a href="user-guide/decoherence.html"><strong aria-hidden="true">2.5.</strong> Decoherence</a></li></ol></li><li><a href="algorithms/algorithms.html"><strong aria-hidden="true">3.</strong> Algorithms</a></li><li><ol class="section"><li><a href="algorithms/bernstein-vazirani.html"><strong aria-hidden="true">3.1.</strong> Bernstein-Vazirani</a></li><li><a href="algorithms/deutsch-jozsa.html"><strong aria-hidden="true">3.2.</strong> Deutsch-Jozsa</a></li><li><a href="algorithms/grover.html"><strong aria-hidden="true">3.3.</strong> Grovers</a></li><li><a href="algorithms/shor.html" class="active"><strong aria-hidden="true">3.4.</strong> Shors</a></li><li><a href="algorithms/super-dense.html"><strong aria-hidden="true">3.5.</strong> Super Dense Coding</a></li></ol></li></ol>
  74. </nav>
  75. <div id="page-wrapper" class="page-wrapper">
  76. <div class="page">
  77. <div id="menu-bar" class="menu-bar">
  78. <div id="menu-bar-sticky-container">
  79. <div class="left-buttons javascript-only">
  80. <button id="sidebar-toggle" class="icon-button" type="button" title="Toggle Table of Contents" aria-label="Toggle Table of Contents" aria-controls="sidebar">
  81. <i class="fa fa-bars"></i>
  82. </button>
  83. <button id="theme-toggle" class="icon-button" type="button" title="Change theme" aria-label="Change theme" aria-haspopup="true" aria-expanded="false" aria-controls="theme-list">
  84. <i class="fa fa-paint-brush"></i>
  85. </button>
  86. <ul id="theme-list" class="theme-popup" aria-label="Themes" role="menu">
  87. <li role="none"><button role="menuitem" class="theme" id="light">Light <span class="default">(default)</span></button></li>
  88. <li role="none"><button role="menuitem" class="theme" id="rust">Rust</button></li>
  89. <li role="none"><button role="menuitem" class="theme" id="coal">Coal</button></li>
  90. <li role="none"><button role="menuitem" class="theme" id="navy">Navy</button></li>
  91. <li role="none"><button role="menuitem" class="theme" id="ayu">Ayu</button></li>
  92. </ul>
  93. <button id="search-toggle" class="icon-button" type="button" title="Search. (Shortkey: s)" aria-label="Toggle Searchbar" aria-expanded="false" aria-keyshortcuts="S" aria-controls="searchbar">
  94. <i class="fa fa-search"></i>
  95. </button>
  96. </div>
  97. <h1 class="menu-title">QCGPU User Guide</h1>
  98. <div class="right-buttons">
  99. <a href="print.html" title="Print this book" aria-label="Print this book">
  100. <i id="print-button" class="fa fa-print"></i>
  101. </a>
  102. </div>
  103. </div>
  104. </div>
  105. <div id="searchbar-outer" class="searchbar-outer">
  106. <input type="search" id="searchbar" name="searchbar" placeholder="Search this book ..." aria-controls="searchresults-outer" aria-describedby="searchresults-header">
  107. </div>
  108. <div id="searchresults-outer" class="searchresults-outer">
  109. <div class="searchresults-header" id="searchresults-header"></div>
  110. <ul id="searchresults">
  111. </ul>
  112. </div>
  113. <!-- Apply ARIA attributes after the sidebar and the sidebar toggle button are added to the DOM -->
  114. <script type="text/javascript">
  115. document.getElementById('sidebar-toggle').setAttribute('aria-expanded', sidebar === 'visible');
  116. document.getElementById('sidebar').setAttribute('aria-hidden', sidebar !== 'visible');
  117. Array.from(document.querySelectorAll('#sidebar a')).forEach(function(link) {
  118. link.setAttribute('tabIndex', sidebar === 'visible' ? 0 : -1);
  119. });
  120. </script>
  121. <div id="content" class="content">
  122. <main>
  123. <a class="header" href="algorithms/shor.html#shors-algorithm" id="shors-algorithm"><h1>Shors Algorithm</h1></a>
  124. <p>This algorithm finds the prime factors (\)u\) and \(v\)) of an odd, composite integer \(n\),
  125. that is not a prime power.</p>
  126. <a class="header" href="algorithms/shor.html#the-algorithm" id="the-algorithm"><h2>The Algorithm</h2></a>
  127. <p>(pseudo code)</p>
  128. <pre><code class="language-pseudo">Repeat
  129. Randomly choose \\(a \in \{ 2, \dots, n - 1 \}\\)
  130. Compute \\(d = gcd(a, n)\\)
  131. If \\(d \geq 2\\) then
  132. Return \\(u = d\\) and \\(v = n/d\\)
  133. Else // We know \\(a \in \mathbb{Z}^*_N\\)
  134. Let \\(r\\) be the order of \\(a\\) in \\(\mathbb{Z}^*_N\\) // Order finding algorithm
  135. If \\(r\\) is even then
  136. Compute \\(x = a^{r/2} - 1 (\mod n)\\)
  137. Compute \\(d = \gcd(x, n)\\)
  138. If \\(d \geq 2\\) then
  139. Return \\(u = d\\) and \\(v = n/d\\)
  140. Until a value is returned.
  141. </code></pre>
  142. <p>See <a href="https://cs.uwaterloo.ca/%7Ewatrous/LectureNotes.html">https://cs.uwaterloo.ca/~watrous/LectureNotes.html</a></p>
  143. <pre><pre class="playpen"><code class="language-rust">extern crate qcgpu;
  144. extern crate rand;
  145. use rand::{thread_rng, Rng};
  146. use qcgpu::{gcd, State};
  147. use qcgpu::gates::h;
  148. fn main() {
  149. let n = 15; // Number to factor
  150. println!(&quot;Factoring {}.&quot;, n);
  151. // Here we should check if the number is even, or if it is a power factor
  152. let mut rng = thread_rng();
  153. loop {
  154. let mut a = rng.gen_range(2, n); // Randomly choose \\(a \in \{ 2, \dots, n - 1 \}\\)
  155. let mut d = gcd(a, n);
  156. if d &gt;= 2 {
  157. // Found the factors; No Quantum needed
  158. println!(
  159. &quot;Factors are {} and {} (No quantum algorithm used)&quot;,
  160. d,
  161. n / d
  162. );
  163. break;
  164. } else {
  165. // We know \\(a \in \mathbb{Z}^*_N\\)
  166. let r = find_order(a, n);
  167. if r % 2 == 0 {
  168. let x = (a.pow((r as f32 / 2 as f32) as u32) - 1) % n;
  169. d = gcd(x, n);
  170. if d &gt;= 2 {
  171. println!(&quot;Factors are {} and {}&quot;, d, n / d);
  172. break;
  173. }
  174. } else {
  175. println!(&quot;Period is odd&quot;);
  176. }
  177. }
  178. }
  179. }
  180. </code></pre></pre>
  181. <a class="header" href="algorithms/shor.html#order-finding" id="order-finding"><h2>Order Finding</h2></a>
  182. <p>Given a positive integer \(n \geq 2\) and an element \(a \in \mathbb{Z}_n^* \),
  183. Find the order of \(a\) in \(\mathbb{Z}_n^*\).</p>
  184. <p>Order finding is the only quantum part in shors algorithm.</p>
  185. <p>\(\mathbb{Z}_n\) is the set of integers from \({0,\dots, n - 1}\) or \(\mathbb{z} \mod n\).
  186. The set \(\mathbb{Z}_n^* = {a \in \mathbb{Z}_n : \gcd(a,n) = 1}\)</p>
  187. <p>The set \( \mathbb{Z}_n^* \) forms a group wth the multiplication modulo \(n\) operation.
  188. Thus, for \( a \in \mathbb{Z}_n^* \) , \(\exists b \in \mathbb{Z}_n^*\) that uniquely satisfies</p>
  189. <p>\[
  190. ab \equiv 1 \mod n
  191. \]</p>
  192. <p>The order of a given element \(a \in \mathbb{Z}_n^*\) is the smallest positive integer \(r\) such that</p>
  193. <p>\[
  194. a^r \equiv 1 \mod n
  195. \]</p>
  196. <pre><pre class="playpen"><code class="language-rust">
  197. # #![allow(unused_variables)]
  198. #fn main() {
  199. fn find_order(a: i32, n: i32) -&gt; i32 {
  200. unimplemented!()
  201. }
  202. #}</code></pre></pre>
  203. </main>
  204. <nav class="nav-wrapper" aria-label="Page navigation">
  205. <!-- Mobile navigation buttons -->
  206. <a rel="prev" href="algorithms/grover.html" class="mobile-nav-chapters previous" title="Previous chapter" aria-label="Previous chapter" aria-keyshortcuts="Left">
  207. <i class="fa fa-angle-left"></i>
  208. </a>
  209. <a rel="next" href="algorithms/super-dense.html" class="mobile-nav-chapters next" title="Next chapter" aria-label="Next chapter" aria-keyshortcuts="Right">
  210. <i class="fa fa-angle-right"></i>
  211. </a>
  212. <div style="clear: both"></div>
  213. </nav>
  214. </div>
  215. </div>
  216. <nav class="nav-wide-wrapper" aria-label="Page navigation">
  217. <a href="algorithms/grover.html" class="nav-chapters previous" title="Previous chapter" aria-label="Previous chapter" aria-keyshortcuts="Left">
  218. <i class="fa fa-angle-left"></i>
  219. </a>
  220. <a href="algorithms/super-dense.html" class="nav-chapters next" title="Next chapter" aria-label="Next chapter" aria-keyshortcuts="Right">
  221. <i class="fa fa-angle-right"></i>
  222. </a>
  223. </nav>
  224. </div>
  225. <!-- Local fallback for Font Awesome -->
  226. <script>
  227. if (getComputedStyle(document.querySelector(".fa")).fontFamily !== "FontAwesome") {
  228. var link = document.createElement('link');
  229. link.rel = 'stylesheet';
  230. link.type = 'text/css';
  231. link.href = '_FontAwesome/css/font-awesome.css';
  232. document.head.insertBefore(link, document.head.firstChild)
  233. }
  234. </script>
  235. <script src="searchindex.js" type="text/javascript" charset="utf-8"></script>
  236. <script src="elasticlunr.min.js" type="text/javascript" charset="utf-8"></script>
  237. <script src="mark.min.js" type="text/javascript" charset="utf-8"></script>
  238. <script src="searcher.js" type="text/javascript" charset="utf-8"></script>
  239. <script src="highlight.js"></script>
  240. <script src="book.js"></script>
  241. <!-- Custom JS script -->
  242. </body>
  243. </html>