2501.11896v2

# Systematic Abductive Reasoning via Diverse Relation Representations in Vector-symbolic Architecture **Authors**: Zhong-Hua Sun, Ru-Yuan Zhang, Zonglei Zhen, Da-Hui Wang, Yong-Jie Li, Xiaohong Wan, Hongzhi You > Corresponding author: Hongzhi You, email: hongzhi-you@uestc.edu.cnZhong-Hua Sun, Yong-Jie Li and Hongzhi You are with School of Life Science and Technology, University of Electronic Science and Technology of China (UESTC), Chengdu, China. Ru-Yuan Zhang is with Brain Health Institute, National Center for Mental Disorders, Shanghai Mental Health Center, Shanghai Jiao Tong University School of Medicine and School of Psychology, Shanghai, China. Zonglei Zhen and Xiaohong Wan are with State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University, Beijing␣ China. Da-Hui Wang is with School of Systems Science, Beijing Normal University, Beijing, China. Abstract In abstract visual reasoning, monolithic deep learning models suffer from limited interpretability and generalization, while existing neuro-symbolic approaches fall short in capturing the diversity and systematicity of attributes and relation representations. To address these challenges, we propose a Systematic Abductive Reasoning model with diverse relation representations (Rel-SAR) in Vector-symbolic Architecture (VSA) to solve Raven’s Progressive Matrices (RPM). To derive attribute representations with symbolic reasoning potential, we introduce not only various types of atomic vectors that represent numeric, periodic and logical semantics, but also the structured high-dimentional representation (SHDR) for the overall Grid component. For systematic reasoning, we propose novel numerical and logical relation functions and perform rule abduction and execution in a unified framework that integrates these relation representations. Experimental results demonstrate that Rel-SAR achieves significant improvement on RPM tasks and exhibits robust out-of-distribution generalization. Rel-SAR leverages the synergy between HD attribute representations and symbolic reasoning to achieve systematic abductive reasoning with both interpretable and computable semantics. Index Terms: Abstract visual reasoning, relation representation, vector-symbolic architecture. I Introduction Raven’s Progressive Matrices (RPM) are a family of psychological intelligence tests widely used for the assessment of abstract reasoning [1, 2]. From a cognitive psychology perspective, abstract visual reasoning in RPM tests involves constructing high-level representations from images and deriving potential relations from these representations [1, 3]. Endowing artificial intelligence with such capabilities is now regarded as a crucial step toward achieving human-level intelligence. However, many recent monolithic deep learning models, which do not explicitly separate perception and reasoning [4, 5, 6, 7, 8, 9], face inherent challenges, such as poor interpretability, limited robustness and generalization, and difficulties in module reuse [10]. Neuro-symbolic architecture, which combines neural visual perception with symbolic reasoning, offers a promising approach to overcoming these challenges and achieving human-level interpretability and generalization [11, 10, 12]. In neuro-symbolic architectures (NSA), Marcus argues that symbol-manipulation in cognition involves representing relations between variables [11]. For RPM tests, object attributes serve as the variables, while potential rules involve the relations. Nevertheless, due to incomplete attribute and relation representations, achieving systematic abduction and execution is still a critical challenge for NSA when performing RPM tests. From the perspective of attributes, recent models such as PrAE [10], the ALANS learner [13], and NVSA (neuro-vector-symbolic architecture) [12] construct attribute representations through neural perception frontends. Notably, the NVSA model achieves hierarchically structured VSA representations of image panels, capturing multiple objects with multiple attributes [12]. Regarding relation representations, PrAE and NVSA achieve abstract reasoning through probabilistic abduction and execution [10] and distributed vector-symbolic architecture (VSA) [12], respectively. Both models rely on predetermined multiple rule templates, each specialized for distinct individual RPM rules. To address the limitations in rule expressiveness, the ALANS learner utilizes learnable rule operators in the abstract algebraic structure, without manual definition for every rules [13]. Additionally, the ARLC model adopts a more expressive VSA-based rule template, operating in the rule parameter space [14]. Both models offer improved interpretabiltiy and generalizability. Despite their advances, previous models fall short in capturing the diversity and systematicity of attribute and relation representations. In contrast, human cognition demonstrates rich and flexible internal representations [15, 16], including arithmetic and logic, and rule-based reasoning systems in cognition are productive and systematic [17]. Therefore, the abstract visual reasoning performance of these models remains open to further improvement. Previous research indicates that Vector Symbolic Architecture (VSA), a form of high-dimensional (HD) distributed representation, possesses algebraic properties for mathematical operations and can also achieve structured symbolic representations of data [18, 19, 20]. In this work, to achieve comprehensive relation representations, we introduce various types of VSA-based atomic HD vectors with distinct semantic representations, including numeric values, periodic values, and logical values. Given that reasoning in RPM problems involves the overall attributes of multiple objects, we further introduce the structured HD representation (SHDR) for the nxn Grid. They serve as attribute representations necessary for abductive reasoning. Meanwhile, we propose numerical and logical relation functions as relation representations that take multiple HD attribute representations as input and define relations among them. Unlike rule templates designed for individual rules, the two proposed relation functions are specifically tailored to numerical and logical types, providing strong rule expressiveness. Here, we propose a Systematic Abductive Reasoning model with diverse relation representations (Rel-SAR) for solving RPM, inspired by the original NVSA model [12]. In the Rel-SAR model, visual attribute extraction and rule inference are implemented within a fully unified computational framework in the VSA machinery. The model comprises a neuro-vector frontend for perceiving object attributes of all raw images in RPM problems and a generic vector-symbolic backend for achieving symbolic reasoning. The perception frontend operates on scene-based SHDR of each image panel, which contains multiple objects, each with various attributes, and predicts HD attribute representations by VSA-based symbolic manipulations. The reasoning backend implements the core idea of systematic abductive reasoning: if the given attributes in an RPM adhere to a specific numerical or logical rule, then the relation representations of all attribute pairs can be defined using the corresponding relation functions with identical parameters. These diverse relation representations are involved in both rule abduction and execution phases, enhancing interpretability and improving the capacity for systematic abductive reasoning. II Related Work <details> <summary>x1.png Details</summary>  ### Visual Description ## Diagram: Context and Answer Panels for Geometric Pattern Matching ### Overview The image presents two 3x3 grid layouts labeled "Context Panels" (Panel A) and "Answer Panels" (Panel B). Each panel contains geometric shapes with positional relationships, testing pattern recognition and spatial reasoning. The task appears to involve matching context rules (Panel A) to corresponding answer configurations (Panel B). ### Components/Axes **Panel A (Context Panels):** - **Top Row:** 1. "Center" (black pentagon centered) 2. "Left-Right" (hexagon/pentagon split vertically) 3. "2x2Grid" (light gray pentagons in 2x2 arrangement) - **Middle Row:** 1. "Up-Down" (gray triangle above black triangle) 2. "3x3Grid" (dark gray triangles in 3x3 arrangement) 3. "Out-InCenter" (black pentagon inside light gray pentagon) - **Bottom Row:** 1. "Out-InGrid" (circles inside larger circle) **Panel B (Answer Panels):** - **Top Row:** 1. "Center" (black pentagon) 2. "Left-Right" (hexagon/pentagon split vertically) 3. "2x2Grid" (light gray pentagons in 2x2 arrangement) - **Middle Row:** 1. "Up-Down" (gray triangle above black triangle) 2. "3x3Grid" (dark gray triangles in 3x3 arrangement) 3. "Out-InCenter" (black pentagon inside light gray pentagon) - **Bottom Row:** 1. "Out-InGrid" (circles inside larger circle) ### Detailed Analysis **Panel A (Context Panels):** - **Top Row:** - "Center" shows a single black pentagon centered in the panel. - "Left-Right" displays a hexagon on the left and pentagon on the right, separated by a vertical line. - "2x2Grid" contains four light gray pentagons arranged in a 2x2 grid. - **Middle Row:** - "Up-Down" features a gray triangle above a black triangle. - "3x3Grid" has nine dark gray triangles arranged in a 3x3 grid. - "Out-InCenter" shows a black pentagon centered within a larger light gray pentagon. - **Bottom Row:** - "Out-InGrid" contains three circles (light gray, light blue, black) inside a larger circle. **Panel B (Answer Panels):** - **Top Row:** Mirrors Panel A's top row configurations. - **Middle Row:** Matches Panel A's middle row configurations. - **Bottom Row:** Includes "Out-InGrid" with circles inside a larger circle, matching Panel A's bottom row. ### Key Observations 1. **Symmetry in Labeling:** Both panels share identical labels, suggesting a direct mapping between context and answer options. 2. **Red Highlight:** The first answer panel in Panel B (top-left) is highlighted in red, indicating it may be the selected or correct answer. 3. **Shape Consistency:** Shapes (pentagons, triangles, circles) are reused across panels, emphasizing pattern recognition over novel geometry. 4. **Spatial Relationships:** Terms like "Out-InCenter" and "Out-InGrid" describe hierarchical positioning (e.g., shapes nested within larger shapes). ### Interpretation The diagram tests the ability to apply spatial rules from context panels to select matching answer configurations. The red-highlighted "Center" panel in Panel B suggests it is the correct match for the context provided in Panel A's top-left panel. This implies a task where users must identify which answer panel replicates the spatial arrangement or rule demonstrated in the corresponding context panel. The consistent use of geometric shapes across panels reinforces the focus on positional logic rather than abstract pattern recognition. </details> Figure 1: Illustrations for RAVEN dataset. (a) An example of RPM test from RAVEN [21] dataset. In an RPM test, there are 8 context panels and 8 candidate panels. Participants are required to identify the underlying rules governing various attributes within the context panels. Subsequently, participants use these rules to infer the attributes of the missing panel (represented by ”?”) and choose the most appropriate option (highlighted with a red box) from the answer panels. (b) The RAVEN dataset includes seven configurations: Center, 2x2Grid, 3x3Grid, Left-Right (L-R), Up-Down (U-D), Out-InCenter (O-IC) and Out-InGrid (O-IG) [21]. Four types of rules, i.e., Constant, Progression, Arithmetic, and Distribute Three, are applied to five attributes, i.e., Position, Number, Type, Size, and Color, in a row-wise manner. The I-RAVEN dataset [8] is a variant of RAVEN, where answer sets are generated using an attribute bisection tree. The Raven Progressive Matrices (RPM) is a widely used nonverbal intelligence test designed to assess abstract reasoning. To explore the limitations of current machine learning approaches in solving abstract reasoning tasks, two automatically generated RPM-based datasets—RAVEN [21] and I-RAVEN [8] —have been introduced (Figure 1). Early efforts on RPM primarily employed Relation Network (RN) [22] and their variants [4, 23, 7, 9] to extract relations between context panels. Concurrently, CoPINet [6], MLCL [24], and DCNet [25] integrate contrastive learning in their models. Approaches like MRNet [9] and DRNet [26] aimed to enhance perception capabilities, while SRAN [8] and PredRNet [27] abstract relations using stratified models and prediction errors, respectively. In addition, several methods have focused on scene decomposition and feature disentanglement [28, 29, 30]. Although these monolithic deep learning models achieve high accuracy, they often suffer from limited interpretability and systematic generalization capabilities. Another branch for solving RPM is based on neuro-symbolic architectures, which explicitly distinguish between perception and reasoning. PrAE [10] employs an object CNN to generate probabilistic scene representations and uses predetermined rule templates for probabilistic abduction and execution. Inspired by abstract algebra and representation theory, ALANS [13], which shares the same perception frontend as PrAE, transforms probabilistic scene distributions into matrix-based algebraic representations. The algebraic reasoning backend of ALANS induces potential rules through trainable operator matrices, eliminating the need for manual rule definitions. In abstract reasoning, Vector Symbolic Architectures (VSA) serve as a bridge between perception and reasoning modules by leveraging its structured distribution representations and algebraic properties. NVSA [12] projects each RPM panel into a high-dimensional vector using a trainable CNN and derives probability mass functions (PMFs) by querying an external codebook. Its reasoning backend embeds these PMFs into distributed VSA representations and performs rule abduction and execution using templates based on VSA algebraic operations. NVSA provides a differentiable and transparent implementation of probabilistic abductive reasoning by leveraging VSA representations and operators. However, its perception frontend requires searching a large external codebook, and its reasoning backend still relies on predetermined rule templates. In contrast, Learn-VRF [31], focuses on reasoning by learning VSA rule formulations, eliminating the need for predetermined templates. ARLC [14] further enhances reasoning by incorporating context augmentation and extending rule templates to accommodate more diverse rules. While ARLC and Learn-VRF implement systematic rule learning, they still struggle to process all RPM rules due to limitations in attribute representation. Recently, a class of methods known as relational bottlenecks has been proposed to enable efficient abstraction, but their capacity to handle complex relations remains uncertain [32, 33, 34, 35]. To address this limitation, Rel-SAR transforms perceptual inputs into high-dimensional attribute representations with symbolic reasoning potential and abducts both logical and numerical rules within a unified framework. III Preliminaries III-A VSA models utilized in this study VSAs are a class of computational models that utilize high-dimensional distributed representations [20]. VSA models used in this study are Holographic Reduced Representations (HRR) and its form in the frequency domain, referred to as Fourier Holographic Reduced Representations (FHRR) [36]. A random FHRR atomic vector, denoted as $\boldsymbol{\theta}:=\left\{\theta_{i}\right\}_{i=1}^{d}$ , is composed of elements $\theta_{i}$ that are independently sampled from a uniform distribution, specifically $\theta_{i}\sim\mathcal{U}(-\pi,\pi)$ [36]. The corresponding HRR atomic vector, $\boldsymbol{x}$ , is then obtained by applying the Inverse Fast Fourier Transform (IFFT) to $\boldsymbol{\theta}$ : $$ \boldsymbol{x}=\mathcal{F}^{-1}\left(e^{j\boldsymbol{\theta}}\right) \tag{1} $$ Here, $\mathcal{F}$ and $\mathcal{F}^{-1}(·)$ represent the Fast Fourier Transform (FFT) and Inverse FFT (IFFT), respectively. When the dimension $d$ is sufficiently large, these randomly generated vectors exhibit pseudo-orthogonality, making them suitable for representing distinct symbols or concepts. The similarity between any two vectors is a crucial metric for evaluating the distributed representations in VSAs. In FHRR and HRR, cosine similarity is employed to measure the similarity between two vectors [20]: $$ \begin{split}sim(\boldsymbol{\theta},\boldsymbol{\phi})&=\frac{1}{d}\sum_{i=1}% ^{d}{\cos\left(\theta_{i}-\varphi_{i}\right)}\\ sim(\boldsymbol{x},\boldsymbol{y})&=\frac{\boldsymbol{x}\cdot\boldsymbol{y}}{% \left|\boldsymbol{x}\right|\left|\boldsymbol{y}\right|}\end{split} \tag{2} $$ where $\boldsymbol{\theta}$ and $\boldsymbol{\phi}$ denote two FHRR vectors, and $\boldsymbol{x}$ and $\boldsymbol{y}$ two HRR vectors. The similarity $sim(·,·)$ ranges from -1 to +1, and above two similarity measures are equivalent. The pseudo-orthogonality refers to the case where the similarity $sim(·,·)≈ 0$ . III-B Basic operations and structured symbolic representations All computations within VSAs are composed of several basic vector algebraic operations, with the primary ones being binding ( $\circ$ ), bundling ( $+$ ) and unbinding ( $\oslash$ ) (Table I). The binding operation ( $\circ$ ) is employed to form a representation of an object that contains information about the context in which it was encountered [20]. The bundling operation ( $+$ ), also known as superposition, generates a composite high dimensional vector that combines several lower-level representations. In calculation, binding has a higher priority than bundling. The unbinding operation ( $\oslash$ ), which is the inverse of binding, extracts a constituent from the compound data structure. Binding and bundling are referred to as composition operations, while unbinding is considered a decomposition operation. All operations do not change the vector dimensionality. Through the combination of these operations, VSAs can effectively achieve structured symbolic representations [20]. For instance, consider a scene $\boldsymbol{s}$ in which a triangle $\boldsymbol{t}$ is positioned on the left $\boldsymbol{p}_{L}$ and a circle $\boldsymbol{c}$ on the right $\boldsymbol{p}_{R}$ . This scene can be represented as $\boldsymbol{s}=\boldsymbol{p}_{L}\circ\boldsymbol{t}+\boldsymbol{p}_{R}\circ% \boldsymbol{c}$ by the role-filler pair [37]. By applying the inverse vector of the left position $\boldsymbol{p}_{L}^{-1}$ to unbind $\boldsymbol{s}$ , we can retrieve an approximate vector representing the content at the left position, i.e., $\boldsymbol{p}_{L}^{-1}\otimes\boldsymbol{s}≈\boldsymbol{t}$ . Moreover, the triangle $\boldsymbol{t}$ can itself be a compositional scene, where attributes such as color and size are combined into a triangle scene in a similar manner. This decomposable, structure-sensitive, high-dimensional distributed representation has the potential to disentangle complex scenes while maintaining the advantages of traditional connectionist approaches [12]. III-C The fractional power encoding method In this study, the rules in RPM are primarily numerical. We introduce the VSA representation of numerical values using the fractional power encoding method (FPE-VSA) [18, 19]. Let $x∈\mathbb{R}$ be a real number and $X∈\mathbb{R}^{d}$ a randomly sampled base vector. The VSA representation $\boldsymbol{v}(x)∈\mathbb{R}^{d}$ for any value $x$ is obtained by repeatedly binding the base vector $X$ with itself $x$ times, as follows: $$ \boldsymbol{v}\left(x\right):=\left(X\right)^{\left(\circ x\right)} \tag{3} $$ The FPE method maps arbitrary real numbers to corresponding HD vector, and has the following properties: $$ \boldsymbol{v}\left(x_{1}+x_{2}\right)=\boldsymbol{v}\left(x_{1}\right)\circ% \boldsymbol{v}\left(x_{2}\right) \tag{4} $$ This demonstrates that addition $+$ in the real number domain can be represented by the binding operation $\circ$ in the vector domain. TABLE I: Basic oprations of FHRR and HRR. | Binding( $\boldsymbol{x}\circ\boldsymbol{y}$ ) Bundling( $\boldsymbol{x}+\boldsymbol{y}$ ) Inverse( $\boldsymbol{x}^{-1}$ ) | $\left(\boldsymbol{\theta}+\boldsymbol{\phi}\right)mod\,\,2\pi$ $angle\left(e^{j\boldsymbol{\theta}}+e^{j\boldsymbol{\phi}}\right)$ $\left(-\boldsymbol{\theta}\right)\,\,mod\,\,2\pi$ | $\mathcal{F}^{-1}\left(\mathcal{F}\left(\boldsymbol{x}\right)·\mathcal{F}% \left(\boldsymbol{y}\right)\right)$ $\boldsymbol{x}+\boldsymbol{y}$ $\mathcal{F}^{-1}\left(1/\mathcal{F}\left(\boldsymbol{x}\right)\right)$ | | --- | --- | --- | | Unbinding( $\boldsymbol{x}\oslash\boldsymbol{y}$ ) | $\boldsymbol{\theta^{-1}}\circ\boldsymbol{\phi}$ | $\boldsymbol{x^{-1}}\otimes\boldsymbol{y}$ | IV Methodology IV-A Atomic HD vectors with semantic representations <details> <summary>x2.png Details</summary>  ### Visual Description ## Diagram: Vector Types and Relational Operations ### Overview The image presents a technical diagram categorizing four vector types (Random, Numeric, Circular, Boolean) and two relational operation systems (Binary/Ternary Relations and Numerical/Logical Transformations). Visual elements include colored squares, directional arrows, and mathematical notation. ### Components/Axes **Top Row (Vector Types):** - **a. Random Vector**: Six colored squares (red, blue, green, orange, yellow, purple) with varying opacity, arranged asymmetrically. - **b. Numeric Vector**: Six blue squares aligned vertically with three dots above and below. - **c. Circular Vector**: Eight blue squares arranged in a perfect circle with a faint circular outline. - **d. Boolean Vector**: Two green squares connected by a light green circular outline. **Bottom Row (Relational Systems):** - **e. Binary/Ternary Relations**: Three gray squares labeled *v₁, v₂, v₃* with bidirectional arrows indicating relations. - **f. Numerical Transformation**: - Left: *v₁:N* → Middle: *R_Num/Lgc,N* → Right: *r* - Arrows labeled *OP₁:M* above input/output. - **g. Logical Transformation**: - Left: *v₁:N-1* → Middle: *R⁻¹_Num/Lgc,N* → Right: *v_N* - Arrows labeled *OP₁:M* and *r*. ### Detailed Analysis **Vector Types (a-d):** - **Random Vector (a)**: Colors distributed without pattern; red (top-left), blue (top-right), green (top-center), orange (center), yellow (bottom-center), purple (bottom-right). - **Numeric Vector (b)**: Uniform blue squares with implied numerical ordering via vertical alignment. - **Circular Vector (c)**: Equidistant blue squares forming a closed loop. - **Boolean Vector (d)**: Symmetrical green squares with connecting outline suggesting logical equivalence. **Relational Systems (e-g):** - **Binary/Ternary Relations (e)**: - *v₁* ↔ *v₂* (bidirectional arrow) - *v₂* ↔ *v₃* (bidirectional arrow) - *v₁* ↔ *v₃* (dashed arrow, implied ternary relation). - **Numerical Transformation (f)**: - Input: *v₁:N* (vector of length N) - Process: *R_Num/Lgc,N* (matrix operation) - Output: *r* (scalar result). - **Logical Transformation (g)**: - Input: *v₁:N-1* (vector of length N-1) - Process: *R⁻¹_Num/Lgc,N* (inverse matrix operation) - Output: *v_N* (single value). ### Key Observations 1. **Vector Symmetry**: - Circular Vector (c) shows perfect radial symmetry. - Boolean Vector (d) exhibits bilateral symmetry. 2. **Color Coding**: - Random Vector (a) uses distinct colors for each element. - Numeric/Circular Vectors (b,c) use uniform blue. 3. **Arrow Directionality**: - Binary relations (e) use bidirectional arrows. - Transformations (f,g) use unidirectional flow. 4. **Mathematical Notation**: - *R_Num/Lgc,N* and its inverse *R⁻¹_Num/Lgc,N* suggest linear algebra operations. - *OP₁:M* likely denotes a 1-to-M mapping function. ### Interpretation This diagram illustrates a framework for vector classification and transformation in computational systems: 1. **Vector Representation**: - Random vectors (a) model unstructured data. - Numeric (b) and Circular (c) vectors represent ordered/structured data. - Boolean (d) vectors encode binary states. 2. **Relational Operations**: - Binary/ternary relations (e) define interactions between vector elements. - Numerical transformation (f) reduces vectors to scalars via matrix operations. - Logical transformation (g) reconstructs vectors through inverse operations, implying error correction or data compression. 3. **System Integration**: - The *OP₁:M* notation suggests these operations handle variable-length inputs/outputs. - The inverse matrix *R⁻¹* in (g) implies reversibility, critical for lossless transformations. **Notable Anomalies**: - The Random Vector (a) lacks a clear pattern, contrasting with the structured Numeric/Circular Vectors. - The Boolean Vector (d) uses green exclusively, while other vectors use blue or mixed colors. This framework likely applies to machine learning (vector embeddings) or database systems (relational algebra), emphasizing structured data manipulation through linear transformations. </details> Figure 2: Atomic HD representations and relation functions. (a-d) The Rel-SAR model utilizes four types of atomic HD vectors. Random Vectors (RVs), sampled independently, are used to represent distinct and unrelated symbols or concepts. Numeric Vectors (NVs) are used to represent real numbers and support VSA-based addition-type arithmetic operations. Circular Vectors (CVs) represent periodic values and enable addition-type arithmetic operations with periodicity. Boolean Vectors (BVs), representing logical values of False and True, support VSA operations for logical reasoning. (e) In the RAVEN dataset, for a given attribute, the HD attribute representations in a row of three image panels involve binary or ternary relations. (f) Relation functions describe the numerical or logical relations between multiple HD vector representations $\boldsymbol{v}_{1:N}$ , where $N=2$ for binary and $N=3$ for ternary relations. These relations are governed by the operator powers $OP_{1:M}$ and the output $\boldsymbol{r}$ . (g) For a given relation defined by $OP_{1:M}$ and $\boldsymbol{r}$ , inverse relation functions infer the last HD vector representation $\boldsymbol{v}_{N}$ according to the first $N-1$ representations $\boldsymbol{v}_{1:N-1}$ . In neuro-vector-symbolic systems, atomic HD vector representations with meaningful semantics are essential for perception and reasoning. We introduce four types of atomic HD vectors used in our model (Figure 2): Random Vectors (RVs), Numeric Vectors (NVs), Circular Vectors (CVs), and Boolean Vectors (BVs). The definitions and properties of these vectors are universal within the VSA framework. IV-A 1 Random Vector RVs are sampled from specific distributions according to the VSA models, as mentioned in the preliminary section. Due to the absence of numerical or logical relations among RVs and their pseudo-orthogonality in the HD vector space (Figure 2 a), they are often used to represent symbols and concepts assumed to be independent and dissimilar. IV-A 2 Numeric Vector NVs, generated using the fractional power encoding (FPE-VSA, Equation 3) [18], are employed to represent real numbers (Figure 2 b). NVs $\boldsymbol{v}(r)∈\mathbb{R}^{d}$ can be used to perform addition-type arithmetic operations through the binding (Equation 4) [19]. IV-A 3 Circular Vector CVs are a special class of NVs used to represent periodic values (Figure 2 c). Given a base vector $P$ , where each phase of its elements $\rho_{i}$ is sampled from a discrete distribution ( e.g., for FHRR, $\rho_{i}\sim\mathcal{U}\left(2\pi j/L,∀ j∈\left\{1,·s,L\right\}\right)$ , with $L$ being an even number), CVs are defined as $\boldsymbol{p}(r):=(P)^{(\circ r)}$ . These CVs are pseudo-orthogonal to one another and exhibit periodicity with a period of $L$ [19]: $$ \boldsymbol{p}\left(r+L\right)=\boldsymbol{p}\left(r\right) \tag{5} $$ If $L$ is odd, the corresponding CVs with period $L$ can be obtained by selecting every other CV from those with period $2L$ . IV-A 4 Boolean Vector BVs are a specific type of CVs with a period of $L=2$ , used to represent Boolean values (Figure 2 d). Following a similar generation method as for CVs, we can generate vectors with a period of $L=2$ , $\boldsymbol{e}\left(0\right)$ and $\boldsymbol{e}\left(1\right)$ , to represent False and True, respectively. Basic logic operations using BVs are implemented as shown in Table II, where $\boldsymbol{a},\boldsymbol{b}∈\{\boldsymbol{e}(0),\boldsymbol{e}(1)\}$ represent arbitrary Boolean values. TABLE II: Logic operations implemented by BV $$ \lnot\boldsymbol{a}=\boldsymbol{a}\circ\boldsymbol{e}\left(1\right) \boldsymbol{a}\oplus\boldsymbol{b}=\boldsymbol{a}\circ\boldsymbol{b} \boldsymbol{a}\land\boldsymbol{b}=\boldsymbol{a}^{\left(\circ sim\left(% \boldsymbol{a},\boldsymbol{b}\right)\right)} \begin{aligned} \boldsymbol{a}\lor\boldsymbol{b}&=\left(\boldsymbol{a}\oplus% \boldsymbol{b}\right)\circ\left(\boldsymbol{a}\land\boldsymbol{b}\right)\\ \,\,&=\boldsymbol{a}\circ\boldsymbol{b}\circ\boldsymbol{a}^{\left(\circ sim% \left(\boldsymbol{a},\boldsymbol{b}\right)\right)}\\ \end{aligned} \tag{1} $$ IV-B Relation functions based on atomic HD representations The rules for abductive reasoning in RPM involve binary and ternary relations among the attributes of corresponding objects in each row of three panels (Figure 2 e and Figure 1 a), as well as numerical and logical relations. In this work, we design general relation functions based on VSA algebra, utilizing the aforementioned atomic vector representations, to be used for rule abductions. IV-B 1 Relation functions Relation functions, which describe the relations between multiple HD vector representations, are categorized into two types: numerical and logical. Among the atomic HD representations, Numeric Vectors (NVs) and Circular Vectors (CVs) are involved in numerical relations, while Boolean Vectors (BVs) are involved in logical relations. The numerical relation function, $R_{Num}$ , is defined as follows (Figure 2 f): $$ \boldsymbol{r}_{Num}=R_{Num}\left(\boldsymbol{v}_{1:N},OP_{1:M}\right)=\circ_{% i=1}^{N}\boldsymbol{v}_{i}^{\left(\circ op_{i}\right)} \tag{6} $$ where $N$ represents the arity of the relation function, and $\boldsymbol{v}_{1:N}:=\left\{\boldsymbol{v}_{i}\right\}_{i=1}^{N}$ denotes the input set of HD vector representations. $M$ is the number of operator powers and $OP_{1:M}:=\left\{op_{i}\right\}_{i=1}^{M}$ represents the operator powers, which can be considered as parameters of the relation function. The notation $\circ_{i=1}^{N}$ denotes the sequential binding operation applied to the $N$ HD vector representations. $\boldsymbol{r}_{Num}$ is the output HD representation. For the binary numerical relation function, $N=2$ and $M=2$ , while for the ternary numerical relation function, $N=3$ and $M=3$ . Based on the arithmetic properties of NVs and CVs, $R_{Num}$ can describe the additive relations of these two types of HD vector representations. The combination of $OP_{1:M}$ and $\boldsymbol{r}_{Num}$ determines the specific numerical relation in this vector-symbolic method. Similarly, the simplified logical relation function, $R_{Lgc}$ , is defined as follows (Figure 2 f): $$ \boldsymbol{r}_{Lgc}=R_{Lgc}\left(\boldsymbol{v}_{1:N},OP_{1:M}\right)=\left(% op_{1}\boldsymbol{v}_{1}\land op_{2}\boldsymbol{v}_{2}\right)\circ op_{3}% \boldsymbol{v}_{3} \tag{7} $$ where $\boldsymbol{v}_{1:N}:=\left\{\boldsymbol{v}_{i}\right\}_{i=1}^{N}∈\{% \boldsymbol{e}(0),\boldsymbol{e}(1)\}$ denotes the input set of BVs. The full version of logical relation function is described in Appendix A. Here, we consider only the ternary logical relation, so $N=3$ and $M=3$ . The parameter $OP_{1:M}:=\left\{op_{i}\right\}_{i=1}^{M}$ , where $op_{i}∈\{0,1\}$ determines whether to negate $\boldsymbol{v}_{i}$ , with negation ( $\lnot$ ) applied when $op_{i}=1$ and no negation applied when $op_{i}=0$ (see Appendix A). The symbol $\land$ denotes the AND operation, as shown in Table II. Based on the computational properties of BVs detailed in Table II, $R_{Lgc}$ can describe the ternary logical relations involved in RPM. The combination of the operator $OP_{1:M}$ and the output $\boldsymbol{r}_{Lgc}$ determines the specific logical relation in this vector-symbolic method. <details> <summary>x3.png Details</summary>  ### Visual Description ## Diagram: Scene Understanding and Object Positioning Model ### Overview The image presents a technical diagram illustrating a computational model for scene understanding and object positioning. It is divided into two main sections: - **Part a**: Focuses on entity attributes, scene representation, and spatial relationships. - **Part b**: Demonstrates a 3x3 grid system for object positioning and cyclic shifts. ### Components/Axes #### Part a 1. **Out-InGrid**: - A circular boundary labeled "Scene" containing three gray triangles. - A dashed red box highlights two triangles with a blue dashed arrow pointing to the "Entity" box. 2. **Entity Box**: - **Labels**: Type, Size, Color, Existence. - **RVS (Representation Values)**: - **Type**: Orange square (`v_j^type`), circle (`o_j`). - **Size**: Blue square (`v_j^size`), circle (`o_j`). - **Color**: Purple square (`v_j^color`), circle (`o_j`). - **Existence**: Green square (`v_j^exist`), circle (`o_j`). - **Summation**: `∑_Attribute` leads to a green square (`o_j`) and blue square (`o_j`). 3. **Scene Section**: - **SHDR for Image Panel**: A summation (`∑_Object`) leading to a black square (`S`). - **HD Representation**: Blue square (`o_j`). - **Role-Filler Binding**: Circle (`o_j`). 4. **Layout**: - **Position**: Labeled `p_j-1`, `p_j`, `p_j+1` with a blue square at `p_j`. #### Part b 1. **3x3 Grid**: - Labeled positions 1–9 in a 3x3 grid. - **Cyclic Shifts**: - `C1` (black pentagons) → `C2` (white circles). - `C2` → `C3` (gray pentagons). - `C3` → `C1` (implied by cyclic notation). 2. **Position and Existence**: - **Position**: Blue square (`p_j`). - **Existence**: Green square (`v_j^exist`). 3. **SHDR for 3x3 Grid**: - Summation (`∑_Object`) leading to a black square (`C^3x3`). ### Detailed Analysis - **Part a**: - The "Entity" box defines attributes (Type, Size, Color, Existence) with RVS values represented by colored squares (HD representation) and circles (role-filler binding). - The "Scene" section aggregates object data via SHDR, combining positional and attribute information. - **Part b**: - The 3x3 grid uses cyclic shifts (`p(1)`) to reposition objects (e.g., `C1` → `C2`). - Positions are mapped to grid cells (1–9), with existence confirmed by green squares. ### Key Observations 1. **Color-Coded Attributes**: - Each attribute (Type, Size, Color, Existence) has distinct colors (orange, blue, purple, green) for clarity. - Squares represent HD values, circles represent role-filler bindings. 2. **Cyclic Shifts**: - Objects (`C1`, `C2`, `C3`) are repositioned cyclically, suggesting a dynamic scene update mechanism. 3. **Summation Notation**: - `∑_Object` and `∑_Attribute` indicate aggregation operations for scene and entity representations. ### Interpretation This diagram models a scene understanding system where: - **Entities** are defined by attributes (Type, Size, Color, Existence) with HD and role-filler representations. - **Scene Composition**: Objects are positioned in a 3x3 grid, with cyclic shifts enabling dynamic updates. - **SHDR**: Serves as a computational tool to aggregate object data into a unified scene representation (`S` or `C^3x3`). The model emphasizes spatial relationships (`p_j`, `p_j±1`) and attribute-based entity classification, likely for applications in computer vision or robotics. The cyclic shifts suggest a mechanism for tracking object movements or transformations over time. </details> Figure 3: Structured HD representations (SHDR) for the image panel and the nxn Grid. (a) SHDR for the image panel (Equation 10). Taking the Out-InGrid configuration as an example, an image panel contains multiple objects, each with four entity attributes: type, size, color, and existence. Through the first layer of role-filler binding, these attributes are combined to form a SHDR for each object. Additionally, at the layout level, each object has a position attribute. By applying a second layer of role-filler binding, the SHDR for the entire image panel is constructed. (b) SHDR for the nxn Grid (Equation 12). Taking the 3x3Grid configuration as an example, the position vectors for all objects are represented using circular vectors (CVs) with a period of $3× 3=9$ . The SHDR $\mathcal{C}^{3\mathrm{x}3}$ for this 3x3Grid is obtained by performing role-filler binding between the corresponding position vectors $\boldsymbol{p}_{j}$ and existence vectors $\boldsymbol{v}_{j}$ . Due to the periodic nature of the position vectors, as all objects shift positions cyclically, the SHDR undergoes a binding operation with the position vectors corresponding to the magnitude of the shift. IV-B 2 Inverse relation functions Rule execution in RPM requires inferring the third attribute value based on the first two attribute values in a row of panels, given a known relation. It represents an inverse problem of rule abduction. In the vector-symbolic method, given the operator power ${OP_{1:M}}$ and the output $\boldsymbol{r}$ , the last vector representation $\boldsymbol{v}_{N}$ can be inferred from the first $N-1$ inputs $\boldsymbol{v}_{1:N-1}$ using the inverse of the relation functions (Figure 2 g). According to Equation 6, the inverse numerical relation function is defined as follows: $$ \displaystyle\boldsymbol{v}_{N} \displaystyle=R_{Num}^{-1}\left(\boldsymbol{v}_{1:N-1},OP_{1:M},\boldsymbol{r}\right) \displaystyle=\left(\circ_{i=1}^{N-1}\boldsymbol{v}_{i}^{\left(\circ\left(-op_% {i}/op_{M}\right)\right)}\right)\circ\boldsymbol{r}^{\left(\circ\left(-1/op_{M% }\right)\right)} \tag{8} $$ Similarly, according to Equation 7, the inverse logical relation function is defined as follows: $$ \displaystyle\boldsymbol{v}_{N} \displaystyle=R_{Lgc}^{-1}\left(\boldsymbol{v}_{1:N-1},OP_{1:M},\boldsymbol{r}\right) \displaystyle=op_{3}\left(op_{1}\boldsymbol{v}_{1}\land op_{2}\boldsymbol{v}_{% 2}\right) \tag{9} $$ . IV-C Structured high-dimensional representation and its attribution decomposition VSA can create structured symbolic representations using atomic HD vector representations and decouple them directly from these structures through algebraic operations [12]. This subsection presents the process of constructing a structured HD representation (SHDR) for an image panel and its decomposition to retrieve individual attribute representations. Additionally, an SHDR for the nxn Grid ( $n=2,3$ ) at the component level is also introduced. IV-C 1 SHDR for the image panel In RAVEN dataset, each image panel consists of objects, with each object characterized by multiple attributes. Consequently, the structured HD representation (SHDR) for each image panel can be obtained through two layers of role-filler bindings (Figure 3 a). First, the bundling operation is used to construct an SHDR for each object at the entity level by combining its attributes. Then, another bundling operation aggregates these object-level representations to construct a SHDR of the image panel at the scene level. Therefore, each image panel $\mathcal{X}∈\mathbb{R}^{r× r}$ , with a resolution $r× r$ , can be represented by an SHDR $\mathcal{S}∈\mathbb{R}^{d}$ as follows: $$ \displaystyle\mathcal{S} \displaystyle=\sum_{j=1}^{N_{pos}}{\boldsymbol{p}_{j}\circ\boldsymbol{O}_{j}} \displaystyle=\sum_{j=1}^{N_{pos}}{\boldsymbol{p}_{j}\circ\left(\sum_{attr\in ATTR% }{\boldsymbol{k}_{attr}\circ\boldsymbol{v}_{j}^{attr}}\right)} \tag{10} $$ Here, $\boldsymbol{O}_{j}$ represents the SHDR of the $j$ th object with different attributes at the entity level, incorporating attributes such as type, size, color, and existence. The attribute set is $ATTR=\left\{type,size,color,exist\right\}$ . At the entity level, the key vector $\boldsymbol{k}_{attr}$ denotes the class of a specific attribute $attr∈ ATTR$ , while the value vector $\boldsymbol{v}_{j}^{attr}$ indicates the attribute’s value at the position $j$ . At the scene level, the position vector $\boldsymbol{p}_{j}$ specifies the location of the $j$ -th object. IV-C 2 Representation decomposition Given an estimated SHDR $\hat{\mathcal{S}}∈\mathbb{R}^{d}$ of an image panel, all SHDRs of objects $\hat{\boldsymbol{O}}_{j}$ at the entity level, along with the corresponding attribute representations $\hat{\boldsymbol{v}}_{j}^{attr}$ , can be derived through a series of unbinding operations [20]. The decomposition process is shown as follows: $$ \begin{cases}\hat{\boldsymbol{O}}_{j}=\boldsymbol{p}_{j}\oslash\hat{\mathcal{S% }}=\boldsymbol{p}_{l}^{-1}\circ\hat{\mathcal{S}}\\ \hat{\boldsymbol{v}}_{j}^{attr}=\boldsymbol{k}_{attr}\oslash\hat{\boldsymbol{O% }}_{j}=\boldsymbol{k}_{attr}^{-1}\circ\hat{\boldsymbol{O}}_{j}\\ \end{cases} \tag{11} $$ It is important to note that due to inaccuracies in the estimated SHDR $\hat{\mathcal{S}}$ and the noise introduced by the unbinding operation, the estimated attribute representations $\hat{\boldsymbol{v}}_{j}^{attr}$ may not fully match the original $\boldsymbol{v}_{j}^{attr}$ used in Equation 10. IV-C 3 SHDR for the nxn Grid component In the RAVEN dataset, three figure configurations— 2x2Grid, 3x3Grid, and Out-InGrid —include components where objects are arranged in an nxn grid pattern at the layout level [21]. Since the positions in the nxn Grid involve component-level rule reasoning, the SHDR for the nxn Grid component ( $n=2,3$ ), focusing only on positions and object existence, is introduced as follows (Figure 3 b): $$ \mathcal{C}^{\mathrm{nxn}}=\sum_{j=1}^{n\times n}{\boldsymbol{p}_{j}\circ% \boldsymbol{v}_{j}^{exist}} \tag{12} $$ IV-D Rules from the perspective of relation functions TABLE III: Attribute representations and the relation functions involved in rule abductions {tblr} hlines = , vlines = , colspec=ccccccc, cells=mode=text, cell13 = c = 2halign = c, cell15 = c = 2halign = c, cell21 = r = 3valign = m, cell51 = r = 5valign = m, cell62 = r = 4valign = m, cell23 = r = 4valign = m, cell63 = r = 3valign = m, cell24 = r = 4valign = m, cell64 = r = 3valign = m, cell26 = r = 7valign = m, cell27 = r = 2valign = m, cell47 = r = 2valign = m, cell67 = r = 2valign = m Level & Attributes HD representations Rules in RAVEN and their types Relation functions Entity Type $\begin{array}[]{c}\boldsymbol{v}^{type}\\ \boldsymbol{v}^{size}\\ \boldsymbol{v}^{color}\\ \boldsymbol{v}^{num}\\ \end{array}$ Atomic (NVs) Constant Numerical rules Binary + Numerical Size Progression Color Distribute Three Ternary + Numerical Layout Number Arithmetic Position $\mathcal{C}^{\mathrm{nxn}}=\sum_{j=1}^{n× n}{\boldsymbol{p}_{j}\circ% \boldsymbol{v}_{j}^{exist}}$ $\begin{array}[]{c}\mathrm{SHDR}\\ \boldsymbol{p}_{j}:\mathrm{CVs}\\ \boldsymbol{v}_{j}^{exist}:\mathrm{RVs}\\ \end{array}$ Constant Binary + Numerical Progression Distribute Three Ternary + Numerical $\boldsymbol{v}^{exist}$ Atomic (BVs) Arithmetic Logical rules Ternary + logical TABLE IV: Rules and Corresponding combinations of $OP_{1:M}$ and $r_{Num}$ in Relation functions {tblr} hlines = , vlines = , colspec=cccccccc, cells=mode=text, cell11 = r = 2, c = 2halign = c, valign = m, cell17 = r = 2valign = m, cell18 = r = 2valign = m, cell31 = r = 6valign = m, cell32 = c = 2halign = l, cell42 = r = 2valign = m,halign = l, cell62 = r = 2valign = m,halign = l, cell82 = c = 2halign = c,halign = l, cell91 = r = 2valign = m, cell92 = r = 2valign = m,halign = l Rules in RAVEN & Ternary $op_{1}$ $op_{2}$ $op_{3}$ $r$ Examples (Rule $→$ Relation function) Binary $op_{1}$ $op_{2}$ Numerical rules Constant 0 $-1$ $+1$ 0 $v_{1}=v_{2}→\boldsymbol{v}\left(0\right)=\boldsymbol{v}_{1}^{\left(% \circ\left(-1\right)\right)}\circ\boldsymbol{v}_{2}^{\left(\circ 1\right)}$ Progression $+$ 0 $-1$ $+1$ $+1$ , $+2$ $v_{1}+1=v_{2}→\boldsymbol{v}\left(+1\right)=\boldsymbol{v}_{1}^{% \left(\circ\left(-1\right)\right)}\circ\boldsymbol{v}_{2}^{\left(\circ 1\right)}$ $-$ 0 $-1$ $+1$ $-1$ , $-2$ $v_{1}-2=v_{2}→\boldsymbol{v}\left(-2\right)=\boldsymbol{v}_{1}^{% \left(\circ\left(-1\right)\right)}\circ\boldsymbol{v}_{2}^{\left(\circ 1\right)}$ Arithmetic $+$ $-1$ $-1$ $+1$ $0 0$ $v_{1}+v_{2}=v_{3}→\boldsymbol{v}\left(0\right)=\boldsymbol{v}_{1}^{% \left(\circ\left(-1\right)\right)}\circ\boldsymbol{v}_{2}^{\left(\circ(-1)% \right)}\circ\boldsymbol{v}_{3}^{\left(\circ 1\right)}$ $-$ $-1$ $+1$ $+1$ $0 0$ $v_{1}-v_{2}=v_{3}→\boldsymbol{v}\left(0\right)=\boldsymbol{v}_{1}^{% \left(\circ\left(-1\right)\right)}\circ\boldsymbol{v}_{2}^{\left(\circ 1\right% )}\circ\boldsymbol{v}_{3}^{\left(\circ 1\right)}$ Distribute Three $+1$ $+1$ $+1$ Any $v_{1}+v_{2}+v_{3}=Any→\boldsymbol{Any}=\boldsymbol{v}_{1}^{\left(% \circ 1\right)}\circ\boldsymbol{v}_{2}^{\left(\circ 1\right)}\circ\boldsymbol{% v}_{3}^{\left(\circ 1\right)}$ Logical rules Arithmetic $+$ $+1$ $+1$ $+1$ 0 $\boldsymbol{e}\left(0\right)=\left(\boldsymbol{e}\left(1\right)\circ% \boldsymbol{v}_{1}\land\boldsymbol{e}\left(1\right)\circ\boldsymbol{v}_{2}% \right)\circ\left(\boldsymbol{e}\left(1\right)\circ\boldsymbol{v}_{3}\right)$ $-$ $0 0$ $+1$ $0 0$ 0 $\boldsymbol{e}\left(0\right)=\left(\boldsymbol{e}\left(0\right)\circ% \boldsymbol{v}_{1}\land\boldsymbol{e}\left(1\right)\circ\boldsymbol{v}_{2}% \right)\circ\left(\boldsymbol{e}\left(0\right)\circ\boldsymbol{v}_{3}\right)$ The RAVEN dataset contains $4$ rules— Constant, Progression, Arithmetic, and Distribute Three —which operate on $5$ rule-governing attributes [21]. These $5$ attributes include $3$ entity-level attributes: Type, Size, and Color, as well as $2$ layout-level attributes: Number and Position. In this study, the HD representations of these attribute values during rule reasoning and relations between rules and relation functions are shown in Table III. For the attributes Type, Size, Color, and Number, the four involved rules follow additive arithmetic operations, meaning the attribute values $\boldsymbol{v}^{attr}$ ( $attr∈\left\{type,size,color,number\right\}$ ) are represented using Numeric Vectors (NVs). Therefore, these rules can be defined using the numerical relation function (Equation 6): Constant and Progression correspond to binary relation functions, while Arithmetic and Distribute Three correspond to ternary relation functions. Each rule is associated with specific combinations of $OP_{1:M}$ and $\boldsymbol{r}_{Num}$ , and corresponding details are shown in Table IV. For the attribute Position, the rules Constant, Progression, and Distribute Three primarily refer to an nxn Grid with multiple objects, which, in an overall sense, follow additive arithmetic operations. Therefore, we use the SHDR $\mathcal{C}^{\mathrm{nxn}}$ for the nxn Grid (Equation 12) to represent the attributes required by these three rules. Since Progression involves a cyclic left or right shift of all objects (Figure 3 b), the position vectors $\boldsymbol{p}_{j}$ in set $\mathcal{C}^{\mathrm{nxn}}$ during rule reasoning are represented by Circular Vectors (CVs). The object existence vectors $\boldsymbol{v}_{j}^{exist}$ are represented by Random Vectors (RVs). These three rules can also be described using numerical relation functions (Equation 6). Take the rule Progression (+1) as an example, where the positions of objects undergo a cyclic right shift. In Figure 3 b, the SHDRs of the 3x3 Grid across a row of three panels exhibit two numerical relations: $\mathcal{C}_{1}\circ\boldsymbol{p}\left(1\right)=\mathcal{C}_{2}$ and $\mathcal{C}_{2}\circ\boldsymbol{p}\left(1\right)=\mathcal{C}_{3}$ , which can be defined using a binary numerical relation function. In addition, the rule Arithmetic on the attribute Position is belong to the logical rule [21]. The attribute values $\boldsymbol{v}_{j}^{exist}$ are also represented using Boolean Vectors (BVs) that can be operated as shown in Table II. Therefore, the rule Arithmetic for Position corresponds to the ternary logical relation function (Equation 7), and corresponding details about the combinations of $OP_{1:M}$ and $\boldsymbol{r}_{Num}$ in the relation function are shown in Table IV. IV-E The Systematic Abductive Reasoning model In this section, we present the Systematic Abductive Reasoning model with diverse relation representations (Rel-SAR), inspired by the NVSA [12]. An overview of Rel-SAR is depicted in Figure 4 a. Similar to previous neuro-symbolic models for abstract visual reasoning, Rel-SAR combines a neural visual perception frontend with a symbolic reasoning backend, both utilizing VSA representations with meaningful semantics to facilitate systematic reasoning. The perception frontend employs a neural network to extract the SHDR $\mathcal{S}$ of each image panel $\mathcal{X}$ in the RPM and achieves feature disentanglement from the SHDR using representation decomposition to obtain the HD representations of attributes ( $\boldsymbol{v}$ , $\boldsymbol{p}$ and $\mathcal{C}$ : Table III) required for reasoning in the backend. The reasoning backend consists of three main modules: the rule abduction module, the rule execution module, and the answer selection module. The rule abduction module extracts the corresponding rules ( $OP_{1:M}$ and $\boldsymbol{r}$ : Table IV) for each attribute representation according to appropriate relation function (Equation 6 and 7, Table III). Subsequently, the rule execution module uses these rules to predict the representations of the missing panel’s attributes according to corresponding inverse relation functions (Equation 8 and 9). Finally, the answer selection module compares the predicted attribute representations of the missing panel with the available options in the answer panels and selects the answer. <details> <summary>x4.png Details</summary>  ### Visual Description ## System Architecture Diagram: Perception, Reasoning, and Rule Execution Framework ### Overview The image presents a three-part technical diagram illustrating a cognitive system architecture. Part A shows the system's frontend-backend structure, Part B details attribute processing components, and Part C demonstrates rule-based reasoning mechanisms. The diagrams use standardized flowchart notation with color-coded components and directional arrows indicating information flow. ### Components/Axes **Part A (Frontend-Backend Structure):** - **Perception Frontend:** - Input: Contexts (X₁₁, X₁₂, X₃₂), Candidates (X₁, X₈) - Processing: SHDR (Spatial Hierarchical Density Representation), Feature disentanglement - Output: Object attributes (Type, Size, Color, Position) - **Reasoning Backend:** - Components: Rule abduction, Rule execution - Output: Answer selection (X̂_y) **Part B (Attribute Processing):** - Input: SHDR (Ŝ) - Components: - Frontend Codebook (C_Front) - Backend Codebook (C_Back) - Query & Attention mechanisms - Output: High-dimensional (HD) attribute representations **Part C (Rule-Based Reasoning):** - **Rule Abduction:** - Input: Numerical rules (f_φ), Attribute sets (V², V³) - Components: - 2-ary Relation (R_Num,2) - 3-ary Relation (R_Num,3) - Output: Unnormalized rule probabilities (s_Num², s_Num³) - **Rule Execution:** - Input: Logical rules (f_Lgc) - Components: - 3-ary Relation (R_Lgc,3) - Rule execution module (R⁻¹_Num, R⁻¹_Lgc) - Output: Final attribute values (v₃₁, v₃₂, v₃₃) ### Detailed Analysis **Part A Flow:** 1. Contextual information (X) → SHDR processing → Feature disentanglement 2. Disentangled features → Attribute extraction (Type, Size, Color, Position) 3. Attributes → Rule abduction → Rule execution → Final answer selection **Part B Flow:** 1. SHDR (Ŝ) → Frontend Codebook (C_Front) → Query generation 2. Query → Attention mechanism → Backend Codebook (C_Back) 3. HD attribute representations generated through this pipeline **Part C Mechanisms:** - **Numerical Rules:** - Attribute sets (V², V³) → 2-ary/3-ary relations - Probability calculation: argmax(s_Num², s_Num³) - **Logical Rules:** - Attribute sets (V³) → 3-ary relations - Probability calculation: s_Lgc³ - Both pathways converge through rule execution modules ### Key Observations 1. **Modular Design:** Clear separation between perception (frontend) and reasoning (backend) components 2. **Attribute Hierarchy:** Attributes progress from basic features (Type, Size) to complex representations (HD attributes) 3. **Rule Integration:** Dual pathways for numerical and logical rule processing 4. **Temporal Flow:** Information flows sequentially from perception through reasoning to final output 5. **Color Coding:** - Blue: Core processing modules - Yellow: Frontend components - Red: Backend components - Orange: Rule-related elements ### Interpretation This architecture demonstrates a hybrid cognitive system combining: 1. **Perceptual Processing:** Initial context analysis and feature extraction 2. **Attribute Representation:** Multi-level feature disentanglement and HD representation 3. **Rule-Based Reasoning:** Dual-path integration of numerical and logical inference 4. **Decision Making:** Final answer selection through rule execution The system appears designed for complex pattern recognition and decision-making tasks, with explicit mechanisms for: - Contextual awareness (multiple context inputs) - Feature abstraction (from raw data to HD attributes) - Multi-modal rule application (numerical + logical) - Probabilistic inference (unnormalized rule probabilities) Notable design choices include the use of both 2-ary and 3-ary relational processing, suggesting capability to handle both binary and ternary relationships in reasoning tasks. The attention mechanism in Part B implies dynamic feature weighting during processing. </details> Figure 4: The Systematic Abductive Reasoning model with diverse relation representations (Rel-SAR). (a) Overall architecture of Rel-SAR model. Our model consists of a visual perception frontend, which processes object attributes for $8$ context and $8$ candidate image panels in a RPM test, and a reasoning backend that performs symbolic arithmetic and logical reasoning. The perception frontend utilizes a neural network, $f_{\theta}$ , to obtain the SHDR of each image panel, and then perceives attributes in the form of HD representations required by the downstream reasoning. In the reasoning backend, the rule abduction module extracts rules for each attribute representation using relation functions. The rule execution module then predicts the missing panel’s attribute representations based on inverse relation functions. Finally, the answer selection module compares the predicted attributes of the missing panel with those in the candidate panels and selects the option with the highest similarity. (b) Given the predicted SHDR for each panel, the SHDR of all objects and their corresponding HD attribute representations can be obtained via representation decomposition. Subsequently, the estimated HD attribute representations are refined in two steps: querying the frontend codebook and applying attention based on the backend codebook. This process produces HD attribute representations suitable for backend reasoning, including attributes such as type, size, color, number, and position. (c) In the rule abduction module, the rule learners $f_{\phi}^{Num}$ and $f_{\varphi}^{Lgc}$ predict the operator powers $\widehat{OP}_{1:M}$ for numerical and logical relation functions based on attributes in the context panels. These predicted $\widehat{OP}_{1:M}$ ensure that all binary or ternary relation input pairs ( $\mathbb{V}^{N},N=2,3$ ) produce the same output $\boldsymbol{\hat{r}}$ when processed through their respective relation functions. Therefore, the rule defined by $\widehat{OP}_{1:M}$ and $\boldsymbol{\hat{r}}$ with the highest overall $\boldsymbol{\hat{r}}$ similarity, also viewed as unnormalized probability, is considered the underlying rule. The rule execution module then predicts the attributes of the missing panel using inverse relation functions with the estimated rules. IV-E 1 Perception frontend The perception frontend operates independently on each of the 16 image panels to extract the HD representations of attributes required for abductive reasoning (Figure 4 a and Figure 4 b). For a given image panel $\mathcal{X}^{ind}∈\mathbb{R}^{r× r}$ , where $ind∈\left\{\left(1,1\right),\left(1,2\right),·s,\left(3,2\right)\right\}$ for 8 contexts and $ind∈\left\{1,2,·s,8\right\}$ for 8 candidates, the frontend uses a trainable neural network (ResNet-50) to map the image panel to its estimated SHDR $\hat{\mathcal{S}}^{ind}∈\mathbb{R}^{d}$ : $f_{\theta}:\mathcal{X}→\hat{\mathcal{S}}$ , where $\theta$ represents the trainable parameters of the network. Theoretically, the expected SHDR $\mathcal{S}$ for each panel should be organized from the corresponding attribute representations as described by Equation 10. Therefore, the learning objective of $f_{\theta}$ is to minimize the difference between its output $\hat{\mathcal{S}}$ and the theoretical SHDR $\mathcal{S}$ , formulated as: $$ \underset{\theta}{\min}\left\|f_{\theta}\left(\mathcal{X};\theta\right)-% \mathcal{S}\right\| \tag{13} $$ Subsequently, the estimated SHDR $\hat{\mathcal{S}}^{ind}$ for each panel undergoes representation decomposition (Figure 4 b), as described in Equation 11, to obtain the estimated HD attribute representations for each object, including Type ( $\hat{\boldsymbol{v}}^{type}_{j}$ ), size ( $\hat{\boldsymbol{v}}^{size}_{j}$ ), color ( $\hat{\boldsymbol{v}}^{color}_{j}$ ), and existence ( $\hat{\boldsymbol{v}}^{exist}_{j}$ ), where $j$ denotes the position index of the corresponding object. The HD attribute representations are expected to be selected from a set of frontend codebooks for the available attributes of interested in the RAVEN dataset (Figure 4 b). These frontend codebooks include $C_{Num}^{Front}:=\left\{\boldsymbol{v}(r)\right\}_{r=0}^{9}\cup\left\{% \boldsymbol{v}_{null}\right\}$ and $C_{Lgc}^{Front}:=\left\{\boldsymbol{e}(r)\right\}_{r=0}^{1}$ , which represent the numerical value and logic, respectively. $\boldsymbol{v}_{null}$ represents the null attribute representation when there is no object. To improve the neural network’s performance in encoding the SHDR of an image panel, all hypervectors in these frontend codebooks are randomly and independently generated as RVs, rather than using NVs, CVs, or BVs. However, the estimated HD attribute representations for each object, $\hat{\boldsymbol{v}}^{attr}_{j}$ ( $attr∈\{type,size,color,exist\}$ ), cannot be directly applied to the reasoning backend. First, these representations contain noise introduced by the bundling operation in the form SHDR. Second, as they are expected to be derived from the frontend codebooks of RVs, there are no intrinsic arithmetic or logical relations between the $\hat{\boldsymbol{v}}^{attr}_{j}$ s, which hinders effective reasoning. To address these issues, we adopt an approach similar to the attention mechanism [38] to obtain the HD attribute representations suitable for the reasoning backend (Figure 4 b). In the query stage, we use the estimated HD attribute representations ( $\hat{\boldsymbol{v}}^{attr}_{j}$ ) as query vectors to compute their similarity with all possible vectors for the corresponding attributes in the frontend codebooks. In the attention stage, these similarity scores are then used as attention weights to perform a weighted summation of the corresponding vectors from the backend codebooks, in which all hypervectors are generated according to their attribute type as shown in Table III. The backend codebooks consist of $C_{Num}^{Back}$ (NVs), $C_{Lgc,BV}^{Back}$ (BVs), $C_{Lgc,RV}^{Back}$ (RVs), and $C_{Pos,\mathrm{nxn}}^{Back}:=\left\{\boldsymbol{p}_{r}\right\}_{r=1}^{n^{2}}$ (CVs for positions in nxn Grid). The updated HD attribute representations obtained after the weighted summation can be utilized in the reasoning backend. The details are provided below. The query stage: For each attribute $attr∈\left\{type,size,color\right\}$ , we compute the cosine similarity between the estimated HD attribute representation $\hat{\boldsymbol{v}}^{attr}_{j}$ and all possible vectors of $C_{Num}^{Front}$ in the frontend codebooks. $$ W_{j}^{attr}=sim\left(\hat{\boldsymbol{v}}_{j}^{attr},C_{Num}^{Front}\right) \tag{14} $$ where $W_{j}^{attr}(r)$ ( $r∈\{0,1,...,9,null\}$ ) represents the attention weights corresponding to the value $r$ of attribute $attr$ at the $j$ th position, based on the query similarity. Similarly, the attention weights for the attribute $attr∈\left\{exist\right\}$ can be obtained by querying the logic codebook $C_{Lgc}^{Front}$ as follows: $$ W_{j}^{exist}=softmax\left(\beta\cdot sim\left(\hat{\boldsymbol{v}}_{j}^{exist% },C_{Lgc}^{Front}\right)\right) \tag{15} $$ where $W_{j}^{exist}(r)$ ( $r∈\{0,1\}$ ) corresponds to the presence and absence of the object at $j$ th position, respectively. Here, we use the $softmax$ function to normalize the weights, and $\beta$ denotes the inverse softmax temperature. The attention stage: The HD attribute representations required by the reasoning backend involve the entity-level attributes Type, Size, and Color, as well as layout-level attributes Number, Position (Table III). For the numerical attribute $attr∈\left\{type,size,color\right\}$ , the corresponding updated HD representation $\boldsymbol{v}_{j}^{attr}$ can be obtained through the weighted summation on the numerical backend codebook $C_{Num}^{Back}$ as follows: $$ \boldsymbol{v}_{j}^{attr}=\sum_{r\in\left\{0,...,9,null\right\}}{W_{j}^{attr}% \left(r\right)\cdot\boldsymbol{v}\left(r\right)},\boldsymbol{v}\left(r\right)% \in C_{Num}^{Back} \tag{16} $$ For the logical existence attribute $attr∈\left\{exist\right\}$ , its updated HD representation $\boldsymbol{v}_{j}^{exist}$ are obtained through the weighted summation on the backend codebook $C_{Lgc,RV}^{Back}$ and $C_{Lgc,BV}^{Back}$ , respectively, as follows: $$ \boldsymbol{v}_{j}^{exist,VT}=\sum_{r\in\left\{0,1\right\}}{W_{j}^{exist}\left% (r\right)\cdot\boldsymbol{e}\left(r\right)},\boldsymbol{e}\left(r\right)\in C_% {Lgc,VT}^{Back} \tag{17} $$ where $VT∈\{RV,BV\}$ represents the type of atomic HD vectors. Additionally, we introduce the overall HD attribute representations for the Type, Size, and Color attributes within the nxn Grid. These representations $\boldsymbol{v}_{\mathrm{nxn}}^{attr}$ ( $attr∈\left\{type,size,color\right\}$ ) can be obtained by bundling corresponding HD attribute representations of all objects in the Grid with their attention weights of existence as follows: $$ \boldsymbol{v}_{\mathrm{nxn}}^{attr}=\sum_{j=1}^{n^{2}}{W_{j}^{exist}\left(1% \right)\cdot\boldsymbol{v}_{j}^{attr}} \tag{1} $$ For the layout-level attribute Number, its HD attribute representation $\boldsymbol{v}^{number}$ is obtained by projecting the sum of the attention weights of presence to FPE-VSA as follows: $$ \boldsymbol{v}^{number}=\boldsymbol{v}^{\left(\otimes\sum_{j=1}^{n^{2}}{W_{j}^% {exist}\left(1\right)}\right)} \tag{1} $$ where $\boldsymbol{v}$ is the base vector of the numerical backend codebook $C_{Num}^{Back}$ . The layout-level attribute Position within the nxn Grid involves both numerical and logical rules (Table III). Therefore, its HD attribute representations correspond to two distinct rules: the logical representation of each individual object $\boldsymbol{v}_{j}^{position}$ and the overall HD position representation $\mathcal{C}^{position}$ of the entire nxn Grid. The former is an HD existence representation with logical computational properties, that is, $\boldsymbol{v}_{j}^{position}=\boldsymbol{v}_{j}^{exist,BV}$ . Inspired from SHDR for the nxn Grid in Equation 12, the overall HD position representation $\mathcal{C}^{position}$ can be obtained as follows: $$ \mathcal{C}^{position}_{\mathrm{nxn}}=\sum_{j=1}^{n\times n}{W_{j}^{exist}% \left(1\right)\cdot\boldsymbol{p}_{j}\circ\boldsymbol{v}_{j}^{exist,RV}}\,,% \boldsymbol{p}_{j}\in C_{Pos,\mathrm{nxn}}^{Back} \tag{1} $$ IV-E 2 Reasoning backend The Rel-SAR model efficiently implements systematic abductive reasoning by leveraging HD attribute representations and VSA-based relation functions. HD attribute representations from the frontend are transformed into the HD vector space, enabling VSA operations on both numerical and logical relation functions (Equation 6 - 9). Consequently, the reasoning backend can perform systematic rule abduction and execution based on these relational functions, without requiring extensive use of explicit rule templates. TABLE V: Attribute sets for n-ary relation. | $\mathbb{V}^{2}$ | $\left(\boldsymbol{v}_{\left(1,1\right)},\boldsymbol{v}_{\left(1,2\right)}\right)$ , $\left(\boldsymbol{v}_{\left(1,2\right)},\boldsymbol{v}_{\left(1,3\right)}\right)$ , $\left(\boldsymbol{v}_{\left(2,1\right)},\boldsymbol{v}_{\left(2,2\right)}\right)$ , | | --- | --- | | $\left(\boldsymbol{v}_{\left(2,2\right)},\boldsymbol{v}_{\left(2,3\right)}\right)$ , $\left(\boldsymbol{v}_{\left(3,1\right)},\boldsymbol{v}_{\left(3,2\right)}\right)$ ; $\left(\boldsymbol{v}_{\left(3,2\right)},\boldsymbol{v}_{\left(y\right)}\right)$ | | | $\mathbb{V}^{3}$ | $\left(\boldsymbol{v}_{\left(1,1\right)},\boldsymbol{v}_{\left(1,2\right)},% \boldsymbol{v}_{\left(1,3\right)}\right)$ , $\left(\boldsymbol{v}_{\left(2,1\right)},\boldsymbol{v}_{\left(2,2\right)},% \boldsymbol{v}_{\left(2,3\right)}\right)$ ; | | $\left(\boldsymbol{v}_{\left(3,1\right)},\boldsymbol{v}_{\left(3,2\right)},% \boldsymbol{v}_{\left(y\right)}\right)$ | | Rule Abduction. Attributes in the RAVEN dataset follow row-major binary or ternary relations [21]. All possible binary $\mathbb{V}^{2}$ and ternary $\mathbb{V}^{3}$ relation pairs in the RPM test are presented in Table V, where $\boldsymbol{v}_{(i,j)}$ denotes the HD attribute representation for a given attribute in the context panel at row $i$ and column $j$ , and $\boldsymbol{v}_{\left(y\right)}$ represents the corresponding HD attribute representation of the target answer panel. Consequently, rule abduction can be formulated as an optimization problem: For both numerical and logical rules, the rule abduction module must identify a set of operator powers $OP_{1:M}$ such that all $N$ -ary ( $N=2,3$ ) relation pairs yield the same output $\boldsymbol{r}_{Num/Lgc}$ when processed through their respective relation functions $R_{Num/Lgc}$ . Formally: $$ \underset{OP_{1:M}}{\max}s^{N}=\prod_{\mathbb{V}_{i}^{N},\mathbb{V}_{j}^{N}\in% \mathbb{V}^{N}}^{i\neq j}{sim\left(R\left(\mathbb{V}_{i}^{N},OP_{1:M}\right),R% \left(\mathbb{V}_{j}^{N},OP_{1:M}\right)\right)} \tag{21} $$ where $sim$ denotes cosine similarity, and $R$ represents either the numerical relation function $R_{Num}$ or the logical relation function $R_{Lgc}$ . $\mathbb{V}_{i}^{N},\mathbb{V}_{j}^{N}∈\mathbb{V}^{N}$ ( $i≠ j$ ) refers to any two relation pairs selected from $\mathbb{V}^{N}$ ( $N=2,3$ ). $s^{N}_{Num/Lgc}$ represents the overall similarity between the outputs $r$ of all corresponding relation pairs and can be interpreted as an unnormalized probability of the corresponding rule. Based on the above idea, in the rule abduction module (Figure 4 c Left), the HD attribute representations of $8$ context panels for a given numerical attribute are input into a trainable neural network $f_{\phi}^{Num}$ as the rule learner to predict the operator powers $\widehat{OP}_{1:M}$ , which are expected to achieve the optimization of the objective defined in Equation 21. Subsequently, all attribute sets for binary and ternary relations (Table V) are input into the corresponding numerical relation functions using the predicted $\widehat{OP}_{1:M}$ , and their outputs $\hat{\boldsymbol{r}}_{Num}$ are obtained. Based on the outputs $\hat{\boldsymbol{r}}_{Num}$ from the relation functions, the unnormalized probability $s^{2}_{Num}$ and $s^{3}_{Num}$ for binary and ternary relations, respectively, can be computed (Equation 21). The operator powers $\widehat{OP}_{1:M}$ and the averaged output $\bar{\boldsymbol{r}}_{Num}$ with larger $s^{N}_{Num}$ are then defined as the underlying numerical rule. Similarly, another trainable rule learner, $f_{\varphi}^{Lgc}$ , is used to predict the operator powers $\widehat{OP}_{1:M}$ for the logical rules associated with the attribute Position (Figure 4 c Right). Since the logical rules in the Raven dataset only involve ternary relations, all logical representations are organized according to the attribute sets of ternary relations outlined in Table V, and are then input into the logical ternary relation functions with the predicted parameters $\widehat{OP}_{1:M}$ . Based on the outputs $\hat{\boldsymbol{r}}_{Lgc}$ from the relation functions, the unnormalized probability $s^{3}_{Lgc}$ for logical relations can be computed (Equation 21). Subsequently, $s^{3}_{Lgc}$ is compared with $s^{N}_{Num}$ , which corresponds to numerical relations for the attribute Position. If $s^{3}_{Lgc}$ is larger, the operator powers $\widehat{OP}_{1:M}$ and the averaged output $\bar{\boldsymbol{r}}_{Lgc}$ for logical relations can be interpreted as the underlying logical rule. Rule Execution. After obtaining the predicted operator powers ${\widehat{OP}_{1:M}}$ and the outputs $\bar{\boldsymbol{r}}_{Num/Lgc}$ that represent the rules, we apply these rules to infer the HD attribute representations of the missing panel. For a given attribute, the corresponding attribute representations from the first two panels in the third row of the RPM test are input into the inverse numerical and logical relation functions using the predicted ${\widehat{OP}_{1:M}}$ and $\bar{\boldsymbol{r}}_{Num/Lgc}$ , resulting in the retrieval of the missing HD attribute representation $\hat{\boldsymbol{v}}_{\left(3,3\right)}$ : $$ \hat{\boldsymbol{v}}_{\left(3,3\right)}=R_{Num/Lgc}^{-1}\left(\mathbb{V},% \widehat{OP}_{1:M},\bar{\boldsymbol{r}}_{Num/Lgc}\right) \tag{22} $$ where $\mathbb{V}=(\boldsymbol{v}_{\left(3,2\right)})$ for binary relation rules and $(\boldsymbol{v}_{\left(3,1\right)},\boldsymbol{v}_{\left(3,2\right)})$ for ternary relation rules. The final answer selection. Finally, we calculate the similarity between all HD attribute representations of the missing panel ( $\boldsymbol{v}^{position}_{\left(3,3\right)},\boldsymbol{v}^{number}_{\left(3,% 3\right)},\boldsymbol{v}^{type}_{\left(3,3\right)},\boldsymbol{v}^{size}_{% \left(3,3\right)},\boldsymbol{v}^{color}_{\left(3,3\right)}$ ) and the corresponding attribute representations of each candidate panel $\left(y\right)$ . The predicted answer panel $\hat{y}$ is the one with the highest total similarity score. IV-F Model training IV-F 1 End-to-end training During end-to-end training, the Rel-SAR model utilizes the rule labels provided by the RAVEN dataset and the answer labels to optimize the objectives of visual perception (Equation 13) and rule abduction (Equation 21). Based on the rule labels, the corresponding ground-truth $OP_{1:M}^{gt}$ and $\boldsymbol{r}_{Num/Lgc}^{gt}$ , which represent the rules, can be obtained from Table IV. To facilitate the learning of $OP_{1:M}$ , we design the loss function $\mathcal{L}_{op}$ , which constrains the rule learners $f_{\phi}^{Num}$ and $f_{\varphi}^{Lgc}$ to optimize the objective described in Equation 21: $$ \mathcal{L}_{op}=MSE\left(\widehat{OP}_{1:M},OP_{1:M}^{gt}\right) \tag{23} $$ which is a mean square error ( $MSE$ ) loss between the predicted operator powers $\widehat{OP}_{1:M}$ and the corresponding ground-truth. Additionally, we introduce the loss function $\mathcal{L}_{\boldsymbol{r}}$ to ensure consistent outputs when the inputs to the relation function follow a given rule. This is formulated as follows: $$ \mathcal{L}_{\boldsymbol{r}}=\sum_{i}{\left(1-sim\left(\hat{\boldsymbol{r}}_{i% },\boldsymbol{r}^{gt}\right)\right)} \tag{24} $$ where $\widehat{\boldsymbol{r}}_{i}$ denotes the output of the relation function for the $i$ -th relation pair in $\mathbb{V}^{2}$ and $\mathbb{V}^{3}$ . The overall loss function $\mathcal{L}$ for end-to-end training is constructed as: $$ \mathcal{L}=\sum{\mathcal{L}_{op}}+\sum{\mathcal{L}_{\boldsymbol{r}}} \tag{25} $$ where $\sum$ represents the sum of the loss functions across all attributes, for both binary and ternary relation functions, and for both numerical and logical rule types. By minimizing $\mathcal{L}_{op}$ and $\mathcal{L}_{\boldsymbol{r}}$ simultaneously during training, the optimization objective of the perception network $f_{\theta}$ , as described in Equation 13, can be achieved. This is because, under the constraints of numerous ground-truth rules ( $OP_{1:M}^{gt}$ and $\boldsymbol{r}_{Num/Lgc}^{gt}$ ), the expected theoretical SHDR constructed from HD attribute representations in the codebooks will be a competitive representation, guiding the reasoning process toward optimality ( $\mathcal{L}_{op}→ 0$ and $\mathcal{L}_{\boldsymbol{r}}→ 0$ ). IV-F 2 End-to-End Training with auxiliary attribute labels Following previous work, we assess the performance of the Rel-SAR model in end-to-end training using both auxiliary attribute labels and answer labels. Here, a cosine similarity loss is employed as the perception loss function to enhance the similarity between the estimated SHDR $\hat{\mathcal{S}}^{ind}$ of the perception frontend and the theoretical SHDR $\mathcal{S}$ (Equation 10) derived from attribute labels, thereby optimizing the trainable perception network $f_{\theta}$ (Equation 13). The perception loss function $\mathcal{L}_{p}$ is defined as follows: $$ \mathcal{L}_{p}=1-sim\left(\hat{\mathcal{S}}^{ind},\mathcal{S}\right) \tag{26} $$ Meanwhile, to achieve rule learning through the optimization objective described in Equation 21, we introduce the loss function $\mathcal{L}_{rs}$ to increase the overall similarity between the outputs $\boldsymbol{r}$ of all corresponding relation pairs in $\mathbb{V}^{2}$ and $\mathbb{V}^{3}$ (Table V), respectively. This is formulated as follows: $$ \mathcal{L}_{rs}=1-\prod_{i,j}^{i\neq j}{sim\left(\hat{\boldsymbol{r}}_{i},% \hat{\boldsymbol{r}}_{j}\right)} \tag{27} $$ where $\hat{\boldsymbol{r}}_{i}$ and $\hat{\boldsymbol{r}}_{j}$ denotes the outputs of relation functions for any two distinct sets of relation pairs in $\mathbb{V}^{2}$ and $\mathbb{V}^{3}$ , including attribute sets from the answer panel $y$ . Therefore, the overall loss function $\mathcal{L}$ for end-to-end training with auxiliary attribute labels is constructed as: $$ \mathcal{L}=\sum{\mathcal{L}_{p}}+\sum{\mathcal{L}_{rs}} \tag{28} $$ where the former $\sum$ represents the sum of the loss functions across all context panels and the answer panel, while the latter $\sum$ represents the sum of the loss functions across all attributes. Additionally, similar to other neuro-symbolic approaches, the perception frontend and reasoning backend in Rel-SAR can also be trained independently using the loss functions $\mathcal{L}_{p}$ and $\mathcal{L}_{rs}$ , respectively. V Experiments V-A Datasets & Experiment setup We evaluate our model on the RAVEN [21] and I-RAVEN [8] datasets. Each dataset contains 70,000 RPM samples, which are divided into training, validation, and test sets with a 6:2:2 ratio. We use a ResNet-50 as the encoder ( $f_{\theta}$ ) to map the image panels $\mathcal{X}$ to their SHDR $\mathcal{S}∈\mathbb{R}^{d}$ , where the dimension $d$ of all vectors is set to $3000$ . Two 5-layer fully connected networks as rule learners ( $f_{\phi}^{Num}$ and $f_{\varphi}^{Lgc}$ ) are employed to extract $OP_{1:M}$ from updated Position representation and other attribute (Number, Type, Size and Color) representation, respectively. The existence vectors of the backend codebook are set to HRR because the normalization of FHRR in the superposition operation will invalidate the logical reasoning. All other vectors are FHRR. We utilize the AdamW optimizer with a learning rate of $1× 10^{-4}$ and a weight decay of $1× 10^{-4}$ . TABLE VI: Test accuracy (%) on RAVEN and I-RAVEN Dataset | Method | RAVEN | I-RAVEN | | | | | | | | | | | | | | | | | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | | Avg | Center | 2X2 | 3X3 | L-R | U-D | O-IC | O-IG | Avg | Center | 2X2 | 3X3 | L-R | U-D | O-IC | O-IG | | | | Deep Learning Model | WReN [4] | 14.7 | 13.1 | 28.6 | 28.3 | 7.5 | 6.3 | 8.4 | 10.6 | 23.8 | 29.4 | 26.8 | 23.5 | 21.9 | 21.4 | 22.5 | 21.5 | | LEN [7] | 72.9 | 80.2 | 57.5 | 62.1 | 73.5 | 81.2 | 84.4 | 71.5 | 41.4 | 56.4 | 31.7 | 29.7 | 44.2 | 44.2 | 52.1 | 31.7 | | | CoPINet [6] | 91.4 | 95.1 | 77.5 | 78.9 | 99.1 | 99.7 | 98.5 | 91.4 | 46.1 | 54.4 | 36.8 | 31.9 | 51.9 | 52.5 | 52.2 | 42.8 | | | MXGNet [39] | 84.0 | 94.3 | 60.5 | 64.9 | 96.6 | 96.4 | 94.1 | 81.3 | 33.1 | 40.7 | 27.9 | 24.7 | 35.8 | 34.5 | 36.4 | 31.6 | | | SCL [28] | 91.6 | 98.1 | 91.0 | 82.5 | 96.8 | 96.5 | 96.0 | 80.1 | 95.0 | 99.0 | 96.2 | 89.5 | 97.9 | 97.1 | 97.6 | 87.7 | | | SRAN [8] | 54.3 | — | — | — | — | — | — | — | 60.8 | 78.2 | 50.1 | 42.4 | 70.1 | 70.3 | 68.2 | 46.3 | | | Rel-Base [30] | 91.7 | 97.6 | 85.9 | 86.9 | 93.5 | 96.5 | 97.6 | 83.8 | 91.1 | — | — | — | — | — | — | — | | | MRNet [9] | 74.4 | 96.2 | 49.1 | 45.9 | 93.7 | 94.2 | 92.5 | 51.3 | 75.0 | 96.8 | 45.6 | 39.6 | 95.7 | 95.9 | 95.6 | 55.5 | | | DCNet [25] | 93.6 | 97.8 | 81.7 | 86.7 | 99.8 | 99.8 | 99.0 | 91.5 | 49.4 | 57.8 | 34.1 | 35.5 | 58.5 | 60.0 | 57.0 | 42.9 | | | STSN [29] | 89.7 | — | — | — | — | — | — | — | 95.7 | 98.6 | 96.2 | 88.8 | 98.0 | 98.8 | 97.8 | 92.0 | | | PredRNet [27] | 95.8 | — | — | — | — | — | — | — | 96.5 | — | — | — | — | — | — | — | | | DRNet [26] | 96.9 | — | — | — | — | — | — | — | 97.6 | — | — | — | — | — | — | — | | | Neuro-Symbolic Model | PrAE [10] | 65.0 | 76.5 | 78.6 | 28.6 | 90.1 | 90.9 | 48.1 | 42.6 | 77.0 | 90.5 | 85.4 | 45.6 | 96.3 | 97.4 | 63.5 | 60.7 | | ALANS [13] | 74.4 | 69.1 | 80.2 | 75.0 | 72.2 | 73.3 | 76.3 | 74.9 | 78.5 | 72.3 | 79.5 | 72.9 | 79.2 | 79.6 | 85.9 | 79.9 | | | NVSA [12] | 87.7 | 99.7 | 93.5 | 57.1 | 99.8 | 99.7 | 98.6 | 65.4 | 88.1 | 99.8 | 96.2 | 54.3 | 100 | 99.9 | 99.6 | 67.1 | | | Rel-SAR (Ours) | 96.5 | 99.1 | 95.7 | 96.2 | 99.6 | 99.6 | 99.1 | 86.2 | 98.0 | 99.8 | 97.1 | 98.1 | 100 | 100 | 99.9 | 90.9 | | | Neuro-Symbolic Model (attribute labels) | NVSA [12] | 98.5 | 100 | 99.4 | 96.3 | 100 | 100 | 100 | 93.9 | 99.0 | 100 | 99.5 | 97.1 | 100 | 100 | 100 | 96.4 | | Xu et al. [40] | 92.9 | 98.8 | 91.9 | 93.1 | 99.2 | 99.1 | 98.2 | 70.1 | 93.2 | 99.5 | 89.6 | 89.7 | 99.7 | 99.5 | 99.6 | 74.7 | | | ALANS-V [13] | 94.4 | 98.4 | 91.5 | 87.0 | 97.3 | 96.4 | 97.3 | 93.2 | 93.5 | 98.9 | 85.0 | 83.2 | 90.9 | 98.1 | 99.1 | 89.5 | | | Rel-SAR (Ours) | 96.6 | 97.9 | 94.3 | 96.6 | 99.0 | 98.7 | 97.9 | 92.2 | 98.1 | 98.5 | 96.7 | 97.8 | 99.4 | 99.8 | 99.2 | 95.6 | | V-B End-to-end learning results Here we first evaluate the Rel-SAR model on both RAVEN [21] and I-RAVEN [8] datasets using end-to-end learning. The results, summarized in Table VI, compare our model with both deep neural network methods and neuro-symbolic AI methods. Rel-SAR achieves an average accuracy of $96.5\%$ on RAVEN and $98.0\%$ on I-RAVEN, comparable to the previous best-performing deep network, DRNet ( $96.9\%$ on RAVEN and $97.6\%$ on I-RAVEN). Compared to the previous neuro-symbolic method NVSA, Rel-SAR demonstrates a significant accuracy improvement on configurations involving rules based on the position attribute, including 2x2Grid ( $+2.2\%$ ), 3x3Grid ( $+39.1\%$ ), and Out-InGrid ( $+20.8\%$ ). This results in an average accuracy improvement of $8.8\%$ on RAVEN and $9.9\%$ on I-RAVEN. These improvements are attributed to the effectiveness of our circular and logical HD representations for the position attribute, as well as the numerical and logical relation functions used in rule abduction and execution. Additionally, Rel-SAR utilizes only the answer panel during training, without leveraging information from other candidate panels. This ensures that Rel-SAR does not exploit potential defects in the RAVEN dataset for shortcut learning [8], resulting in similar accuracy on both RAVEN and I-RAVEN datasets ( $96.5\%$ vs. $98.0\%$ ). In addition, since neuro-symbolic models are hybrid architectures that disentangle perception and reasoning, they are often trained with auxiliary attribute labels. We therefore evaluate the Rel-SAR model on both datasets when trained with additional attribute labels. The results, presented in the bottom section of Table VI, show that our model achieves an average accuracy of $96.6\%$ on RAVEN and $98.1\%$ on I-RAVEN. Among the neuro-symbolic models compared, ALANS Learner is the most similar to Rel-SAR, as both models incorporate learnable parameters in their reasoning backends. Compared to ALANS, our model demonstrates an average accuracy improvement of $2.2\%$ on RAVEN and $4.6\%$ on I-RAVEN. In contrast, the reasoning backends in NVSA and Xu’s model rely on pre-designed rule templates or algebraic invariance modules without learnable parameters. NVSA, with its comprehensive rule template library, achieves near-perfect accuracy when trained with auxiliary attribute labels. Although our model slightly underperforms NVSA when trained with auxiliary attribute labels, it exhibits an important advantage: when trained without attribute labels, both NVSA and ALANS Learner suffer significant performance drops, whereas our model maintains nearly identical performance. This consistency highlights the robust synergy between the perception and reasoning modules in our model. V-C Perception results learned with the attribute labels TABLE VII: The visual perception accuracy on RAVEN dataset. | PrAE | 85.27 | 88.65 | 93.56 | 73.95 | 100 | 100 | 94.23 | 46.25 | | --- | --- | --- | --- | --- | --- | --- | --- | --- | | Xu et al. | 96.10 | 100 | 100 | 99.99 | 100 | 100 | 99.96 | 72.78 | | NVSA | 99.76 | 100 | 99.83 | 98.61 | 99.97 | 99.96 | 99.97 | 99.95 | | Rel-SAR | 99.99 | 100 | 100 | 99.99 | 100 | 100 | 100 | 99.94 | Similar to other neuro-symbolic methods, the perception frontend in Rel-SAR can also be independently trained with attribute labels using the perception loss function (Equation 26). By querying the estimated HD attribute representations for the highest similarity with the attribute vectors in the frontend codebooks, we can retrieve the predicted entity-level attributes. A panel’s features are considered correctly extracted only if the attributes of all objects in the panel match those provided by the dataset. As shown in Table VII, the evaluation results demonstrate that the perception frontend in our model achieves an average panel accuracy of $99.99\%$ on the RAVEN dataset, maintaining nearly perfect performance, consistent with NVSA. Notably, NVSA reports a resolution issue in the RAVEN dataset, where some objects in the inner regions of the O-InGrid configuration have a different size attribute but the same image representation [12]. To address this, we adopt a similar solution to NVSA, merging classes with different sizes but identical panel representation. To evaluate the generalization capability of the perception frontend of our model, we adopt the experimental settings outlined in [12] to test the model on unseen combinations of attribute values. Specifically, we focus on the single-object case in the 2x2Grid configuration. As shown in Table VIII, we select two attributes (e.g. Position and Type) along with partial value sets for each (e.g. Position $∈\{0,3\}$ and Type $∈\{0,2\}$ ) as the target attributes. Panels containing these target attribute values are included in the training set, while panels lacking them are designated for the test set. This approach ensures that the attribute sets of the training and test datasets are entirely disjoint, enabling a rigorous evaluation of the module’s ability to generalize beyond the observed data. In this experiment, we observe that our perception frontend struggles to identify unseen combinations of attribute values when the value vectors $\boldsymbol{v}$ in Equation 10 are Random Vectors (RVs). We attribute this issue to the orthogonality of RVs, which hinders the model’s ability to recognize adjacency concepts. To address this, we replace the RVs with fractional power encoding (FPE) vectors based on a Gaussian kernel, which introduces progressive similarity between vectors [18, 19]. As shown in Table VIII, our model demonstrates better generalization on unseen attribute combinations compared to NVSA [12] when using FPE vectors with progressive similarity (NVs). However, the type-size combination still shows low accuracy, likely because there is no clear continuous progression between adjacent attribute values of type (e.g., triangle and square). We also conduct an experiment to evaluate the generalization ability of our perception frontend when applied to an unseen number of objects [12]. The training set consists of panels with a finite number of objects (e.g., $k_{train}=1$ ), while the test set consists of panels with a larger number of objects (e.g., $k_{test}=2$ ). As shown in Table IX, when $k_{train}=1$ , our model fails to accurately predict the attributes of panels containing more objects. By analyzing the predicted attribute values, we observe that our model exhibits a consistent distribution of predicted attribute attention weights across different test sets, attributable to the limited diversity (Figure 5 a). However, as the number of objects in the training panels increases, the model learns more complex patterns and can correctly allocate attention to different positions (Figure 5 b). Our model achieves perfect generalization performance when $k_{train}=2$ for 2x2Grid and when $k_{train}=3$ for 3x3Grid. TABLE VIII: Accuracy of attribute-value generalization on 2x2Grid containing $k=1$ object | Position $∈\{0,3\}$ OR | Position $∉\{0,3\}$ AND | 6720 | 2880 | 0.0 | 0.0 | 26.8 | | --- | --- | --- | --- | --- | --- | --- | | Type $∈\{0,2\}$ | Type $∉\{0,2\}$ | | | | | | | Position $∈\{0,3\}$ OR | Position $∉\{0,3\}$ AND | 6400 | 3200 | 15.1 | 0.0 | 31.4 | | Size $∈\{1,5\}$ | Size $∉\{1,5\}$ | | | | | | | Position $∈\{0,3\}$ OR | Position $∉\{0,3\}$ AND | 6720 | 2880 | 34.8 | 0.0 | 73.5 | | Color $∈\{0,3,6,8\}$ | Color $∉\{0,3,6,8\}$ | | | | | | | Type $∈\{0,2\}$ OR | Type $∉\{0,2\}$ AND | 5760 | 3840 | 0.0 | 0.0 | 13.9 | | Size $∈\{1,5\}$ | Size $∉\{1,5\}$ | | | | | | | Type $∈\{0,2\}$ OR | Type $∉\{0,2\}$ AND | 6144 | 3456 | 72.0 | 0.0 | 89.5 | | Color $∈\{0,3,6,8\}$ | Color $∉\{0,3,6,8\}$ | | | | | | | Size $∈\{1,5\}$ OR | Size $∉\{1,5\}$ AND | 5760 | 3840 | 29.3 | 0.0 | 82.6 | | Color $∈\{0,3,6,8\}$ | Color $∉\{0,3,6,8\}$ | | | | | | TABLE IX: Generalization to a growing number of unseen objects in the RAVEN panel | 2x2Grid | 9600 | 1 | 0.3 | 0.0 | 0.0 | — | — | — | — | — | 0.1 | 0.1 | | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | | 2x2Grid | 19200 | 2 | 100.0 | 100.0 | 100.0 | — | — | — | — | — | 100.0 | 100.0 | | 3x3Grid | 21600 | 1 | 0.2 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | | 3x3Grid | 43200 | 2 | 100.0 | 100.0 | 99.4 | 70.3 | 21.6 | 0.4 | 0.0 | 0.0 | 49.0 | 41.7 | | 3x3Grid | 64800 | 3 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | | 3x3Grid | 86400 | 4 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | 100.0 | <details> <summary>x5.png Details</summary>  ### Visual Description ## Density Plots: Sum of Existence Weights vs. k_train and k_test ### Overview The image contains two density plots comparing the distribution of "Sum of Existence Weights" for different combinations of `k_train` (training parameter) and `k_test` (testing parameter). Each plot uses colored curves to represent distinct `k_test` values, with density on the y-axis and sum of weights on the x-axis. --- ### Components/Axes - **Chart (a)**: `k_train=1` - **X-axis**: "Sum of Existence Weights" (range: 0.7–1.2) - **Y-axis**: "Density" (range: 0–8) - **Legend**: - Red: `k_test=2` - Green: `k_test=3` - Purple: `k_test=4` - **Title**: "(a) k_train=1" - **Chart (b)**: `k_train=2` - **X-axis**: "Sum of Existence Weights" (range: 2.0–4.0) - **Y-axis**: "Density" (range: 0–6) - **Legend**: - Red: `k_test=2` - Green: `k_test=3` - Purple: `k_test=4` - **Title**: "(b) k_train=2" --- ### Detailed Analysis #### Chart (a): `k_train=1` - **Trends**: - All curves peak near **0.9–1.0**, with overlapping distributions. - `k_test=2` (red) has the highest peak density (~6) and the widest spread. - `k_test=3` (green) peaks slightly lower (~5) but overlaps significantly with `k_test=2`. - `k_test=4` (purple) has the lowest peak (~4) and a narrower distribution. - Tails extend to the right, with `k_test=2` having the longest tail (up to 1.2). - **Key Values**: - `k_test=2`: Peak density ~6 at ~0.95. - `k_test=3`: Peak density ~5 at ~0.98. - `k_test=4`: Peak density ~4 at ~1.0. #### Chart (b): `k_train=2` - **Trends**: - Curves are more distinct and separated compared to Chart (a). - `k_test=2` (red) peaks sharply at **2.0** with a density of ~5. - `k_test=3` (green) peaks at **3.0** with a density of ~4. - `k_test=4` (purple) peaks at **3.5** with a density of ~3. - No significant overlap between curves. - **Key Values**: - `k_test=2`: Peak density ~5 at 2.0. - `k_test=3`: Peak density ~4 at 3.0. - `k_test=4`: Peak density ~3 at 3.5. --- ### Key Observations 1. **Overlap vs. Separation**: - For `k_train=1`, distributions overlap heavily, suggesting similar behavior across `k_test` values. - For `k_train=2`, distributions are distinct, indicating divergence in behavior as `k_train` increases. 2. **Peak Shifts**: - As `k_train` increases from 1 to 2, the peak of the sum of existence weights shifts rightward (from ~0.9–1.0 to 2.0–3.5). 3. **Density Magnitude**: - Maximum densities decrease from Chart (a) to (b), suggesting reduced concentration of weights as `k_train` increases. --- ### Interpretation - **Model Behavior**: The plots likely represent a machine learning context (e.g., k-nearest neighbors), where `k_train` and `k_test` influence data weighting. The overlap in Chart (a) implies that small changes in `k_test` have minimal impact when `k_train` is low. In Chart (b), increased `k_train` leads to clearer separation, suggesting stronger sensitivity to `k_test` adjustments. - **Practical Implications**: Larger `k_train` values may improve model robustness by reducing ambiguity in weight distributions, while smaller `k_train` values risk overfitting due to overlapping weight distributions. - **Anomalies**: The abrupt separation in Chart (b) could indicate a threshold effect, where `k_train=2` triggers a qualitative shift in how `k_test` values are weighted. --- ### Spatial Grounding & Verification - **Legend Placement**: Both legends are positioned in the upper-right corner of their respective charts, ensuring clear association with curve colors. - **Color Consistency**: Red (`k_test=2`), green (`k_test=3`), and purple (`k_test=4`) are consistently used across both charts, with no mismatches observed. - **Axis Alignment**: X-axis labels and tick marks are centered below each plot, with y-axis labels aligned to the left. --- ### Content Details - **No Text Blocks or Tables**: The image contains only graphical data with no embedded text or tables. - **Uncertainty**: Values are approximate (e.g., "~0.95" for peaks) due to the absence of exact numerical annotations on the curves. </details> Figure 5: Attention weights for existence on 2X2Grid. For the existence attribute, the presence weights ( $W_{j}^{exist}\left(1\right)$ in Equation 15) represent the probability of an object being present at the $j$ th position. The sum of these weights provides an approximate indication of the total number of objects in a panel. (a) When $k_{train}=1$ , the model consistently predicts the presence of only one object in the panel, regardless of the actual test scenario. (b) When $k_{train}=2$ , the model learns to distinguish between panels with varying numbers of objects, effectively adapting its feature extraction process. V-D Reasoning results utilizing the attribute labels We analyze the performance of the reasoning backend of our model on the I-RAVEN dataset. Following the assumption of perfect perception in prior works [14, 31], we also utilize the ground truth attribute labels provided by the I-RAVEN dataset to generate HD attribute representations. The model is trained using the reasoning loss function (Equation 27), and the corresponding evaluation results are shown in Table X. We compare our model with LLM based GPT-3 [41], deep neural network SCL [28], neuro-symbolic based method PrAE [10] and NVSA [12], VSA based method Learn-VRF [31] and ARLC [14]. The performance of GPT-3 is reported in [31]. Our method achieves an accuracy of $99.2\%$ , representing a $6.8\%$ improvement over the previous state-of-the-art model, ARLC. Notably, our approach significantly outperforms prior VSA-based methods ARLC on the 2X2Grid ( $+14.9\%$ ), 3X3Grid ( $+17.7\%$ ), and O-InGrid ( $+9.4\%$ ) configurations, which involve rules on the Position attribute. This improvement demonstrates that the circular and logical HD representations of Position, along with the corresponding relation functions, effectively handle the reasoning of Arithmetic and Progression rules of Position. TABLE X: In-distribution accuracy on the I-RAVEN dataset | GPT-3 [41] | 86.5 | 86.4 | 83.2 | 81.8 | 83.4 | 84.6 | 92.8 | 93.0 | | --- | --- | --- | --- | --- | --- | --- | --- | --- | | SCL [28] | 84.3 | 99.9 | 68.9 | 43.0 | 98.5 | 99.1 | 97.7 | 82.6 | | PrAE [10] | 71.1 | 83.8 | 82.9 | 47.4 | 94.8 | 94.8 | 56.6 | 37.4 | | NVSA [12] | 88.1 | 99.8 | 96.2 | 54.3 | 100 | 99.9 | 99.6 | 67.1 | | Learn-VRF [31] | 79.5 | 97.7 | 56.3 | 49.9 | 94.0 | 95.6 | 98.3 | 64.8 | | ARLC [14] | 92.4 | 98.4 | 83.4 | 80.0 | 98.7 | 98.4 | 98.8 | 89.4 | | Rel-SAR | 99.2 | 99.9 | 98.3 | 97.7 | 100.0 | 99.9 | 100.0 | 98.8 | TABLE XI: Out-of-distribution accuracy on unseen rule-attribute pairs on I-RAVEN | | Type | Size | Color | | | | | | | | | | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | | Const. | Progr. | Dist.3 | Const. | Progr. | Arth. | Dist.3 | Const. | Progr. | Arth. | Dist.3 | | | GPT-3 [41] | 88.5 | 86.0 | 88.6 | 93.6 | 93.2 | 92.6 | 71.6 | 94.2 | 94.7 | 94.3 | 65.8 | | Learn-VRF [31] | 100 | 100 | 99.7 | 100 | 100 | 99.8 | 99.8 | 100 | 98.8 | 100 | 100 | | ARLC [14] | 100 | 98.6 | 99.7 | 100 | 100 | 99.6 | 99.6 | 100 | 100 | 100 | 99.8 | | Rel-SAR | 99.7 | 99.6 | 100 | 99.3 | 99.6 | 97.9 | 100 | 96.8 | 89.5 | 86.4 | 70.7 | We also evaluate the out-of-distribution (OOD) generalization ability of our reasoning backend, following the experimental setup described in [31]. In this evaluation, a specific rule (e.g., Type) for an attribute (e.g., Const.) is designated as the target rule. The model is trained and validated using samples that exclude the target rule and is then tested on samples containing only the target rule. This experiments is conducted on the Center configuration, while training, validation, and test sets are filtered from the I-RAVEN dataset. The corresponding evaluation results are shown in Table XI. For attributes Type and Size, our model, like Learn-VRF [31] and ARLC [14], demonstrates a near-perfect ability to generalize to unseen attribute rules. This capability arises from the use of unified HD vectors to represent values across different attributes, facilitating rule transfer between attributes. However, our model exhibits relatively limited performance on the Color attribute. This could be due to the fact that Color has a broader range of attribute values, making it difficult to transfer rules learned from attributes with narrower value ranges, such as Type (0-4) and Size (0-5), to the Color (0-9) attribute. VI Conclusion and Future Directions In this work, we propose Rel-SAR, a novel model that leverages VSA algebra to facilitate systematic rule abduction and execution. Rel-SAR adopts a neuro-symbolic architecture, where the perception frontend extracts diverse high-dimensional attribute representations with intrinsic algebraic properties, and the reasoning backend systematically derives a variety of rules based on relation functions. Extensive experiments demonstrate that Rel-SAR achieves superior performance while offers better interpretability and transparency on RPM tasks. The perception frontend of Rel-SAR effectively extracts object-level attributes while preserving the structure information of the image through the binding and bundling operations of VSA. By utilizing fractional power encoding vectors with progressive similarity, our model demonstrates its ability to generalize beyond the observed data. However, Rel-SAR exhibits relatively low accuracy on the Out-InGrid configuration, which may be attributed to the small object size and resolution issue. As noted in [27], CNNs with large kernel sizes or more stacked blocks are less effective at extracting features from RPM images with smaller objects. NVSA further shows that reducing the stride from 2 to 1 in the first convolutional block and removing the max-pooling layer in ResNet-18 can improve accuracy [12]. We believe that using CNNs with smaller kernel sizes and strides in the shallow layers will enhance the accuracy on the O-IG configuration. Consistent with other neuro-symbolic methods, we also introduce auxiliary rule labels during training. This is because, without precise rule-driven guidance, the model struggles to learn meaningful structured high-dimensional representations (SHDR) of attributes. We note that slot attention [42] enables unsupervised scene decomposition, while VQ-VAE [43] learns discrete latent representations, disentangling different concepts (i.e., attributes in RPM) from raw images. Therefore, we suggest employing a learnable frontend codebook, combined with slot attention to enable unsupervised extraction of SHDR from raw images in the future work. The reasoning backend of Rel-SAR implements sytematic abductive reasoning based on diverse HD attribute representation and relation function. Our model has a significant improvement on configurations with rules on Position. This demonstrates the effectiveness of circular and logical representations of Position attribute. The reasoning backend of the Rel-SAR exhibits limited out-of-distribution (OOD) generalization ability. This limitation may stem from the rule learner, a multi-layer fully connected neural network, which lacks the capacity to generalize to OOD attribute values. A growing body of recent work emphasizes reasoning based on relations between perceptual inputs rather than the features of individual inputs [32, 33, 34]. This trend is encapsulated by the ”relational bottleneck” principle [35], which aims to mitigate the influence of the absolute size of attribute values on relational reasoning. However, relational bottleneck may struggle with more complex relations. Future studies can focus on combining relational bottlenecks with VSA algebra to handle complex relations while maintaining robust out-of-distribution generalization. References - [1] P. A. Carpenter, M. A. Just, and P. Shell, “What one intelligence test measures: a theoretical account of the processing in the raven progressive matrices test.” Psychological review, vol. 97, no. 3, p. 404, 1990. - [2] W. B. Bilker, J. A. Hansen, C. M. Brensinger, J. Richard, R. E. Gur, and R. C. Gur, “Development of abbreviated nine-item forms of the raven’s standard progressive matrices test,” Assessment, vol. 19, no. 3, pp. 354–369, 2012. - [3] M. Mitchell, “Abstraction and analogy-making in artificial intelligence,” Annals of the New York Academy of Sciences, vol. 1505, no. 1, pp. 79–101, 2021. - [4] D. Barrett, F. Hill, A. Santoro, A. Morcos, and T. Lillicrap, “Measuring abstract reasoning in neural networks,” in International conference on machine learning. PMLR, 2018, pp. 511–520. - [5] F. Hill, A. Santoro, D. G. Barrett, A. S. Morcos, and T. Lillicrap, “Learning to make analogies by contrasting abstract relational structure,” arXiv preprint arXiv:1902.00120, 2019. - [6] C. Zhang, B. Jia, F. Gao, Y. Zhu, H. Lu, and S.-C. Zhu, “Learning perceptual inference by contrasting,” Advances in neural information processing systems, vol. 32, 2019. - [7] K. Zheng, Z.-J. Zha, and W. Wei, “Abstract reasoning with distracting features,” Advances in Neural Information Processing Systems, vol. 32, 2019. - [8] S. Hu, Y. Ma, X. Liu, Y. Wei, and S. Bai, “Stratified rule-aware network for abstract visual reasoning,” in Proceedings of the AAAI Conference on Artificial Intelligence, vol. 35, no. 2, 2021, pp. 1567–1574. - [9] Y. Benny, N. Pekar, and L. Wolf, “Scale-localized abstract reasoning,” in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, 2021, pp. 12 557–12 565. - [10] C. Zhang, B. Jia, S.-C. Zhu, and Y. Zhu, “Abstract spatial-temporal reasoning via probabilistic abduction and execution,” in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, 2021, pp. 9736–9746. - [11] G. F. Marcus, The algebraic mind: Integrating connectionism and cognitive science. MIT press, 2003. - [12] M. Hersche, M. Zeqiri, L. Benini, A. Sebastian, and A. Rahimi, “A neuro-vector-symbolic architecture for solving raven’s progressive matrices,” Nature Machine Intelligence, vol. 5, no. 4, pp. 363–375, 2023. - [13] C. Zhang, S. Xie, B. Jia, Y. N. Wu, S.-C. Zhu, and Y. Zhu, “Learning algebraic representation for systematic generalization in abstract reasoning,” in European Conference on Computer Vision. Springer, 2022, pp. 692–709. - [14] G. Camposampiero, M. Hersche, A. Terzić, R. Wattenhofer, A. Sebastian, and A. Rahimi, “Towards learning abductive reasoning using vsa distributed representations,” in International Conference on Neural-Symbolic Learning and Reasoning. Springer, 2024, pp. 370–385. - [15] F. A. Mansouri, D. J. Freedman, and M. J. Buckley, “Emergence of abstract rules in the primate brain,” Nature Reviews Neuroscience, vol. 21, no. 11, pp. 595–610, 2020. - [16] G. Marcus and E. Davis, “Insights for ai from the human mind,” Communications of the ACM, vol. 64, no. 1, pp. 38–41, 2020. - [17] S. A. Sloman, “The empirical case for two systems of reasoning.” Psychological bulletin, vol. 119, no. 1, p. 3, 1996. - [18] T. A. Plate, Distributed representations and nested compositional structure. Citeseer, 1994. - [19] E. Frady, D. Kleyko, C. Kymn, B. Olshausen, and F. Sommer, “Computing on functions using randomized vector representations.” arXiv: Learning,arXiv: Learning, Sep 2021. - [20] D. Kleyko, D. A. Rachkovskij, E. Osipov, and A. Rahimi, “A survey on hyperdimensional computing aka vector symbolic architectures, part i: Models and data transformations,” ACM Computing Surveys, vol. 55, no. 6, pp. 1–40, 2022. - [21] C. Zhang, F. Gao, B. Jia, Y. Zhu, and S.-C. Zhu, “Raven: A dataset for relational and analogical visual reasoning,” in Proceedings of the IEEE/CVF conference on computer vision and pattern recognition, 2019, pp. 5317–5327. - [22] A. Santoro, D. Raposo, D. G. Barrett, M. Malinowski, R. Pascanu, P. Battaglia, and T. Lillicrap, “A simple neural network module for relational reasoning,” Advances in neural information processing systems, vol. 30, 2017. - [23] M. Jahrens and T. Martinetz, “Solving raven’s progressive matrices with multi-layer relation networks,” in 2020 International Joint Conference on Neural Networks (IJCNN). IEEE, 2020, pp. 1–6. - [24] M. Małkiński and J. Mańdziuk, “Multi-label contrastive learning for abstract visual reasoning,” IEEE Transactions on Neural Networks and Learning Systems, vol. 35, no. 2, pp. 1941–1953, 2022. - [25] T. Zhuo and M. Kankanhalli, “Effective abstract reasoning with dual-contrast network,” arXiv preprint arXiv:2205.13720, 2022. - [26] K. Zhao, C. Xu, and B. Si, “Learning visual abstract reasoning through dual-stream networks,” in Proceedings of the AAAI Conference on Artificial Intelligence, vol. 38, no. 15, 2024, pp. 16 979–16 988. - [27] L. Yang, H. You, Z. Zhen, D. Wang, X. Wan, X. Xie, and R.-Y. Zhang, “Neural prediction errors enable analogical visual reasoning in human standard intelligence tests,” in International Conference on Machine Learning. PMLR, 2023, pp. 39 572–39 583. - [28] Y. Wu, H. Dong, R. Grosse, and J. Ba, “The scattering compositional learner: Discovering objects, attributes, relationships in analogical reasoning,” arXiv preprint arXiv:2007.04212, 2020. - [29] S. S. Mondal, T. Webb, and J. D. Cohen, “Learning to reason over visual objects,” arXiv preprint arXiv:2303.02260, 2023. - [30] S. Spratley, K. Ehinger, and T. Miller, “A closer look at generalisation in raven,” in Computer Vision–ECCV 2020: 16th European Conference, Glasgow, UK, August 23–28, 2020, Proceedings, Part XXVII 16. Springer, 2020, pp. 601–616. - [31] M. Hersche, F. Di Stefano, T. Hofmann, A. Sebastian, and A. Rahimi, “Probabilistic abduction for visual abstract reasoning via learning rules in vector-symbolic architectures,” arXiv preprint arXiv:2401.16024, 2024. - [32] T. W. Webb, I. Sinha, and J. D. Cohen, “Emergent symbols through binding in external memory,” arXiv preprint arXiv:2012.14601, 2020. - [33] A. Altabaa, T. Webb, J. Cohen, and J. Lafferty, “Abstractors and relational cross-attention: An inductive bias for explicit relational reasoning in transformers,” arXiv preprint arXiv:2304.00195, 2023. - [34] G. Kerg, S. Mittal, D. Rolnick, Y. Bengio, B. Richards, and G. Lajoie, “On neural architecture inductive biases for relational tasks. arxiv,” arXiv preprint arXiv:2206.05056, 2022. - [35] T. W. Webb, S. M. Frankland, A. Altabaa, S. Segert, K. Krishnamurthy, D. Campbell, J. Russin, T. Giallanza, R. O’Reilly, J. Lafferty et al., “The relational bottleneck as an inductive bias for efficient abstraction,” Trends in Cognitive Sciences, 2024. - [36] T. A. Plate, “Holographic reduced representations,” IEEE Transactions on Neural networks, vol. 6, no. 3, pp. 623–641, 1995. - [37] P. Kanerva, “Hyperdimensional computing: An introduction to computing in distributed representation with high-dimensional random vectors,” Cognitive computation, vol. 1, pp. 139–159, 2009. - [38] A. Vaswani, “Attention is all you need,” Advances in Neural Information Processing Systems, 2017. - [39] D. Wang, M. Jamnik, and P. Lio, “Abstract diagrammatic reasoning with multiplex graph networks,” arXiv preprint arXiv:2006.11197, 2020. - [40] J. Xu, T. Vaidya, Y. Wu, S. Chandra, Z. Lai, and K. F. E. Chong, “Abstract visual reasoning: An algebraic approach for solving raven’s progressive matrices,” in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, 2023, pp. 6715–6724. - [41] T. B. Brown, “Language models are few-shot learners,” arXiv preprint arXiv:2005.14165, 2020. - [42] F. Locatello, D. Weissenborn, T. Unterthiner, A. Mahendran, G. Heigold, J. Uszkoreit, A. Dosovitskiy, and T. Kipf, “Object-centric learning with slot attention,” Advances in neural information processing systems, vol. 33, pp. 11 525–11 538, 2020. - [43] A. Van Den Oord, O. Vinyals et al., “Neural discrete representation learning,” Advances in neural information processing systems, vol. 30, 2017. Appendix A The full version of logical relation function The RPM-style logic rules are expressed as that the set of attribute values in the third panel of a row or column corresponds to the logic operation applied to the first two panels. The simplified version ( $M=3$ ) of the logical relation function, $R_{Lgc}$ , cannot handle all RPM-style logic rules, such as XOR. To address this limitation, we developed a full version of $R_{Lgc}$ with $OP_{1:M}$ ( $M=5$ ) to describe all RPM-style logical rules that are meaningful, defined as follows: $$ \displaystyle\boldsymbol{r}_{Lgc} \displaystyle=R_{Lgc}\left(\boldsymbol{v}_{1:N},OP_{1:M}\right) \displaystyle=\left(\left(op_{1}\boldsymbol{v}_{1}\land op_{2}\boldsymbol{v}_{% 2}\right)\lor\left(op_{3}\boldsymbol{v}_{1}\land op_{4}\boldsymbol{v}_{2}% \right)\right)\circ op_{5}\boldsymbol{v}_{3} \tag{29} $$ where $op_{i}∈\{0,1\}$ determines whether to negate $\boldsymbol{v}_{i}$ , with negation ( $\lnot$ ) applied when $op_{i}=1$ and no negation applied when $op_{i}=0$ . The implementation of $op_{i}\boldsymbol{v}_{i}$ is as follows: $$ op_{i}\boldsymbol{v}_{i}=\left(\boldsymbol{e}\left(1\right)\right)^{\left(% \circ op_{i}\right)}\circ\boldsymbol{v}_{i} \tag{1} $$ where $\boldsymbol{e}\left(1\right)$ is BVs representing a Boolean value indicating True. When $op_{i}=1$ , $\left(\boldsymbol{e}\left(1\right)\right)^{\left(\circ op_{i}\right)}$ simplifies to $\boldsymbol{e}\left(1\right)$ , which signifies the negation of HD vector $\boldsymbol{v}_{i}$ (see Table II). Similarly, When $op_{i}=0$ , $\left(\boldsymbol{e}\left(1\right)\right)^{\left(\circ op_{i}\right)}$ simplifies to $\boldsymbol{e}\left(0\right)$ , which does not change the $\boldsymbol{v}_{i}$ . The inverse logical relation function is defined as follows: TABLE XII: Logical rules and Corresponding combinations of $OP_{1:M}$ and $r$ in Logical relation functions | AND OR DIFF | 0 1 0 | 0 1 1 | 0 1 0 | 0 1 1 | 0 1 0 | 0 0 0 | | --- | --- | --- | --- | --- | --- | --- | | XOR | 0 | 1 | 1 | 0 | 0 | 0 | $$ \displaystyle\boldsymbol{v}_{N} \displaystyle=R_{Lgc}^{-1}\left(\boldsymbol{v}_{1:N-1},OP_{1:M},\boldsymbol{r}\right) \displaystyle=op_{5}\left(\left(op_{1}\boldsymbol{v}_{1}\land op_{2}% \boldsymbol{v}_{2}\right)\lor\left(op_{3}\boldsymbol{v}_{1}\land op_{4}% \boldsymbol{v}_{2}\right)\right) \tag{31} $$ Several typical logical rules, along with the corresponding $OP_{1:M}$ and $\boldsymbol{r}$ in the full version of logical relation function, are presented in Table XII. The DIFF logic rule, also known as the Arithmetic - rule for the Position attribute of the RAVEN dataset, specifies that the set of attribute values for the third panel is the difference set of the attribute values from the first two panels. OR rule is the Arithmetic + rule for Position attributes in RAVEN. For the OR, AND, and DIFF rules, $op_{3}$ and $op_{4}$ are equal of $op_{1}$ and $op_{2}$ , respectively, allowing them to be represented using the simplified version $R_{Lgc}$ . However, $OP_{1:2}$ and $OP_{3:4}$ of the XOR rule are different, making it impossible to express using the simplified version $R_{Lgc}$ . Because the RAVEN dataset only involve Arithmetic + and Arithmetic - logical rules, this work employs the simplified logical relation function.