[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"doc-detail-84009-en":3,"doc-seo-84009-105":30,"detail-sidebar-cat-0-en-105":91},{"code":4,"msg":5,"data":6},0,"success",{"doc_id":7,"user_id":8,"nickname":9,"user_avatar":10,"doc_module":4,"category_id":11,"category_name":12,"doc_title":13,"doc_description":14,"doc_content":15,"file_id":16,"file_url":17,"file_type":18,"file_size":19,"view_count":20,"is_deleted":4,"is_public":21,"is_downloadable":21,"audit_status":21,"page_count":22,"language":23,"language_code":24,"site_id":25,"html_lang":24,"table_of_contents":26,"faqs":27,"seo_title":13,"seo_description":14,"update_tm":28,"read_time":29},84009,7971461740909,"Levi","https://ap-avatar.wpscdn.com/davatar_155a257f0dc6eb9ab79c44ca47cae57d",8,"Research & Report","Computing Smith Forms Modulo 2 of Sparse Matrices Faster Than Matrix Multiplication","Computing the Smith normal form over R = Z/p^2Z, the document studies sparse or structured matrices given only through a black-box action. For a prime p, every matrix has Smith factors with three possible multiplicities, and the work focuses on determining the counts of invariant factors equal to 1, p, and 0 efficiently. A randomized Monte Carlo algorithm is presented, using block-Wiedemann-style iterations and structured preconditioners, achieving expected R-operation complexity below the best dense baseline via fast multiplication exponents.","Computing Smith Forms Modulo 􀁰 2 of Sparse Matrices  \nFaster Than Matrix Multiplication  \nMark Giesbrecht  \nCheriton School of Computer Science, Faculty of Mathematics, University of Waterloo, Canada  \n[mwg@uwaterloo.ca](mwg@uwaterloo.ca)  \narXiv :2607 .05800v 1 [ cs . SC] 7 Jul 2026  \nAbstract. Let 􀀿 be a prime and R = Z/􀀿2Z the ring of integers modulo 􀀿2 . Any  \n􀀖 ∈ R􀀽×􀀽 is unimodularly equivalent to its Smith form  \n􀀨 = diag 􀀀, . . . , 1, 􀀿, . . . , 􀀿 , 0, . . . , 0􀀁 ∈ R􀀽×􀀽 ,  \n{z } | {z } | {z }  \n􀁁0 􀁁1 􀁁2  \ni.e., there exist 􀀪, 􀀫 ∈ R􀀽×􀀽 such that 􀀪 􀀖􀀫 = 􀀨, with det 􀀪, det 􀀫 ∈ R∗ (where R∗ is the set of units in R, elements not equivalent to 0 mod 􀀿) . Our goal in this paper is to determine 􀁁0, 􀁁 1, 􀁁 2 quickly when 􀀖 is sparse or structured. By “sparse”  \nwe mean 􀀖 is gi˜ ven by a black box such that for any 􀁅 ∈ R􀀽×1 we can compute  \n􀁅 ↦→ 􀀖􀁅 with 􀀤 (􀀽) operations in R, which captures having few nonzero elements or a multiplicative structure (e.g., Hankel or Toeplitz matrices) . We present a randomized algorithm which requires an expected number of  \n􀀤˜ 􀀐􀀽3−1/( 􀁬−1) 􀀑 operations in R  \nto compute the Smith form, where 􀁬 is the exponent of dense matrix multiplication. Usin˜ g standard cubic matrix multiplication (􀁬 = 3) our algorithm thus requires 􀀤 (􀀽2.5 ) operations in R, while using the current asymptotic˜ally fastest  \nmatrix multiplication, with 􀁬 \u003C 2. 371339, our algorithm requires 􀀤 (􀀽2.270786 ) operations in R. Our algorithm is probabilistic of the Monte Carlo type, meaning it fails on any invocation with controllably small probability. We employ iterative block-Wiedemann-style matrix techniques and structured preconditioners. To our  \nknowledge, this is the first algorith˜m to compute the modular Smith Normal Form  \nmodulo 􀀿2 requiring fewer than 􀀤 (􀀽􀁬 ) operations in R, i.e., faster than any dense algorithm.  \nKeywords: Smith normal form · sparse linear algebra · black-box matrices · block Wiedemann algorithm · local rings · matrix multiplication complexity  \n1 Introduction  \nThe Smith normal form (SNF), introduced by Smith (1861), diagonalizes matrices over Z  \nand has applications in diophantine analysis (Chou and Collins, 1982), integer programming (Hu, 1969), combinatorics (Stanley, 2016), determining the structure of abelian  \nTo appear, Computer Algebra in Scientific Computation (CASC) conference, August 31– September 4, 2026, Bath, UK  \n2 M. Giesbrecht  \ngroups (Newman, 1972) and class groups (Hafner and McCurley, 1989), in control system theory (Kailath, 1980), and especially computing simplicial homology (Dumas et al., 2003) . To deal with large instances of these problems we must exploit sparsity in the matrix (see, e.g., Giusti et al. (2016) in an application to neural data) . This remains challenging.  \nThe problem of computing the Smith form of a sparse matrix over a principal ideal ring presents several challenges. The mathematics goes back to Kaplansky (1949) . One approach is to simply compute the SNF over Z and then reduce the result modulo the power of the prime ideal. The algorithm of Eberly et al. (2007) for SNF of a sparse  \nmatrix over ˜Z could be used, which reduces the cost using standard matrix arithmetic  \nover Z to 􀀤 (􀀽2.66 ) bit operations, but does not seem directly amenable to asymptotic speedup. The best known algorithm f˜or any matrix over any Z/(􀀳), for any 􀀳 ∈ Z, is  \nby Storjohann (1996) and requires 􀀤 (􀀽􀁬 ) operations mod 􀀳, and is not sensitive to sparsity. The “valence” algorithm of Dumas et al. (2003) uses a hybrid of a “global”method over Z and a (dense) 􀀿 -adic elimination method for selected primes 􀀿 . Our method here may be useful as a subroutine in this approach.  \nWe will address a simple but important case here to attempt to break through this barrier: integer matrices modulo 􀀿2 for a prime 􀀿 . For matrices over R = Z/􀀿2Z, there are only three possible Smith factors, and the SNF has shape  \n􀀨 = diag 􀀀, . . . , 1, 􀀿, . . . , 􀀿 , 0, . . . , 0􀀁 = diag 􀀀 1􀁁0 , 􀀿 􀁁1 , 0􀁁2 􀀁","cbCaitZqYkgkDKWG","https://ap.wps.com/l/cbCaitZqYkgkDKWG","pdf",304672,4,1,19,"English","en",105,"# Introduction\n## Smith normal form background and applications\n## Prior work and remaining challenge\n## Goal and algorithmic approach\n## Matrix decomposition and Schur complement strategy\n## Computing multiplicities via efficient solves","[{\"question\":\"What problem does the paper solve for matrices modulo p^2?\",\"answer\":\"It computes the Smith normal form multiplicities for a matrix over R = Z/p^2Z, focusing on the counts of invariant factors equal to p and 0 without constructing expensive dense structures.\"},{\"question\":\"How is the input matrix represented when it is sparse or structured?\",\"answer\":\"The matrix is accessed as a black box: for any vector x in R^n, the algorithm can compute A x using a number of ring operations proportional to the available sparsity or structure.\"},{\"question\":\"What is the main algorithmic idea used to achieve faster complexity?\",\"answer\":\"The method splits off the unit part modulo p after fast rank-profile preconditioning, reducing the remaining task to solving a Schur-complement-like system whose residue is divisible by p and can be handled via structured, sparse block techniques and Hensel lifting.\"}]",1784191990,48,{"code":4,"msg":31,"data":32},"ok",{"site_id":25,"language":24,"slug":33,"title":13,"keywords":34,"description":14,"schema_data":35,"social_meta":86,"head_meta":88,"extra_data":90,"updated_unix":28},"computing-smith-forms-modulo-2-of-sparse-matrices-faster-than-matrix-multiplication","",{"@graph":36,"@context":85},[37,53,68],{"@type":38,"itemListElement":39},"BreadcrumbList",[40,44,48,51],{"item":41,"name":42,"@type":43,"position":21},"https://docshare.wps.com","Home","ListItem",{"item":45,"name":46,"@type":43,"position":47},"https://docshare.wps.com/document/","Document",2,{"item":49,"name":12,"@type":43,"position":50},"https://docshare.wps.com/document/research-report/",3,{"item":52,"name":13,"@type":43,"position":20},"https://docshare.wps.com/document/computing-smith-forms-modulo-2-of-sparse-matrices-faster-than-matrix-multiplication/84009/",{"url":52,"name":13,"@type":54,"author":55,"headline":13,"publisher":57,"fileFormat":60,"inLanguage":24,"description":14,"dateModified":61,"datePublished":62,"encodingFormat":60,"isAccessibleForFree":63,"interactionStatistic":64},"DigitalDocument",{"name":9,"@type":56},"Person",{"url":41,"name":58,"@type":59},"DocShare","Organization","application/pdf","2026-07-25","2026-07-16",true,{"@type":65,"interactionType":66,"userInteractionCount":20},"InteractionCounter",{"@type":67},"ViewAction",{"@type":69,"mainEntity":70},"FAQPage",[71,77,81],{"name":72,"@type":73,"acceptedAnswer":74},"What problem does the paper solve for matrices modulo p^2?","Question",{"text":75,"@type":76},"It computes the Smith normal form multiplicities for a matrix over R = Z/p^2Z, focusing on the counts of invariant factors equal to p and 0 without constructing expensive dense structures.","Answer",{"name":78,"@type":73,"acceptedAnswer":79},"How is the input matrix represented when it is sparse or structured?",{"text":80,"@type":76},"The matrix is accessed as a black box: for any vector x in R^n, the algorithm can compute A x using a number of ring operations proportional to the available sparsity or structure.",{"name":82,"@type":73,"acceptedAnswer":83},"What is the main algorithmic idea used to achieve faster complexity?",{"text":84,"@type":76},"The method splits off the unit part modulo p after fast rank-profile preconditioning, reducing the remaining task to solving a Schur-complement-like system whose residue is divisible by p and can be handled via structured, sparse block techniques and Hensel lifting.","https://schema.org",{"og:url":52,"og:type":87,"og:title":13,"og:site_name":58,"og:description":14},"article",{"robots":89,"canonical":52},"index,follow",{"doc_id":7,"site_id":25},{"code":4,"msg":5,"data":92},[93,97,101,105,110,115,120,123,128,131,135],{"id":21,"doc_module":4,"doc_module_name":46,"category_name":94,"show_sort_weight":95,"slug":96},"Story & Novel",90,"story-novel",{"id":47,"doc_module":4,"doc_module_name":46,"category_name":98,"show_sort_weight":99,"slug":100},"Literature",80,"literature",{"id":20,"doc_module":4,"doc_module_name":46,"category_name":102,"show_sort_weight":103,"slug":104},"Exam",70,"exam",{"id":106,"doc_module":4,"doc_module_name":46,"category_name":107,"show_sort_weight":108,"slug":109},5,"Comic",60,"comic",{"id":111,"doc_module":4,"doc_module_name":46,"category_name":112,"show_sort_weight":113,"slug":114},6,"Technology",50,"technology",{"id":116,"doc_module":4,"doc_module_name":46,"category_name":117,"show_sort_weight":118,"slug":119},7,"Healthcare",40,"healthcare",{"id":11,"doc_module":4,"doc_module_name":46,"category_name":12,"show_sort_weight":121,"slug":122},30,"research-report",{"id":124,"doc_module":4,"doc_module_name":46,"category_name":125,"show_sort_weight":126,"slug":127},9,"Religion & Spirituality",20,"religion-spirituality",{"id":126,"doc_module":4,"doc_module_name":46,"category_name":129,"show_sort_weight":126,"slug":130},"World Cup","world-cup",{"id":132,"doc_module":4,"doc_module_name":46,"category_name":133,"show_sort_weight":132,"slug":134},10,"Lifestyle","lifestyle",{"id":22,"doc_module":4,"doc_module_name":46,"category_name":136,"show_sort_weight":106,"slug":137},"General","general"]