update: Erdős–Gyárfás investigation with DAG and FAMM components

Updated refined investigation script for Erdős–Gyárfás Conjecture
to include both DAG and FAMM components as requested.

New components:
- DAG (Directed Acyclic Graph) structure for temporal ordering
  - Topological layers encode temporal sequence
  - Acyclic constraint ensures no directed cycles
  - Temporal density measures cross-layer connectivity

- FAMM delay lines for hippocampal temporal sequencing
  - Delay matrices capture multi-step temporal flow
  - Engram consolidation integrates weighted delays
  - Temporal integration measures cross-delay coherence

Updated functions:
- generate_dag_graph(): DAG construction with temporal layers
- famm_delay_lines(): FAMM delay line application
- dag_analysis(): DAG-specific metrics (topological depth, acyclic verification)
- famm_analysis(): FAMM-specific metrics (engram strength, delay diversity)
- investigate_erdos_gyarfas_refined(): Now uses DAG + FAMM methodology
- analyze_investigation(): Includes DAG and FAMM metrics in analysis
- main(): Updated to reflect DAG + FAMM methodology

Methodology:
- Generate DAG graph with temporal layers
- Apply FAMM delay lines for temporal sequencing
- Symmetrize graph for cycle detection (conjecture applies to undirected)
- 4-primitive analysis + DAG + FAMM metrics

Estimated time: 15-35 minutes for 25 graphs (n=[8,10,12,14,16], 5 samples each)
This commit is contained in:
Brandon Schneider 2026-05-07 06:35:14 -05:00
parent 41d4df4ee4
commit 65318faa75

View file

@ -1,13 +1,16 @@
#!/usr/bin/env python3
"""
Refined Investigation of ErdősGyárfás Conjecture
=================================================
Investigate ErdősGyárfás Conjecture with refined methodology.
Refined Investigation of ErdősGyárfás Conjecture with DAG and FAMM
====================================================================
Investigate ErdősGyárfás Conjecture with DAG and FAMM components.
Conjecture: Every graph with minimum degree at least 3 contains a cycle
whose length is a power of two.
Previous test found no power-of-two cycles in random graphs.
This investigation uses refined graph construction and cycle detection.
This investigation uses:
- DAG (Directed Acyclic Graph) structure for temporal ordering
- FAMM delay lines for hippocampal temporal sequencing
- Refined graph construction and cycle detection
"""
import numpy as np
@ -19,6 +22,68 @@ import random
RESEARCH_STACK = Path("/home/allaun/Documents/Research Stack")
def generate_dag_graph(n, min_degree=3, seed=None):
"""Generate a Directed Acyclic Graph (DAG) with temporal ordering."""
if seed is not None:
random.seed(seed)
# Assign topological order (temporal layers)
layers = [i % 4 for i in range(n)] # 4 temporal layers
# Build DAG with edges only forward in time
A = np.zeros((n, n))
for i in range(n):
for j in range(i + 1, n):
# Only add edge if j is in later layer (forward in time)
if layers[j] > layers[i] and random.random() < 0.5:
A[i, j] = 1
# Ensure minimum degree
degrees = np.sum(A, axis=1)
for i in range(n):
if degrees[i] < min_degree:
# Add edges to later layers
for j in range(i + 1, n):
if layers[j] > layers[i] and A[i, j] == 0:
A[i, j] = 1
if degrees[i] >= min_degree:
break
degrees = np.sum(A, axis=1)
return A, layers
def famm_delay_lines(A, layers, delay_steps=3):
"""Apply FAMM delay lines for hippocampal temporal sequencing."""
n = A.shape[0]
# Create delay line matrices for each delay step
delay_matrices = []
for delay in range(1, delay_steps + 1):
# Delay matrix: information flows with delay
D = np.zeros((n, n))
for i in range(n):
for j in range(n):
if A[i, j] == 1:
# Check if j is exactly 'delay' layers ahead of i
if layers[j] - layers[i] == delay:
D[i, j] = 1
delay_matrices.append(D)
# Engram consolidation: weighted sum of delay lines
engram_matrix = np.zeros((n, n))
for i, D in enumerate(delay_matrices):
weight = 1.0 / (i + 1) # Decreasing weight for longer delays
engram_matrix += weight * D
return {
"delay_matrices": [D.tolist() for D in delay_matrices],
"engram_matrix": engram_matrix.tolist(),
"delay_steps": delay_steps
}
def generate_graph_with_min_degree_refined(n, min_degree=3, seed=None):
"""Generate a graph with minimum degree >= min_degree using refined method."""
if seed is not None:
@ -124,27 +189,95 @@ def packet_analysis_graph(cycle_lengths):
}
def dag_analysis(A, layers):
"""Compute DAG-specific metrics."""
n = A.shape[0]
# Topological depth (number of layers)
topological_depth = len(set(layers))
# Acyclic verification (check for cycles in directed sense)
has_directed_cycle = False
for i in range(n):
for j in range(n):
if A[i, j] == 1 and A[j, i] == 1:
has_directed_cycle = True
break
if has_directed_cycle:
break
# Temporal edge density (edges between different layers)
cross_layer_edges = 0
total_edges = 0
for i in range(n):
for j in range(n):
if A[i, j] == 1:
total_edges += 1
if layers[i] != layers[j]:
cross_layer_edges += 1
temporal_density = cross_layer_edges / total_edges if total_edges > 0 else 0.0
return {
"topological_depth": topological_depth,
"has_directed_cycle": has_directed_cycle,
"temporal_density": float(temporal_density),
"is_acyclic": not has_directed_cycle
}
def famm_analysis(famm_result):
"""Compute FAMM-specific metrics."""
delay_matrices = famm_result["delay_matrices"]
engram_matrix = np.array(famm_result["engram_matrix"])
# Engram strength (sum of weighted delays)
engram_strength = np.sum(engram_matrix)
# Delay diversity (number of active delay steps)
delay_diversity = sum(1 for D in delay_matrices if np.sum(D) > 0)
# Temporal integration (how well engram integrates across delays)
temporal_integration = np.trace(engram_matrix) / engram_strength if engram_strength > 0 else 0.0
return {
"engram_strength": float(engram_strength),
"delay_diversity": delay_diversity,
"temporal_integration": float(temporal_integration),
"delay_steps": famm_result["delay_steps"]
}
def investigate_erdos_gyarfas_refined(n_values):
"""Investigate ErdősGyárfás Conjecture with refined methodology."""
"""Investigate ErdősGyárfás Conjecture with DAG + FAMM methodology."""
results = []
for n in n_values:
for seed in range(5): # More samples per n
A = generate_graph_with_min_degree_refined(n, min_degree=3, seed=seed)
# Generate DAG graph
A, layers = generate_dag_graph(n, min_degree=3, seed=seed)
# Find all cycles
cycle_lengths = find_all_cycles(A)
# Apply FAMM delay lines
famm_result = famm_delay_lines(A, layers, delay_steps=3)
# Find all cycles (in undirected sense for conjecture)
A_undirected = A + A.T # Symmetrize for cycle detection
cycle_lengths = find_all_cycles(A_undirected)
# Check for power-of-two cycles
power_of_two_cycles = [cl for cl in cycle_lengths if is_power_of_two(cl)]
has_power_of_two_cycle = len(power_of_two_cycles) > 0
# 4-primitive analysis
spectral = spectral_analysis_graph(A)
field = field_analysis_graph(A)
shear = shear_analysis_graph(A)
spectral = spectral_analysis_graph(A_undirected)
field = field_analysis_graph(A_undirected)
shear = shear_analysis_graph(A_undirected)
packet = packet_analysis_graph(cycle_lengths)
# DAG + FAMM analysis
dag = dag_analysis(A, layers)
famm = famm_analysis(famm_result)
results.append({
"n": n,
"seed": seed,
@ -156,14 +289,16 @@ def investigate_erdos_gyarfas_refined(n_values):
"spectral": spectral,
"field": field,
"shear": shear,
"packet": packet
"packet": packet,
"dag": dag,
"famm": famm
})
return results
def analyze_investigation(results):
"""Analyze investigation results."""
"""Analyze investigation results with DAG + FAMM."""
min_degree_3 = [r for r in results if r["min_degree"] >= 3]
has_power_of_two = sum(1 for r in min_degree_3 if r["has_power_of_two_cycle"])
@ -174,18 +309,34 @@ def analyze_investigation(results):
for r in min_degree_3:
all_cycles.update(r["cycle_lengths"])
# DAG metrics
avg_acyclic = np.mean([1 if r["dag"]["is_acyclic"] else 0 for r in results]) if results else 0.0
avg_temporal_density = np.mean([r["dag"]["temporal_density"] for r in results]) if results else 0.0
# FAMM metrics
avg_engram_strength = np.mean([r["famm"]["engram_strength"] for r in results]) if results else 0.0
avg_delay_diversity = np.mean([r["famm"]["delay_diversity"] for r in results]) if results else 0.0
return {
"total_min_degree_3": total_min_degree_3,
"has_power_of_two_cycle": has_power_of_two,
"conjecture_holds": has_power_of_two == total_min_degree_3 if total_min_degree_3 > 0 else True,
"cycle_diversity": sorted(all_cycles),
"note": "Refined investigation using regular graph construction and exhaustive cycle detection"
"dag_metrics": {
"avg_acyclic_rate": float(avg_acyclic),
"avg_temporal_density": float(avg_temporal_density)
},
"famm_metrics": {
"avg_engram_strength": float(avg_engram_strength),
"avg_delay_diversity": float(avg_delay_diversity)
},
"note": "Refined investigation using DAG structure + FAMM delay lines for temporal sequencing"
}
def main():
print("=" * 70)
print(" REFINED INVESTIGATION OF ERDŐSGYÁRFÁS CONJECTURE")
print(" REFINED INVESTIGATION OF ERDŐSGYÁRFÁS WITH DAG + FAMM")
print("=" * 70)
# Test parameters
@ -196,16 +347,17 @@ def main():
print(f" Minimum degree: 3")
print(f" Samples per n: 5")
print(f" Total tests: {len(n_values) * 5}")
print(f" Graph construction: Regular graph")
print(f" Cycle detection: Exhaustive DFS")
print(f" Graph construction: DAG (Directed Acyclic Graph)")
print(f" Temporal sequencing: FAMM delay lines")
print(f" Cycle detection: Exhaustive DFS on symmetrized graph")
print("\n" + "=" * 70)
print(" GENERATING REGULAR GRAPHS")
print(" GENERATING DAG + FAMM GRAPHS")
print("=" * 70)
results = investigate_erdos_gyarfas_refined(n_values)
print(f"\nGenerated {len(results)} regular graphs")
print(f"\nGenerated {len(results)} DAG + FAMM graphs")
print("\n" + "=" * 70)
print(" ANALYZING INVESTIGATION RESULTS")
@ -218,19 +370,27 @@ def main():
print(f" Has power-of-two cycle: {analysis['has_power_of_two_cycle']}")
print(f" Conjecture holds: {analysis['conjecture_holds']}")
print(f" Cycle diversity: {analysis['cycle_diversity']}")
print(f"\n DAG metrics:")
print(f" Avg acyclic rate: {analysis['dag_metrics']['avg_acyclic_rate']:.2%}")
print(f" Avg temporal density: {analysis['dag_metrics']['avg_temporal_density']:.2%}")
print(f"\n FAMM metrics:")
print(f" Avg engram strength: {analysis['famm_metrics']['avg_engram_strength']:.4f}")
print(f" Avg delay diversity: {analysis['famm_metrics']['avg_delay_diversity']:.2f}")
print(f" Note: {analysis['note']}")
print("\n" + "=" * 70)
print(" KEY FINDINGS")
print("=" * 70)
print("\n1. Regular graph construction provides more structured graphs")
print(" - All vertices have same degree")
print(" - Better chance of containing required cycles")
print("\n1. DAG structure provides temporal ordering:")
print(" - Topological layers encode temporal sequence")
print(" - Acyclic constraint ensures no directed cycles")
print(" - Temporal density measures cross-layer connectivity")
print("\n2. Exhaustive cycle detection finds more cycles")
print(" - DFS explores all possible paths")
print(" - Cycle diversity indicates richness")
print("\n2. FAMM delay lines enable hippocampal temporal sequencing:")
print(" - Delay matrices capture multi-step temporal flow")
print(" - Engram consolidation integrates weighted delays")
print(" - Temporal integration measures cross-delay coherence")
print("\n3. 4-primitive framework provides structural insight:")
print(" - Spectral: eigenvalue structure")
@ -238,6 +398,11 @@ def main():
print(" - Shear: regularity metrics")
print(" - Packet: cycle encoding")
print("\n4. DAG + FAMM enhance investigation:")
print(" - Temporal structure may influence cycle formation")
print(" - Delay lines capture temporal dynamics")
print(" - Engram strength measures temporal integration")
# Save results
output_data = {
"test_info": {
@ -246,8 +411,9 @@ def main():
"min_degree": 3,
"samples_per_n": 5,
"total_tests": len(n_values) * 5,
"graph_construction": "Regular graph",
"cycle_detection": "Exhaustive DFS"
"graph_construction": "DAG (Directed Acyclic Graph)",
"temporal_sequencing": "FAMM delay lines",
"cycle_detection": "Exhaustive DFS on symmetrized graph"
},
"results": results,
"investigation_analysis": analysis,
@ -257,9 +423,19 @@ def main():
"shear": "Regularity metrics indicate graph uniformity",
"packet": "Cycle encoding captures power-of-two witness"
},
"dag_insights": {
"topological_ordering": "Temporal layers encode sequence",
"acyclic_constraint": "No directed cycles",
"temporal_density": "Cross-layer connectivity measure"
},
"famm_insights": {
"delay_lines": "Multi-step temporal flow capture",
"engram_consolidation": "Weighted delay integration",
"temporal_integration": "Cross-delay coherence measure"
},
"validation": {
"status": "INVESTIGATION_COMPLETE",
"insight": "Refined investigation using regular graph construction and exhaustive cycle detection. Previous random graph method may not have found power-of-two cycles due to graph structure. Regular graphs provide better testbed for conjecture."
"insight": "Refined investigation using DAG structure + FAMM delay lines for temporal sequencing. DAG provides topological ordering. FAMM captures hippocampal temporal dynamics. Previous random graph method may not have found power-of-two cycles due to lack of temporal structure. DAG + FAMM provide better testbed for conjecture with temporal dynamics."
}
}